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ABSTRACT 

This final report documents the design, implementation, and testing of a 
digital computer program, and the specification of the necessary interface 
hardware for the control of an early-launch, laser tracking system. The con- 
trol computer is to accept angular pointing error information from a star 
tracker tube, phase error information from a laser ranging system, and mode 
control from a reacquisition system. During the tracking mode the computer 
will provide drive signals for the tracking telescope in order to null angular 
pointing errors, and will provide phase shift commands to the digital oscilla- 
tor in the ranging system in order to null the phase and, hence, range error. 
In the reacquisition mode the control computer will output commands to the 
telescope drive and the range system oscillator, based on extrapolated target 
motion, so as to minimize initial errors when the track mode is reentered 
after target reacquisition. In addition to its control functions, the pro- 
gram computes target position, velocity, and acceleration in a launch-site- 
centered coordinate system. Provision is also made for interfacing with 
printing, plotting, and magnetic tape units in order to display and store 
the tracking information. 
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SECTION 1 

INTRODUCTION AND SUMMARY 

This final report documents the design, implementation, and testing of a 
digital computer program, and the specification of the necessary interface 
hardware for the control of an early-launch, laser tracking system. The con- 
trol computer is to accept angular pointing error information from a star 
tracker tube, phase error information from a laser ranging system, and mode 
control from a reacquisition system. During the tracking mode the computer 
will provide drive signals for the tracking telescope in order to null angular 
pointing errors, and will provide phase shift commands to the digital oscilla- 
tor in the ranging system in order to null the phase and, hence, range error. 
In the reacquisition mode the control computer will output commands to the 
telescope drive and the range system oscillator, based on extrapolated target 
motion, so as to minimize initial errors when the track mode is reentered 
after target reacquisition. In addition to its control functions, the pro- 
gram computes target position, velocity, and acceleration in a launch-site- 
centered coordinate system. Provision is also made for interfacing with 
printing, plotting, and magnetic tape units in order to display and store the 
tracking information. 

The main body of the report is organized, as nearly as possible, accord- 
ing to what is considered a typical procedure for the design of a digital 
computer-based, real-time, control system. In detail, the following are con- 
sidered to be the major distinguishable stages of a thorough design process. 

1. An input-output performance-oriented specification leading to 
a functional design in engineering terms (e.g., the setting of 
servo loop bandwidths to achieve desirable position and velocity 
tracking errors) . 

2. The translation of the design into an algorithmic form (or the 
devising of an algorithm to realize the desired operations) suit- 
able to implementation in a digital system. 

3. Coding of the algorithm in a programming language. 

4. Simulation and test of the individual subsystems in order to 
assess performance and to effect any necessary modifications. 
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5. Integration of all subsystems Into the overall system which 
is then computer tested to assure satisfactory data transfer 
and real-time performance. 

In the final program there are four major subroutines, for each of which 
the above five stages are described in Sections 2 through 4, with supplemen- 
tary appendices. Stages 1 and 2 are covered in Section 2. Stages 3 and 4 
are described in Section 3; and stage 5 is the subject of Section 4, Descrip- 
tive listings of all subroutines will be found in Appendix A. 
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SECTION 2 
SYSTEM DESIGN 

2 . 1 PHILOSOPHY OF SYSTEM OPERATION 

The philosophy to be dispensed in this section will be limited to the 
way in which the basic angular control system will function in the three al- 
lowable system modes. These three modes relate to the acquisition, tracking, 
and reacquisition functions of the overall system. It is best to begin with 
a discussion of the track mode, which is the steady-state mode for the system 
program. 

A functional block diagram is presented in Figure 2-1. The salient fea- 
ture of this diagram is the way in which the error signals are obtained for 
the mount angle servo. Because of the presence of the target dynamics esti- 
mator it might be thought desirable to Interpose the estimator between the 
star tracker output (the raw angular tracking error data) and the position 
error input point of the angle servo, in order to smooth the noise on the star 
tracker output signal. However, the time constant of the estimator would in- 
troduce an undesirable lag. In the absence of any filtering by the estimator, 
the proportional plus integral compensation block in Figure 2-2 must include 
any noise filtering which has to be done. Of course this noise filtering de- 
creases the bandwidth of the position loop, usually below a value considered 
satisfactory. The classical tradeoff between noise attenuation and a wide 
bandwidth tracking loop can be avoided by using the estimator to provide a 
velocity feed-forward signal. By this means good dynamic tracking performance 
can be obtained while position error measurement noise is substantially re- 
duced. 

When the return signal from the target is lost, a signal from the reac- 
quisition hardware will cause the tracking system to enter the reacquire mode. 
The basic feature of operation in this mode is the provision of all command sig- 
nals to the angle servos by the target dynamics extrapolator. The philosophy 
here is that loss of signal does not indicate a sudden change in target motion, 
causing a loss of track, but only an obstruction to transmission of the laser 
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Figure 2-1 . Angular Control System Configuration in Track Mode 
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beam (typically a cloud or vapor). Hence while the reacquisition system is 
searching for a return signal, the telescope will be kept pointed at the ex- 
pected target position as provided by the extrapolator. Provided the loss of 
signal condition does not prevail for too long a time period (less than 1.5 
seconds, say), the angular tracking error may not be greater than 25 to 30 
arc seconds, depending upon the target dynamics at the time of loss of signal. 

Finally, it is necessary to mention system operation in the acquire mode. 
Reference to Figure 2-3 shows that basically the telescope mount is following 
a programmed target position- -which could be a single position, as in calibra- 
tion tests and initial launch vehicle acquisition, or a sequence of positions, 
as in acquiring a satellite. The reason for leaving the estimator in is to 
minimize the velocity and acceleration estimation transient following switch- 
ing to the track mode. 

2.2 DESIGN OF SERVOS FOR CONTROL OF MOUNT ANGULAR MOTION 

The hour angle and declination axes Of the tracking telescope are in- 
dividually controlled by identical servos. Because these high-performance 
servos are implemented in a digital computer, much flexibility in design is 
available. In particular it now becomes feasible to implement a procedure 
for compensating for servo saturation during large signal transients. The 
existence of the target dynamics estimator allows the use of velocity feed- 
forward, permitting, in turn, the smoothing of position error signals without 
loss of dynamic performance. In other respects the design is straightforward. 

The full design details will be found in Appendix B. 


2.3 DESIGN OF RANGING SERVO LOOP 

The range hardware external to the computer provides a three-channel 
output for phase error signals and requires, as input, a phase -shifting com- 
mand signal to control the digital local oscillator. It is the function of 
the software portion of the range servo loop to process the phase error sig- 
nals both to close the loop via the input command to the digital oscillator. 


2-4 



F-7 180-1 



2-5 




F-7180-1 


and to accumulate phase error inputs so as to provide the instantaneous range 
value. Hence, the design problems are those of a standard servo (compensation 
to assure stability and reasonable tracking errors) with the exception of 
choosing a word size for the oscillator phase-shift signal. 

As detailed in Appendix C, proportional plus integral compensation was 
chosen to eliminate velocity error and limit the acceleration error to less 
than one centimeter at a maximum range acceleration of 10 m/sec . Choosing 
the feedback word (phase shift command to digital oscillator) to contain 10 
bits (plus sign) keeps the range error to less than 0.5 cm for a maximum range 
rate of 360 m/sec. An internal repetition rate of 256/sec is sufficient to 
implement the design. 
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2.4 PROCEDURES FOR COORDINATE CONVERSION AND REFRACTION CORRECTION 

A detailed description of the coordinate systems and coordinate conver- 
sion equations used is given in Appendix D. The refraction correction equa- 
tions are given in paragraph 3.5. A brief verbal description of the pro- 
cedures will be given here. 

Initial target position information is obtained in an equatorial coordi- 
nate system. The first step is to convert to coordinates relative to a rec- 
tangular topocentric system, at the tracker site. The standard Apollo topo- 
centric system is used with axes directed to the south, to the east, and 
along the local vertical. At this point a second, but partial, conversion 
is performed, solely to allow refraction corrections to be made; the elevation 
angle (above the local horizon) is computed from the rectangular topocentric 
coordinates. This elevation angle is corrected for atmospheric refraction, 
and the rectangular topocentric coordinates are then corrected through a back 
conversion. The next coordinate system change is from rectangular topocentric 
at the tracker to rectangular topocentric at the launch site. Finally, the 
rectangular topocentric coordinates at the launch site are converted to spher- 
ical coordinates. 

It should be apparent that the final data accuracy (for target position 
in spherical topocentric coordinates at the launch site) depends critically 
upon the accuracy of the transformation from the tracker site to the launch 
site. The quantities influencing this accuracy are five in number: they 

are the absolute geodetic latitude of the launch site, the height of the origin 
of coordinates above the reference ellipsoid for both tracker and launch sites, 
and the differences in both geodetic latitude and longitude between the two 
sites. 
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2 . 5 PRINCIPLES OF OPERATION OF TARGET DYNAMICS ESTIMATOR AND EXT'RAPOLATOR 

A detailed explanation of the theory of the estimator algorithm Is given 
In Appendix E. In this section a nonmathematlcal and heuristic discussion Is 
given. In an attempt to set forth the basic operating principles. Before pro- 
ceeding to that discussion, a few remarks should be made about other estimators 
which were studied and the reasons why they were found unsuitable. The most 
general linear least-squares estimator Is the Kalman-Bucy filter which requires 
a full description of the target dynamics and expected noise. For a Saturn V 
launch vehicle the equations of motion are far too complex to allow a real-time 
estimator to be realized on a medium-scale computer. In addition, the straight 
Kalman-Bucy filter has a rather long time constant (due, mainly, to the assumed 
completeness of the signal and noise dynamical models) which would not allow 
accurate enough tracking during a loss-of-englne Incident. A short- time -con- 
stant, or fast responding, estimator must be designed If good tracking data 
Is to be realized following loss of one of the five first-stage engines. Ac- 
cordingly, a modified Kalman-Bucy filter was designed and Implemented; the 
modification was achieved by recycling the error covariance matrix periodically 
In order to Increase the apparent gain of the estimator and so allow any change 
In the data trend to Influence the estimator output. Unfortunately the dis- 
continuity In the velocity signal at the Instant of recycling resulted In too 
high a noise level to provide a useful velocity feed-forward signal for the 
angular control system. 

To achieve fast response and to reduce periodic transients In the esti- 
mator output. It was decided to perform a polynomial fit to short data spans 
(nominally, 90 data points with a 1/64 second spacing for a total period of 
1.4 seconds) and to overlap these fitting polynomials. In choosing the degree 
of the fitting polynomials, consideration must be given to the nature of the 
acceleration estimate. Thus a quadratic would yield a piecewise constant ac- 
celeration, any change occurring only when the output Is derived from the 
next -In-sequence polynomial. From a consideration of the nominal tracking 
data (refer to paragraph 3.2) a quadratic Is eliminated so that a cubic poly- 
nomial must be used. The question still remains as to how many overlapping 
polynomials should be used; two constltue the minimum number with no clear 
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limit to the maximum. An answer to this question must await an elucidation 
of the intended purpose of each of the overlapping polynomials. 

To fit a polynomial to a large number of measured data points is but one 
aspect of least-squares parameter estimation. This parameter estimation may 
be accomplished in two ways in the time domain: by waiting until all measure- 

ments have been made and then estimating the parameters via the classical 
Gauss-Markov formula; or by updating the parameter estimate as each new mea- 
surement is obtained, using the more modern recursive equivalent of the Gauss- 
Markov formula. With obvious justification, the former procedure will hence- 
forth be labelled ba t ch pr oce s s ing , and the latter procedure either real-time 
or sequential estimation . To aid further discussion an additional piece of 
terminology must be introduced: because of the use of overlapping fitting 

polynomials, the complete estimator will actually consist of as many subiinits 
as there are polynomials, each subunit being devoted to estimating the param- 
eters,' or coefficients, of one of the polynomials. Because the overall esti- 
mator is to operate in real time, it is necessary for at least one of the sub- 
units to be performing sequential estimation; the remainder of the subunits 
can be doing batch processing. It is at this point that a decision can be 
made on the number of subunits or polynomials which must be used. The key 
question to be answered is how good must the estimate of the polynomial coef- 
ficients be before Initiation of the sequential estimation phase? The reason 
for asking this question is that the current estimates of position, velocity, 
and acceleration will be derived from that subunit performing sequential esti- 
mation. Hence it is necessary that a reasonable number of data points will 
have been batch processed by the subunit before the beginning of the sequen- 
tial estimation phase. With only two polynomials, or subunits, sequential 
estimation begins after half the data span has been batch processed; experi- 
ence has indicated that the resulting sequential estimation produces large 
errors until more data points have been incorporated. Consequently, three or 
more polynomials must be used. 

The final estimator structure is composed of three subunits; the use of 
four subunits was excluded as requiring too much computation time. The choice 
of 90 samples (with 1/64 second spacing) results from the necessity that the 
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data span be divisible into three equal segments; 60 samples proved inadequate 
from an estimation error standpoint, and 120 samples would result in a time 
constant approaching 1.0 second- -also deemed undesirable because of the re- 
sults of the following analysis. 

The time constant of the estimator will be taken as one half of the time 
duration of a subunit's processing cycle: for the 120 sample cycle the time 

constant would become 60/64 = 0.938 second; for the 90 sample cycle it would 
become 45/64 = 0.703 second. Using the nominal tracking data (see paragraph 

3.2) the peak accelerations for range, hour angle, and declination angle are, 

4 2 2 2 

respectively, 10 cm/sec , 457 arc sec/sec , and 230 arc sec/sec . If one 

of the five first-stage engines were to be shut down, there would be a nomi- 
nal 20 percent decrease in acceleration. As a worst case situation, it will 
be assumed that this change in acceleration does not affect the estimator out- 
put until one time constant after its occurrence. In this case the use of the 
120-sample cycle would result in net errors of 435 cm, 39.8 arc seconds, and 
20 arc seconds, respectively, for range, hour angle, and declination. 

Using the 90 sample cycle gives range, hour angle, and declination errors 
of 247 cm, 22.5 arc seconds, and 11.3 arc seconds, respectively. Thus use of 
120 samples instead of 90 results in an almost twofold increase in tracking 
error in the event of engine malfunction, without a commensurate decrease in 
estimation error. 

A timing diagram is shown in Figure 2-4; it is not intended to be indica- 
tive of the actual estimator program structure, only of its operating prin- 
ciples. In the actual implementation of the triply overlapping polynomial 
fit, only one set of estimator equations is programmed; this set of equations 
is time -shared by the three subunits to update the three parameter sets (see 
paragraph 3.6). 

A detailed presentation of the mathematics underlying the estimator de- 
sign is given in Appendix E. 
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2.6 SPECIFICATION OF SYSTEM HARDWARE 

In the initial stages of this contract it was thought that a small-scale 
digital processor using fixed point arithmetic and a 16- or 18-bit word length 
would be adequate for handling the computational load. However, after the 
coordinate transformations were developed it became apparent that floating 
point arithmetic was extremely desirable both to avoid the awkward scaling 
problems and the use of double precision arithmetic in certain calculations. 
Accordingly, attention focussed on 32- or 36-bit machines that were supplied 
with floating point hardware. 

Both the SIGMA 5 (Scientific Data Systems) and the PDP-10 (Digital Equip- 
ment Corporation) are considered adequate to handle the real-time computa- 
tional and formatting requirements imposed by the system design. Both machines 
can be supplied with the complete line of interface and peripheral equipment 
necessary to perform all A/D and D/A conversions, data sampling functions, and 
control of printer, plotter, and magnetic tape units. The manufacturer- 
supplied system configurations and equipment specifications for the PDP-10 and 
SIGMA 5 are given in Figures 2-5 and 2-6. 

The input/output data signal characteristics are specified below: 

Input Signal Characteristics 

A. Digital Signals 

1. 2 angle encoder signals^ — 22 bits each at 64/second 

2. Clock signal 

B. Analog Signals To Be Digitized 

1. 2 Star tracker signals--14 bits each, including sign, at 
64/ second 

2. 3 range phase discriminator signals--14 bits each, including 
sign at 2 56 /second 

3. 2 tachometer signals--14 bits each, including sign, at 64/second 
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ITEM 

QUANTITY 

DESCRIPTION 

UNIT PRICE 

AMOUNT 

1 

1 

PDP-lO/10 System 

Includes; KAIO Arithmetic Processor 
(with 300 char /sec Paper 
Tape Reader, 50 char /sec 
Paper Tape Punch, Console 
Teleprinter, lOP with 7 
levels of priority 
interrupt 

$113,000 

$113,000 



MAIOA 8,192 Word Core 
Memory (1 u sec cycle 
time) 



2 

1 

KEIO Extended Order Code - includes 
hardware for floating point and 
variable byte manipulation instruc- 
tions) 

12,000 

12, 000 

3 

1 

TMIO Magnetic Tape Control for DEC 
TU20 Tape Transports (controls up to 
8 transports) 

18, 000 

18,000 

4 

' 1 

' TU20 Magnetic Tape Transport (200, 
556,800 bpi at 45 "/sec or 36KC max. 
' character rate) 

12,000 

' 

12,000 


TERMS: NET 30 DAYS, F.O.B. MAYNARD, MASSACHUSETTS 
i DELIVERY SCHEDULE + 


•mis QUOTATION IS BASED UPON ACCEPTANCE WITHIN 

• DAYS FROM THE DATE HEREOF. AND IS 
SUBJECT TO CREDIT APPROVAL AND TO THE TERMS 
AND CONDITIONS HEREON AND ON THE REVERSE SIDE. 
ANY CONTRACT RESULTING FROM THIS QUOTATION 
MUST BE SIGNED IN MAYNARD, MASSAOHUSETTS, 
BY A DULY AUTHORIZED REPRESENTATIVE OF DIGITAL 
EQUIPMENT CORPORATION. 

DIGITAL EQUIPMENT CORPORATION 
.( ■ ' 


By 


c. 




Edward A. Kramer - Applications Engineer 
CUSTOMER 

Figure 2-5. Equipment Specifications for PDP-10. (Sheet 1 of 2) 
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QUOTATION NO- CONTINUED 

G-113-113-298 

PAGE 2 OF 2 PAGES 

PLEASE REFER TO THIS QUOTATION NO. IN 
ALL CORRESPONDENCE AND ORDERS. 

DATE Noveirtoer 20, 1967 


ITEM 

QUANTITY 

DESCRIPTION 

UNIT PRICE 

AMOUNT 

5 

1 

KMIO Fast Registers (16 general 
purpose register accumulators, 150 
n sec access time) 

$ 9,000 

$ 9,000 

6 

1 

AA05 Multiple Digital-to-Analog 
Converter Control (controls up to 
64 DAC's) 

4,500 

4 , 500 

7 

4 

A608 10-bit, 10 u sec D/A Converter 
Module 

350 

1,400 

8 

1 

RTC Real Time Glock - (Selectable 
frequency, generates interrupt on 
overflow) 

900 

900 

9 

1 

LDU-lO 36 bit output unit (type 
LDC) 

3,300 

3,300 

i „J 

1 

IMU-10 3-36 bit word input unit 

3 , 500 

3,500 

11 

1 

GP I/O General Purpose I/O hardware 
interface 

3,000 

3,000 

12 

1 

AD-14 Special A/D Converter (13 bits 
+ sign, includes sample & hold, mul- 
tiplexer control, 200 KC max. con- 
version rate, 8 multiplexer switches) 

17,350 

17,350 



TOTAL 


$197,950 


CUSTOMER 

Figure 2-5. Equipment Specifications for PDP-10 (Sheet 2 of 2) 
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COMPUTER: 
Qty Model 

SDS SIGMA 5 SYSTEM 

Description 

Purchase Price 

1 

8201 

Sigma 5 CPU with integral I/O Processor^ Two 
real time c locks, control panel and power supplies. 

$ 50/000.00 

1 

8214 

Memory Protect 

4,000.00 

1 

8218 

Floating Point Hardware 

25,000.00 

1 

8251 

Memory Module: 4096 Words (32 bits + parity) 

33,000.00 

1 

8252 

Memory Increment: 4096 Words (32 bits + parity) 

17,500.00 

1 8270 

PERIPHERALS: 

External Interface Feature 

Computer Price: 

2,000.00 

$131,500.00 

Qty 

Model 

Description 

Purchase Price 

I 

7010 

Console keyboard/printer and controller (KSR35 
teletype). 

$ 6,000.00 

1 

7060 

Paper tape reader 300 cps; paper tape punch 60 
cps; spooler; cabinet and controller. 

12,000.00 

1 

7361 

Magnetic tape controller: 1 - 8 units 

6,000.00 

1 

7362 

Magnetic tape transports: 7 channels; 556 bpi; 
37-1/2 i ps . 

19,000.00 

1 

7922 

Analog and Digital Adapter 

4,000.00 

1 

7930 

Digital I/O Adapter 

3,600.00 

6 

7950 

Eight (8) Stored digital outputs. 

708.00 

1 

7953 

Eight (8) Pulsed Digital Outputs 

Peripheral Price: 

75.00 
$ 51,383.00 


Figure 2-6* System Configuration and Equipment Specifications for 
Sigma 6 (Sheet 1 of 3) 
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ANALOG SECTION: 
Qty Model 

Description 

Purchase Price 

1 

8920 

Cabinet 

$ 1,000.00 

1 

ZX14 

Blower Assembly 

180.00 

1 

PTIO 

Power Supply 

280.00 

1 

PX12 

Power Supply 

280.00 

1 

PXIO 

Power Supply 

370.00 

2 

ETlO-15 

Cables 

140.00 

1 

ET21 

Cable 

250.00 

1 

ET20 

Cable 

250.00 

1 

EZ28 

Cable 

20.00 

1 

EZ 

Coble 

250.00 

1 

DA35-15 

D/A Controller 

1,800.00 

*4 

DX16 

D/A Converter: 1 1 bits + sign; ±10V 

2,800.00 

1 

AD35 

A/t) Converter: 14 bits + sign 

6,040.00 

1 

MU55 

Multiplexer for up to 32 channels 

2,450.00 

2 

SX40-2 

Four Adjustable analog switches 

360.00 

1 

SSIO 

Sample and Hold Controller 

1,000.00 

7 

HX35 

Sample and Hold Amplifier 

3,500.00 



Analog Section: 

$ 20,970.00 



SYSTEM TOTAL: 

$203,853.00 


*Two of the D/A converters are for a pen recorder. If a SDS model 7530 Graph Plotter (Calcomp) is 
substituted two of the DX16*s can be eliminated. The price of the model 7530 ( Plotter and Controller) 
is $13,000.00. 

Figure 2-6- System Configuration and Equipment Specifications for 
Sigma 5 (Sheet 2 of 3) 
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Output Signal Characteristics 

A. Digital Signals 

1. 3 range servo signals — 10 bits each at 256 /second 

2. 3 strobe signals to range servo--l bit each at 256/second 

3. Digital signals, as required to printer and magnetic tape 

B. Digital Signals To Be Converted 

1. 2 angle servo drive signals^ — 10 bits each at 64/second 

2. 2 signals to X-Y plotter- -10 bits each at 64/ second 

As noted from the above listing, seven analog input signals must be 

sampled and digitized 64 times/second along with the sampled angle encoder 

values. Accuracy considerations, which required an ultimate angular resolu- 
-21 

tion of 2 of a full circle, demanded that a tolerance be established for 
the maximum time skew of the sampled quantities. Assuming a worst case value 
for the angle rates (4 rad/sec), it was found that a time skew of 100 jusec 
for the sampled quantities was a sufficient margin to guard against inaccura- 
cies being introduced at the input stage. In consultation with the computer 
manufacturer's representatives, it was decided to specify sample and hold 
circuitry in conjunction with the A/D converter. 
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SECTION 3 

SOFTWARE DESIGN, IMPLEMENTATION, AND SIMULATION 


3 . 1 OBJECTIVES 

The purpose of the system simulation was to exercise as many aspects of 
the system performance as possible in order to test the accuracy of the algo- 
rithms and the consistency of interfacing between the various subroutines. 
System simulation was only carried out for the tracking mode, for reasons to 
be stated. The flow chart in Figure 3-1 shows system operation, as simulated, 
blocked out in terms of subsystem functions. In some cases these blocks cor- 
respond directly to subroutines used in the simulation, such blocks being 
noted with an asterisk. 

The basis for generating the polynomials to give range, hour angle, and 
declination angle data is explained in paragraph 3.2. In order to obtain 
realistic data for exercising the design algorithms it was necessary to simu- 
late the actual range and angle servo hardware equipment. In this manner it 
was possible to obtain range channel, star tracker, and tachometer Inputs con- 
sistent with the actual trajectory data and the response of the equipment to 
these inputs; thus the inclusion of the blocks "RSERVO" and "ASERVO" as shown 
in Figure 3-1. 

A few words of explanation are in order concerning the blocks enclosed 
in dotted lines. These blocks develop the timing sequence for the operation. 
The outer loop ((^) is incremented every second, the next loop ((^) is in- 
cremented 64 times every second, and the innermost loop ((£)) is Incremented 
256 times every second. Thus the range servo loop is made to operate at four 
times the rate of the angle servo loop, the larger rate being required for 
achieving sufficient ranging accuracy. The entire system was simulated over 
a period of approximately 70 seconds, the length of time over which supplied 
trajectory data is available. 

Each of the blocks labeled with an asterisk was tested individually be- 
fore being combined into the overall system slinulation. In the following 
paragraphs the structure of each of the subsystem subroutines is detailed. 
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NOTE I: THIS BLOCK IS EQUIVALENT TO PREDETERMINED 

POINTING OF ANGLE SERVOS AND STEPPING 
UP OF RANGE SERVO TO INITIAL VALUES 

NOTE 2: THE BLOCKS ENCLOSED WITHIN THE DOTTED 

LINES DEVELOP THE TIMING SEQUENaS FOR 
THE RANGE AND ANGLE SERVO LOOPS 



♦RSERVO 

SUBROUTINE 


*LT05 

SUBROUTINE 


*ASERVO 

SUBROUTINE 


♦LT03 

SUBROUTINE 


*LT04 

SUBROUTINE 


Figure 3-1 . System Simulation Flow Chart 
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The real-time projections for system performance are discussed in paragraph 
4.1. 

The reasoning for restricting the overall simulation testing to the 
tracking mode is as follows. In the track mode all basic subsystem functions 
are exercised; only the logical steps necessary to switch between modes are 
absent. If the various algorithms did not perform properly in the track mode 
then, in any case, they would be unsatisfactory. The only feature of system 
behavior to be tested in the reacquisition mode is the target state estimator. 
This estimator (see paragraph 3.6) has been designed on the basis of MSFC- 
supplied reference trajectory data and a pessimistic noise level. Thus it 
must be assumed that under normal signal dropout considerations, the estimator 
can provide usable output during the reacquisition phase based on extrapola- 
tion of the current best estimates of velocity and acceleration. Abnormal 
signal dropout conditions, that is, loss of signal for periods exceeding one 
or two seconds, cannot be accounted for in the design. In the acquisition 
mode the servos are to respond to pointing commands generated from pre-assigned 
coordinate information. 

For the angle servos the pointing commands can be implemented directly, 
with no "stepping" increment necessary. For the range servos, however, it 
is necessary to use increments that are no greater than one half the coarsest 
range channel resolution to "step" the range servo up to its initial value 
and thereby avoid saturation. 
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3.2 COMPUTATION OF N(MINAL TRACKING DATA 

In order to determine the tracking data characteristies it was necessary 
to write a special computer program to 

1. Convert MSFC-supplied trajectory data to tracking data for a 
particular launch site--tracker site geometry, 

2. Derive accelerations for the tracker-measured quantities, and 

3. Do a least-squares fit of polynomials to the results so that 
simulation program inputs could he obtained. 

The MSFC-supplied trajectory data were position and velocity, at one second 
intervals, in an inertial coordinate system which was earth-fixed and launch- 
derived. That is, the trajectory was known in a geocentric inertial coordin- 
ate system which was aligned with the launch azimuth and the geodetic vertical 
at the launch site. The tracker site was chosen as that one of the allowable 
tracker sites which was closest to the launch site; in this way, the maximum 
velocities and accelerations of the tracking data would be obtained. The re- 
sults are tabulated in Figure 3-2; the acceleration information, which was 
not in the original trajectory data, was obtained by differencing the velocity 
data. The coefficients of sixth-degree, least-squares fitted poljmomials 
are in Tables 3-1 and 3-2. Significant excerpts from the data, pertaining to 
maximum values, are given in Table 3-3. Note that the form of the fitted 
polynomial is 

2 3 4 5 6 

Cq + ^ 3 ^ ^ 6 ^ 

Data on the coordinate conversions performed are in Appendix D. 


3.3 MOUNT ANGLE SERVO SUBROUTINE 

A full discussion of the software, and simulation results, is given 
in Appendix B. 
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TfMF 



HOUR ANRLF 

OECUNATION 


acCFlE5ATtriy 


accflepatiom 

ACC5L58ATinw 

ft 

.3 

7.6554045gr. 

ft 

03 

0 

-1,442395905-02 

0 

8,T40flll86F.ft3 

v.o 

l.«506n9?F- 

02 

-3.441755295-02 

2.084250ft25-fl? 

2.0 

3»?6047?5iF- 

02 

-5.029764875-02 

3.03432389F-02 

3.0 

4.569ft6785r- 

02 

-5.17820124F-02 

3.095599325*02 

4.0 

6.7l«4286Qr- 

0? 

-5.336224865-02 

3*l630lO?7F-ft2 

5.0 

9.69375740F- 

02 

-5.510U083F-02 

3*lR0423l8F.ft? 

6.0 

l. 358621? 1 r- 

01 

-5.7n337o5oe-02 

3.207380775-0? 

7.0 

l.«44^8975^- 

01 

-5e5i53l69lE-02 

3.221680925-0? 

6.0 

2e42356375F- 

01 

-6.145352565-02 

1.222670395-O? 

6. ft 

3,109o866nF- 

■01 

-6.39789362E-02 

3e208?48695-ft? 

10.0 

3.9036?04?F- 

■01 

-6.67731762E-02 

3.180456825-0? 

n.n 

4e«054395lF- 

■01 

-6.98267400E-02 

3.137406515-02 

l?.o 

5,fl5877702F- 

-01 

-7.321855645-02 

3.059187V3F-02 

13.0 

6,580^6«35F. 

■01 

-7.66524770E-02 

2,979376035-0? 

14.0 

fi.10744902F- 

►01 

-7,996253652-02 

2.918796815-02 

15.0 

9,776ftl57qF- 

■01 

-8.32702323E-02 

2.862205115-02 

16.0 

1.ft55flO«63F 

00 

-8.A7458658E-02 

2.788910065-02 

17.0 

l,l97?436ir 

no 

-9,043971845-02 

2,689730575-02 

l«.o 

1 . 352051 89F 

00 

-9.4313U855-02 

2.555601915-02 

15.0 

1.522fl5o55f 

on 

-9.837216355-02 

2.376378395-02 

?o.o 

1 .7lOl564or 

00 

-1.0261V276F-01 

2.148748325-0? 

21.0 

1. 51377331 F 

on 

-1,069499865-01 

1,869319935-02 

22.0 

2.V3203883F 

00 

-1,112686575-01 

1.535909485-0? 

23,0 

2,364ftll37F 

on 

-1.153568035-01 

1,150951325-0? 

24.0 

2,610374«PF 

00 

-1,191504805-01 

7,081775075.03 

25.0 

2.86642484F 

Oft 

.1.22408573F-01 

2.188369755-03 

26.0 

3.13541^77F 

no 

-1.249703635-01 

-3.239215295-03 

27.0 
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Oft 
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2*^.0 
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no 
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25.0 
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no 
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30,0 
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no 
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-2.809829875-0? 

31,0 
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00 

-1,179531095-01 

-3.661113475-0? 

32.0 

4.85872447F 

00 

-1 • 1 ft896483F-01 

-4.08557680F-O? 

33.0 

5.14962mF 

00 

-1.015352365-01 

-4,677748525-0? 

34.0 

5.43534634F 

00 

-8.586439965-02 

-5.211224875-0? 

35,0 

5,709V470«;f 

00 

-7,615261275-02 

-5,653700035-0? 

36.0 

5,079«6105F 

00 

-6.045293305-02 

-6,014875465.02 

37.0 

6.23848494F 

no 

-4.335280435-02 

-6,261488165-0? 

3«.0 

6,48670825F 

00 

-2.532296525-02 

-6,390884555-0? 

35.0 

6.7200447ftF 

00 

-7.067098635-03 

-6,396652725-0? 

40,0 

6,5451 4767F 

no 

1.099473505-02 

-6,294108285-0? 

41.0 

7.l59ft7834F 

00 

2.818224475-02 

-6,086030865-02 

42.0 

7,35752793? 

00 

4.372883115-02 

-5.778929545-0? 

43.0 

7.S5008455F 

no 

5.774602565-02 

-5,403064405-n? 

44.0 

7.728ft8l49F 

00 

6.945258565-02 

.4,96?798575.o? 

45.0 

7.R9861884F 

00 

7.9O881434F-02 

-4,484330205-02 

46,0 

8.o57ft639lF 

00 

8,628125955-02 

-3,979067495-0? 

47.0 

8.21056754F 

00 

9.162121245-02 

-3,470345655-02 

4«,0 

8.35358991F 

Oft 

9.480440605-02 

-2.965666325-0? 

45.0 

8.45148143F 

00 

9.6?9ai3985-02 

-2,678694065-02 

50,0 

8. 62342251 F 

on 

9,6?6319465-Q2 

-2.018749185-0? 

51.0 

8,760Q85nF 

00 

9.516863575-02 

-1.590543445-0? 

52.0 

8.R7046556F 

00 

9.300522555-02 

-1,199204235-02 

53,0 

8,50741R61F 

00 

9,015199555-02 

-8.462697185-03 

54.0 

9,10245035F 

00 

8.680038115-02 

-5.318174205-03 

55.0 

9.21335196F 

00 

8,300829835-02 

-2.571661705-03 

56.0 

9.3233119PF 

00 

7,910352515-02 

-1,699442075-04 

57.0 

9.431ft2032F 

00 

7,498025735-02 
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5».0 
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60.0 
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00 
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62.0 
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00 
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63.0 
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01 
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01 
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65.0 
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01 
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1,044o9UOF 

01 
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01 
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68.0 
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01 
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01 
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01 
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71.0 
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01 
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72.0 
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01 
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Figure 3-2. Tracking Data as Obtained from Trajectory 
Data (Sheet 1 of 2) 
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Figure 3-2. Tracking Data as Obtained from Trajectory 
Data (Sheet 2 of 2) 


3-6 



F-7 180-1 


TABLE 3-1 

COEFFICIENTS OF POLYNOMIAL FIT TO TRACKING DATA (POSITION) 



Range (meters) 

Hour Angle (degrees) 

Declination (degrees 


3.074285390E 03 

1.174125936E 02 

4. 103893 12 IE 01 

^1 

-2.855444864E 00 

-8.129312544E-02 

2. 12 902 66 3 IE- 01 

s 

6.487582373E-0I 

-2.099399294E-02 

-3.227090683E-02 

s 

-4.946418019E-02 

8.142128362E-04 

3. 69409015 6E- 03 

^4 

1.930160010E-03 

-8.858125812E-05 

-1. 181951605E-04 

s 

-1.928922545E-05 

1.941353416E-06 

1.480725462E-06 


6.734844011E-08 

- 1.22 18085 lOE- 08 

-6.537293629E-09 
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TABLE 3-2 

COEFFICIENTS OF POLYNOMIAL FIT TO TRACKING DATA (VELOCITY) 


Range Rate (m/sec) 

Cq -5.053388666E-02 

2.507330732E-02 

Cg 1. 4376448 15E- 02 

Cg -1. 177956236E-03 

C, I.394900353E-04 

4 

C- -2.590478355E-06 

0 

C„ 1.457511935E-08 

b 


Hour Angle Rate (°/sec) 
4.709042989E-02 
-1.330068257E-01 
1.715535798E-02 
-1.296911800E-03 
3.784461609E-05 
-4.654506566E-07 
2. 060032254E-09 


Declination Rate (°/sec) 
7.323323364E-03 
2.449193662E-02 
1.067491980E-04 
1.081695380E-04 
-7. 557305 137E-06 
1.456644425E-07 
-8.779346109E-10 
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TABLE 3-3 

MAXIMUM VALUES OF TRACKING DATA QUANTITIES 

Va lue Time of Occurrence 


r 

10.21 km 

68 sec 

r 

354.42 m/sec 

68 sec 

r 

10.806 m/sec^ 

68 sec 

^HA 

117.32 deg 

0 sec 

^HA 

-3.258 deg /see 

39 sec 

V 

-0.1269 deg/sec 

28 sec 

d 

DEC 

56.75 deg 

40 see 

®DEC 

0.644 deg/sec 

28 sec 

®DEC 

2 

-0.06397 deg/sec 

39 sec 


3-9 



F-7 180-1 


3.4 RANGE SERVO SUBROUTINE 

The block diagram in Figure 3-3 shows the implementation of the Range 
Servo loop as simulated on the digital computer. The linearized and simplified 
representation of the hardware is simulated on the computer for testing purposes 
only; while the electronic controller represents a part of the software for the 
Laser Tracker. 

To achieve the required accuracy of 0.5 cm in the output range word, a 
cycle time (TP) of 1/256 sec for the electronic controller was found necessary; 
a 9 bit A/D converter was found adequate and a bandwidth of 55 rad/sec for 
the servo loop was selected. Depending on the magnitude of the error signal 
in the closed loop, the range servo will track in either the fine, medium, 
or coarse mode. As seen from Figure 3-3 these channels differ in gain values 
only. A digital signal from the peripheral A/D converter, proportional to 
the error in the loop, serves as the input to the controller. For a phase error 
of 1/2 rad, which corresponds approximately to an error in range of 40 cm, the 
fine channel goes into saturation and the software switches the tracking opera- 
tion from channel 3 to channel 2. 

The range rate and the output range word are calculated numerically. 
Proportional plus integral compensation is included to eliminate velocity error, 
reduce acceleration errors, and assure good stability. It can be shown that 
the acceleration errors in such a system are given by 


R„ = 


R 


w. 


IP 




CP 


where co^p is the lead corner frequency of proportional plus integral compensa- 
tion (16.7 rad/sec) and o)^p is the crossover frequency of the servo loop (55 
rad/sec). The value of acceleration at each iteration is calculated by time 
differencing and averaging of available range rate words and past acceleration 
values. This value of acceleration is then used to compensate for the accelera- 
tion error in the output range word. 
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Figure 3“3 . Block Diagram of Range Servo 
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Figure 3-4. Flow Diagram of the Program for Range Servo (Sheet 1 of 2) 
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NOMENCLATURE 


DESCRIPTION , 


RI 

X(0) 

x(-i) 

•^ 1 ’ h.2’ hs 

^B3 

TPl 

TP 

^4 

El, E2, E3 

IDR 

RDD 

PR 

Step 1 

Step 2 

Step 3 

^ 2 * ^3 


Simulated range input from nominal tracking data 

Value of X at this instant, t 

’ n 

Value of X at previous time, 

Phase detector input to channels 1, 2, 3, respectively 
Phase error in the loop of channels 1, 2, 3, respectively 

Feedback phase signal for channels 1, 2, 3, respectively 

Time for one computation in hardware section corresponding 
to 20 kHz: 5x10"^ sec 


Time for one computer iteration- -corresponding to 256 sec 
or 4x10"^ sec 

Filter cutoff frequency = 160 rad/sec 
Gain constant = 10 v/rad 

Output of filter in channels 1, 2, 3 respectively 
Change in range in steps of 1/8 cm 

Computed acceleration — for compensation of acceleration 
error 

Previous value of RR (i.e., 8 ♦Range) 

8 ^ lio ^ = 


-1 


^^2 1 11 . „ 

^c ^ 8 ^ 100 ^ TP = ^*2 X 10 




Step increment in feedback phase in channels 1, 2, 3, 
respectively 


Figure 3-4. Flow Diagram of ihe Program for Range Servo (Sheet 2 of 2) 
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The phase feedback in the loop is proportional to the range word calculated 
at the last iteration. This is updated by generating a digital word proportional 
to the change in range, which controls the nmnber of pulses to be added to the 
train of clock pulses. This operation has been explained in detail in paragraph 

2.3. 

The range output, and hence the digital word proportional to the change in 
range, is calculated every 1/256 sec, while the hardware and the converting 
equipment operates at 20 kHz. The synchronization between these two is main- 
tained in the software. 

The flow chart for the computer program is shown in Figure 3-4. Figure 
3-5 shows the response of the servo during "worst case" tracking of a Saturn V 
launch when error due to acceleration has not been compensated for; while 
Figure 3-6 shows the same case with the inclusion of compensation for accelera- 
tion error. It can be seen that error in range at the output stays well with- 
in allowable limits of 0.5 cm. 


3.5 COORDINATE CONVERSION AND REFRACTION CORRECTION SUBROUTINE 

A detailed explanation of the required coordinate conversions is contained 
in Appendix D. In the subroutine, the tracking data is converted into rectangu- 
lar topocentric coordinates at the tracker via Eqs. (D-13) and (D-14) . The re- 
sult is the position vector X . The velocities and accelerations are con- 
verted to and by expressions which are directly obtained by differentia- 
tion of Eqs. (D-13) and (D-14). The second step is a conversion to rectangular 
topocentric coordinates at the launch site via Eq. (D-10) which is valid for X^. , 
Xj.j., and Xj.^ because Ar and T^^^ are constant (i.e., the launch and tracker sites 
have a fixed relationship) . The final step is to convert from the rectangular 
topocentric coordinates at the launch site to spherical coordinates via Eq. 

(D-12) . Once again, straightforward differentiation supplies the expressions 
for obtaining the first and second derivatives of the spherical coordinates from 
— t£’ — tf* — tC’ position, velocity, and acceleration, respectively, in 

rectangular topocentric coordinates at the launch site. Hence the conversion 
from tracker data to launch-site-centered data is straightforward. The coded 
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coordinate conversion subroutine was checked out by inputting to it the nomi- 
nal tracking data (paragraph 3.2). 

The atmospheric refraction correction is made at an intermediate step 
of the coordinate conversion by using the measured target coordinates in 

the rectangular topocentric system, to compute the elevation and azimuth angles. 
Thus the result of combining Eqs. (D-13) and (D-14) may be written in the alter- 
nate form 


-sin 


"Si“ 

COS 9 cos 9. 
E A 

= r 

^2 

cos 9_ sin 9 

E A-' 


-S 3 - 


(3-1) 


where 9 and 9 are elevation and azimuth in a tracker -centered topocentric sys- 

iL A 

tern. Only the elevation angle is corrected (range correction is negligible). 
Letting the correction be the corrected coordinates are 




+ AS^ 

+ AS2 
+ ASg 


where 


ASj^ = -A0g cos 
AS 2 = A0g sin 
AS3= 



(3-2) 


(3-3) 


This is a first-order correction valid for 


A9 


E 


-4 

<2x10 radians , a bound 


which is not exceeded, 
yield 


Empirical fits to standard atmospheric refraction data 
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r • B • C • 


B . r 

1 + e 


cos 
• sin 


(3-4) 


where 


B = 1/7.61 (km"b 
C = 292 X 10‘® 
r = range in kilometers 
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3.6 ESTIMATOR-EXTRAPOLATOR SUBROUTINE 

Implementation of the triply overlapping polynomial fitting described 
in paragraph 2.5 requires the devising of an algorithm whereby a single set 
of estimator equations may be time-shared among the three subunits. This is 
done so as to eliminate triplication of the estimator equations and thus con- 
serve memory space; an additional advantage is the construction of an algo- 
rithm which could be used for any number of subunits if desired. In addition, 
a strategy must be developed for the extrapolator operation during the reac- 
quisition phase. These matters are discussed in paragraph 3.6.1 

Using the nominal tracking data (see paragraph 3.2) and a noise genera- 
tor, the estimator operation was simulated. A discussion of the simulation 
and the results is found in paragraph 3.6.2. 

3.6.1 Software Description 

The equations used during the batch and sequential estimation phases 
are developed in Appendix E. It should be noted that the computation during 
the batch processing phase amounts to the multiplication of a measurement 
vector by a matrix. In the computer program this matrix -vector multiplica- 
tion is done piecewise by multiplying one column vector of the matrix and 
the corresponding element of the measurement vector, as this element of the 
measurement vector becomes available. It is possible to spread the batch 
processing computation over the entire batch processing interval since the 
parameter vector estimate generated is not used until the end of the batch 
processing interval. The computation during the sequential processing phase 
is similarly straightforward, consisting of the evaluation of a cubic poly- 
nomial to get the predicted measurement, and the modification of the param- 
eter vector estimate by adding the current gain vector multiplied by the dif- 
ference between the actual and predicted measurements. Hence the software 
for the actual mathematics of the estimator is simple and straightforward. 

The bulk of the programming is due to time-sharing of the estimator mathemat- 
ics among the three subunits, using the same estimator for range, hour angle, 
and declination, and implementing an extrapolation function for use during the 
system reacquisition mode. The flow diagrams are given in Figures 3-7 and 
3-8; a detailed FORTRAN listing is in Appendix A. 


3-19 



F-7 180-1 



Figure 3-7. Functional Flow Diagram for Estimator — Extrapolator Subroutine 
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INITIALIZATION: 
NSAM = 90 
KSA = 3 
KSB = I 


ENTER 


CALL LT03 (INDEX, RANGE, 
HOUR, DECLIN, IFLAG) 


IK = INDEX - NSAM* (INDEX/NSAM) 


IK = NSAM 


IK - 0?^ 


JIK = IK - (NSAM/3)*(IK/(NSAM/3)) 


JIK = NSAM/3U _<^K = 0^ 

1 I3 


'INDEX < NSAM? 


YES NO 

< NSAM/3^?>— - 


JIK = 1?^ 


IK < 2*NSAM/3? 
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Figure 3-8. FORTRAN Flow Diagram for Estimator -Extrapolator Subroutine 
(Sheet 2 of 4) 
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Figure 3-8. FORTRAN Flow Diagram for Estimator-Extrapolator Subroutine (Sheet 3 of 4) 
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Z(l) = CALPH (5,K) + CALPH (1,K) + TK*(CALPH (2,K) + TK*(CALPH (3, K) + TK-CALPH (4, K) ) ) 



Figure 3-8. FORTRAN Flow Diagram for Estimator-Extrapolator Subroutine (Sheet 4 of 4) 
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In paragraph 4.1 it will be noted that the present FORTRAN version of 
the estimator-extrapolator subroutine, LT03, has been timed at 17 milliseconds 
on the CDC-3200. The major reason for this long running time is the use of 
inultiply-subscripted arrays; inspection of the assembly language program gen- 
erated by the FORTRAN 3200 compiler discloses the large amount of arithmetic 
computation being done in order to convert, say, a triple subscript on a 
variable to a single address modifier number which then must be loaded into 
an index register before the actual program computation involving that triply 
subscripted variable can be done. A relatively conservative estimate of the 
time which is consumed in this subscript conversion is 6 milliseconds; hence 
actual computation within the LT03 subroutine consumes less than 11 milli- 
seconds; how much less is not easily determined because of the time spent in 
indexing the DO loops on the subscripts plus the time spent in branching cal- 
culations necessitated by sharing of the estimator equations. However, an 
estimate of 9 milliseconds on the CDC-3200 for the estimator computations 
alone is, again, adequately conservative. On the PDP-10 or Sigma 5, the run- 
ning time would be 6 milliseconds, within allowances for real-time operation 
(refer to paragraph 4.1). The key conclusion now is that the estimator would 
run fast enough if no subscripting, and hence no time-sharing of estimator 
equations, were allowed. To achieve this would require the triplicating of 
all groups of computation instructions. Again based on the CDC-3200 compiler 
results, the memory requirements for the estimator subroutine LT03 would in- 
crease from 685 words to 2661 words. Thus, halving the running time requires 
quadrupling the memory requirements. 

Because of the above remarks it is likely that the subscripting pres- 
ently in LT03 will be reduced by simply repeating instruction groups now in- 
dexed by a DO loop parameter. A detailed procedure for accomplishing this 
elimination of subscripts is given in Appendix E. Accordingly, the discussion 
in this section of the detailed flow diagram. Figure 3-8, will be abbreviated. 
Primarily the discussion will refer to the more functional flow diagram of 
Figure 3-7. 

Note that in Figure 3-7 there is only one loop detailed, that one in- 
dexed by K. As K ranges from one to three, the estimator subroutine performs 
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its computations for range, hour angle, and declination angle, respectively. 
The role of K is readily apparent from the detailed branching in Figure 3-8 
(b). In the initial block in Figure 3-7, INDEX is a running index obtained 
from the main program and denotes the number of 1/64 second intervals which 
have elapsed since the main program first entered the track mode. IK and JIK 
are two counters derived from INDEX and range, respectively, from 1 to 90 and 
from 1 to 30. The IK counter is used only during the startup period, i.e., 
when the first ninety samples are being processed. The JIK counter is used 
continuously to control the estimator subunits and to select the proper one 
of the three sets of coefficients for input to the extrapolator section (see 
Figure 3-8(d)). After the initial block, the diagram in Figure 3-7 is largely 
self-explanatory except for the branching associated with the testing of IFLAG. 

IFLAG within the LT03 subroutine corresponds to MODE in the main pro- 
gram; MODE values of -1, 0, +1 indicate, respectively, that the system is in 
the reacquire, acquire, or track mode. There is no estimator function during 
acquisition; during tracking the estimator operation is detailed in Figure 3-7 
by the straight through flow of the diagram. The alternate path indicated in 
Figure 3-7 as that followed during reacquisition is based upon the following 
design philosophy. 

It is assumed that the transition from one mode to another can occur 
only at the beginning of a 1/64 second interval. Hence the estimator will 
have finished processing, in the track mode, the last measurement before 
IFLAG indicates a transition to the reacquisition mode. When the system 
operation switches to target reacquisition, the estimator subroutine flow 
branches to F, which is the entry point for the evaluation of the fitting 
polynomials for range, hour angle, and declination angle. Since the current 
sets of coefficients are based on the previous tracking measurements, the 
net result is that of extrapolating forward by 1/64 second. If the flow now 
were to be an exit from the subroutine then all future entrys into the sub- 
routine, while the system is still in the reacquire mode, would yield further 
extrapolations using the same set of polynomial coefficients. Such a pro- 
cedure, however, is not advisable when extrapolating from a polynomial fit, 
except for extrapolation periods that are short compared with the fitting 
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period. The reason is simply that beyond the region of fit the polynomial 
has no more "zero-crossings" relative to the data trend curve and its value 
will deviate arbitrarily far from the data trend curve. The avoidance of this 
problem is a simple matter: the extrapolated values are fed back into the 
estimator and used to modify the polynomial coefficients. In this way the 
extrapolated values are treated as measurements, the "zero-crossing" region 
of the fitting polynomials moves forward in time, and the extrapolated data 
trend will not deviate arbitrarily far from the actual data trend. In Figure 
3-7 this strategy is indicated by the branching, after F, back to point C 
(to enter the estimator) and then to point G. 

3.6.2 Simulation Results 

The estimator performance has been simulated by processing measure- 
ments obtained from the nominal tracking data and simulated measurement noise. 
The nominal tracking data is discussed in paragraph 3.2, and the measurement 
noise model in Appendix E. The results of the simulation runs are presented 
in Table 3-4 and in Figures 3-9 to 3-14. The simulation process will be de- 
scribed first, followed by a discussion of the estimator performance. 

Tracking data for range, hour angle, and declination were obtained by 
evaluating, at 1/64 second intervals, the sixth-degree polynomials (see 
paragraph 3.2) and adding to these values a measurement noise quantity (see 
Appendix E) . The simulation covered a total tracking interval of approxi- 
mately 69 seconds; the tracking Interval was actually composed of 49 sub- 
intervals of 90 samples each for a total of 4410 measurement points. In 
addition, each 90-measurement interval was broken into three 30 -measurement 
intervals, for each of which the rms estimation was evaluated. That is, since 
each estimator subunit is in the sequential estimation phase for 30/64 second, 
the estimation errors during that 30 -measurement interval were combined to 
form an rms estimation error for the interval. The estimation errors were 
obtained by forming the difference between the position, velocity, and ac- 
celeration estimates and the true position, velocity, and acceleration as 
obtained by evaluating the nominal tracking data polynomials and their first 
and second derivatives. For comparison, the rms value of the measurement 
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TABLE 3-4 


ESTIMATOR PERFORMANCE AS A FUNCTION OF MEASUREMENT 
AND NOISE PARAMETERS FOR FULL TRACKING PERIOD 
(0 to 68 SECONDS) 



Noise 

Measured 

Measured 

Measured 

Measured 

Data 

Correlation 

RMS 

RMS 

RMS 

RMS 

Quantity 

Time Constant 

Noise 

Position 

Velocity 

Acceleration 


(Second) 


Error 

Error 

Error 



1.04 

1.05 

8.87 

39.0 


1/14.76 

2.06 

2.13 

17.3 

75.5 



2.96 

2.94 

24.2 

104.0 

•N 

o 


0.986 

0.838 

7.12 

31.5 

cu 

CO ^ 

1/29.52 

1.95 

1.59 

13.6 

60.0 

<D e O 


3.08 

2 . 58 

21.6 

94.3 

bO U CD 






a CO 






CO 



0.764 

6.48 

28.7 

o b 

1/44.28 


1.51 

12.9 

57.9 




2,17 

18.1 

80.1 



1.0 

0.985 

8.03 

35.0 


1/14.76 


2 . 02 

16.7 

72.8 

a 



2.89 

23.3 

102.0 

<u <u 






iH CO 






C /o<M 


1.0 

0.827 

6.88 

29.8 

<; (D a 
\ co 0) 

1/29,52 

2.0 

1.61 

13.7 

60.2 

u CO 

n /fj /fi 


3.0 

2.41 

20.2 

89.2 

:r! o 0) 



0.736 



\co \cn 


1.0 

6.04 

26.5 


1/44.28 

2.0 

1.49 

12.7 

56.5 



3.0 

2.14 

18.5 

82.9 



1.0 

0.983 

8.11 

35.0 


1/14.76 

2.0 

2.08 

17.6 

76.6 

*\ 

c a 
o <u 


3.0 

2.98 

24.5 

107.0 

•H CO 

cd /oca 


1.0 

0.852 

7.19 

31.7 

a CD o 

•H\cq <U 

1/29.52 

2.0 

1.66 

13.3 


1-1 ' CO 

o , -v;-. 
<u o o 
Olo'nj 
\co \co 


3.0 

2.42 

20.6 

91.4 







1.0 

0.720 

6.15 

27.5 


1/44.28 

2.0 

1.49 

12.8 

56.6 



3.0 

2.17 

18.5 

81.5 


Note: Units of quantities in table are are seconds, centiineters, and seconds. 
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Figure 3-9. Estimal’or Sinriulation for Hour Angle (Sigma = 3, Tau = 14.76) 
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Figure 3-10. Estimator Simulation for Hour Angle (Sigma = I, Tau = 44.28) 
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Figure 3-11. Estimator Simulation for Declination (Sigma = 3, Tau = 14.76) 
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Figure 3-12. EsHmalor Simulation for Declination (Sigma = 1, Tau = 44.28) 
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Figure 3-13 . Esh’mator SimulaHon for Range (Sigma = 3, Tau = 14.76) 
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Figure 3-14, Estimator Simulation for Range (Sigma = I, Tau = 44.28) 
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noise samples for each 30-sample interval was computed. The results of the 
simulation runs, for different measurement noise model assumptions, are in 
Figures 3-9 to 3-14. Note that the continuous curves were obtained by joining 
points spaced at 30/64 second intervals. As described in Appendix E, the 
measurement noise model is a stationary random process with autoeorrelation 
function 

p(t) = cr^ e"^!*"! (unit^) 

and power spectral density 

= 2 °^ ^2 (uni t^ /Hz) 

T + u 

This power spectral density matches the 1/f trend observed in experimental 
determinatons of angle fluctuation spectra (see Appendix F for a treatment 
of the theoretical and experimental results on atmosphere -induced measurement 
noise). The estimation error curves were obtained for different o’ and t 
values. Increasing o values indicate an increasing noise level, but increas- 
ing T values indicate a decreasing correlation time. It is clear that de- 
creasing o and increasing t will decrease the estimation errors. 

The 147 sets of rms estimation errors on 30-sample Intervals were also 
processed to yield rms estimation errors for the entire tracking interval. 
These total rms errors are arrayed in Table 3-4 as a function of o and t. 

From this table it is easy to see the extent to which the estimation errors 
are influenced by the measurement noise characteristics. 

The values of cr = 3 and t = 14.76 correspond to the anticipated worst 
case conditions, when the elevation angle of the laser beam is less than 10 
degrees above the local horizon. The actual noise conditions may be less 
severe. In any event, as the launch vehicle rises, <7 will decrease signifi- 
cantly (see Appendix F, paragraph 3.1) and t is expected to increase. The 
reasons are, simply, that the encountered turbulence is less severe at 
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moderate elevation angles (O’ decreases) and is less correlated (t increases) 
as the target range increases. Hence it is definitely expected that the 
estimation errors will decrease as time from liftoff increases. 

Prediction of the way in which the measurement noise will change with 

elevation angle and target range is a difficult task at present. As indicated 

2 

in Appendix F, Hodara predicts that will decrease exponentially with in- 
creasing altitude; the experimental results of Kurtz and Hayes indicate that 
T will increase roughly proportionally to target range. However, experimental 
results for the tracter site will be required before any reasonable prediction 
of <y and T variations can be made. 
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SECTION 4 

SYSTEM SOFTWARE IMPLEMENTATION 

4.1 REAL-TIME ASPECTS OF SOFTWARE OPERATION 

As noted in paragraph 2.6 the real-time system implementation could be 
accomplished by use of either a PDP-10 or Sigma 5 computer system. The simula- 
tion effort (see Section 3) was aimed primarily at designing the various 
algorithms so as to ensure satisfactory system performance in terms of track- 
ing accuracy requirements, and to guarantee interfacing compatibility among 
the several subroutines. Running time estimates were then made for each 
principal subroutine in order that real-time operation of the system program 
would be assured. The process used in the determination of these estimates 
will now be described. 

Each of the principal subroutines involved in the system simulation 
effort was timed separately on the CDC-3200 by running the subroutine ten 
thousand times and noting the elapsed time by means of the computer's real- 
time clock. The resulting average running times are given in Table 4-1. 

Also included in the table is a projection for the running time on a PDP-10. 

The basis of this projection is the fact that arithmetic operation speeds 
for the PDP-10 are approximately one-third greater than those for the CDC-3200. 
Neglected is the fact that the PDP-10 is, in reality, a much more efficient 
machine than the CDC-3200, in that it has sixteen fast registers and a reper- 
toire of 385 instructions. Thus it is expected that the realizable running 
times are, in fact, less than those indicated. The column labeled "perent- 
age computation time" accounts for the percentage of time that must be devoted 
to each subroutine per clock cycle (1/64 sec). Note that the range servo 
operates at a rate of 256 sec"^ and thus operates four times in every clock 
cycle. 

The overall system flow chart appears in Figure 4-1. This flow chart 
includes the program initialization and calibration steps necessary for proper 
program implementation. The key subroutines tested in the system simulation 
are labeled. Input and output quantities for each of these subroutines are 
given in Figure 4-2. Each of the subroutines is currently FORTRAN coded. 
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In implementing the real-time program it is not expected that any of the 
FORTRAN-coded subroutines need be converted to machine -language coding. In 
particular this is true for the PDP-10 system whose FORTRAN compiler has, on 
Other projects at ARL, generated extremely tight machine coding, so that less 
than 10 percent reduction in running time can be achieved by recoding in as- 
sembly language. 


TABLE 4-1 

RUNNING TIMES FOR PRINCIPAL SUBROUTINES 


Subroutine 

Angle Servo 

Range Servo 

Coordinate Con- 
version 

Estimator/Ex- 

trapolator 


Measured (CDC-3200) 
1.37 milllsec 
1.035 mlllisec 

11.0 millisec 

11.0 mlllisec 


Proiected (PDP-10) 
0.9 milllsec 
0.675 mlllisec 

4.0 millisec 

7.0 milllsec 


Percentage 

5.8fo 

17.3^ 


25.6fo 


48 % 


See paragraph 3.6 for detailed discussion of timing estimate for this 
subroutine. 
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4 . 2 SOFTWARE FOR OFF-LINE DATA REDUCTION 

The real-time data processing is restricted to a cubic polynomial fit 
to a 90/64 second data span. In off-line processing no such restriction is 
necessary. FORTRAN programs have been provided which can accomplish weighted 
least-squares fitting of higher degree polynomials to longer data spans. 

These programs, in subroutine form, are described and listed in Appendix A. 

The major subroutine is POLFITW, which accomplishes the weighted least- 
squares polynomial fit with the aid of GINV (for generalized inverse), 

CHOLSQR (Cholesky decomposition and Inversion of lower triangular matrix) , R 
(yielding the correlation, or weighting, matrix) and POLORT which orthogonal- 
izes cubic polynomials on a particular time interval. The theory behind the 
use of these programs is in Appendix E. 
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WE-WtOG«AM INITIAitZATtON AND CAUBKATION 



Figure 4-1 . System Operation Flow Chart Subroutines Discussed in Text Are Labeled LTOl Through LT05 
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LTOl 


LTD 2 


COORDINATE CONVERSIONS AND REFRACTION CORRECTIONS 


INPUT DATA; 


A A A 

r, r, r 


^HA' ^HA’ ®HA 


FROM ESTIMATOR (LT03) 


^DEC" ®DEC" ^DEC 


SIN 0^^, COS SIN 0^^^, COS 0 




DEC' 


^DEC 


OUTPUT DATA; 


9, 0, 0 

Hv ^EL’ 
^AZ' Kz’ 



LAUNCH- SITE CENTERED 
SPHERICAL TOPOCENTRIC 
COORDINATES 


RANGE SERVO INITIALIZATION AND CORRECTION 
INPUT DATA ; 

SYSTEM MODE 

PROGRAMMED RANGE, OR ESTIMATED RANGE 
OUTPUT DATA ; 

THREE RANGE CHANNEL OUTPUTS (PHASE CONTROL SIGNALS) 


Figure 4-2. Subroutine Input/Output Quantities for Simulation Program (Sheet 1 of 3) 
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LT03 ESTIMATOR/EXTRAPOLATOR 
INPUT DATA : 

SYSTEM MODE 

®HA’ ®DEG 

( EXTRA POLATOR OUTPUT WHEN IN REAGQUISITION MODE) 
OUTPUT DATA; 


POLYNOMIAL 
GOEFFIGIENTS 
FOR GENERATING 
RANGE AND ANGLE 
ESTIMATES 


FROM 

ESTIMATOR 


r, r, 


^HA^ ^HA’ ^HA 


^DEG’ ^DEG’ ^DEG 


FROM 

EXTRA POLATOR 


LT04 ANGLE SERVO GOMPENSATION AND GONTROL 
INPUT DATA ; 

SYSTEM MODE 

A A 

^Ha# ^dec velocity feed-forward) ( TRAGK, AGQUIRE) 

^HA’ ^DEG (reacquire) 

^^HA’ %EG 

^HA" ^DEG C^CK, AGQUIRE, REAGQUIRE) 

^HA’ ^DEG (acquire, REAGQUIRE) 

PREPROGRAMMED 9 ^, (AGQUIRE) 

OUTPUT DATA ; 

HOUR ANGLE AND DEGLINATION ANGLE 
SERVO DRIVE SIGNALS 

Figure 4-2. Subroutine Input/Output Quantities for Simulation Program (Sheet 2 of 3) 
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LT05 RANGE SERVO CONTROL 
INPUT DATA ; 

THREE RANGE CHANNEL INPUTS 
OUTPUT DATA ; 

THREE RANGE CHANNEL OUTPUTS ( PHASE CONTROL SIGNALS) 
r (CALCULATED RANGE) 


Figure 4-2. 


Subroutine Input/Output Quantities for Simulation Program (Sheet 3 of 3) 


4-7 



F-7180-1 


SECTION 5 
CONCLUSIONS 

The control computer program has been completely designed, and the opera- 
tion of all major subsystems tested. Specifications have been set for all 
peripheral equipment. 

All software has been supplied, with the exception of some connective 
commands and some output commands which are best supplied when the actual com- 
puter system is obtained, and the nature of which is fully detailed in the 
flow charts in paragraph 4.1. The design of the angle and range servos has 
been optimized to the point that performance limitations will be due strictly 
to external hardware. In some measure the target dynamics estimator/extra- 
polator is performance-limited by the computer hardware, but it is also con- 
strained by the time constant consideration (maintenance of a low error level 
in the event of an engine failure). More significantly, it is impossible to 
State just what the estimator performance will be, because of the generally 
confused state of knowledge concerning atmospheric noise properties. The 
type of estimator which has been implemented will perform very well in a wide 
variety of noise environments because its design is not tightly wedded to any 
specific measurement noise model. 

When a computer system is obtained it must be expected that further de- 
bugging of the control program will be required. However, short of actually 
building the final system, the simulation and test procedures which have been 
followed in checking subroutine operations constitute the best that can be 
done before actual system implementation. 
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APPENDIX A 

PROGRAM LISTINGS FOR SUBROUTINES 
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program tracker 

laser tracker PR06RAM SIMULATOR 

COMMOK CONZfCON(lOO) tTZtT(50) tXT(3*3) tXL(3t3) (POLAR (3 t 3) 

common IXT(IXTA»IXTb 

COMMON R ( 3 ) ( H A ( 3 } ( DEC ( 3 } t S ( 3 ( 3 ( 2 ) ( OS (3 ( S ) ( RR ( 3 1 2 ) 

COMMON NSAM(KSA«KSB 

COMMOn/daTa/RK(4(50) *HK(4(50) *DK(4#5W) 

COMMON 1HEZ(IHE(3600) (iDEZtlOEtsOO) ilHM(i20) ,idM(30) 
COMMOn/DATa/RQ(7*3) •HQt7*3) (OQ(7t3) 

0ATA(((RQM»J)#l«l»7)»j»l#3)(« 3. 07426539E3( •2.855444664 . 

X. 6487582373c. 0«94641 b 0U( I.930l600lg..3^.I,92e922545E«S( 
X6.73484401lE«8t •. 05053388666 t .02507330732( .Ol437644el5( 
X-1,177956236E^3. I .3949003S3 e- 4( -2.59047835SE-6# l,457Sll935E»8t 

X. 02507330732. .0287S28963( -3.533868708E«3( 5.579601412E-4. 
X-1.295239l775E-5( a.74507l6lE«8( 0,) 

DATA( { (HQ(I.J) (I«l(7) .J«l(3)a 117.4125936. -.08129312544. 

X-. 02099399294. 8. 142128362E-4. •8.858125812E-5. 1.94l3534l6E-6. 
X-l,22l80e5lE-8( 4, 709042989E-2.-. 1330068257. .01715535798. 

X-. 0012969118. 3. 784461 609E»5. •4.6S4506566E-7. 2,060032254E«9« 

X-. 1330068257. 3,43l07l596E*2. -3.8907354E-3( 1.5l37e46436E-4. 

X-2.327253283E.6. 1 .2360l93524E*8. 0.) 

DATA ( ( (00 (1 . J) »1»1 .7) . J«1 .3) * 41.03893121, ,2129026631 . 
X-3.227090683E-2. 3.694090156E-3. •l.l8l95l60SE-4. 1 ,4e0725462E-6. 
X-6.537293629E-9. 7.323323364E-3. 2,449l93662E«2. 1.06749l98E«4. 

XI. 08169538E-4, .7.557305137E-6. 1 .456644425E-7. -8.779346109E-10. 
X. 02449193662. 2. 13498396E-4. 3,2450a6l4E-4. -3. 0229220548E-5 . 
X7.283222125E-7. -5.2676076654E-9. 0.) 

Call setdata(dum) 

Mooe*i 

XKHST*XK0ST»782. 

XkhC1»^K0C1«819.1 
PI«3. 1415926536 
RANGEaRQd.l) 

OEGRAD*PI/l0O. 

RAOSEC«3600,/OEGRAD 
HAFACb2.049236 
OEFAC«. 7162645 
lHEZalDEZ»0 
DO 20 Ial.3600 

20 IHE(I)«0 
00 21 I>1.900 

21 I0E(I)*0 

00 22 Ial.120 

22 iHMdlsO 

00 23 Isl.30 

23 IOMd)sO 
00 5 I«li3 
00 5 ij*1.5 

5 OS ( I fij ) *0 . 

00 6 1*1.3 

00 6 lj*1.2 

6 RR(Iftvi)*0. 

C 

C 

00 1000 ISECS*0.60 
00 1000 lLO0P*^f53 
IN0EX*ISECS«64*IL00P 
DO 1 1*0.3 
Tl*dN0EX*I/4.)/64, 
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T 2 «n*Tl 

T3«Ta*Tl 

T4«T3«TI 

T5*T4#Tl 

T6.T5«T1 

RTa {2» 1) *TI*RQ (3, 1 ) *T2*I^0 (4» 1 )*T3 

« ♦RQ (5 1 1 ) *T4*RQ <6 * I ) #T5*RQ ( Tf 1 ) *T6 

CALL KSERV0(IN0EX,ItRT,Jl,J2,J3,ll,I2,I3) 

CALL LTo5 ( INDEX* I » M » 1 2 * 13 ♦ J1 » J2 » J3 • HANQE * RQNRATE ) 

I CONTINUE 

GO TO (44*46) *SSWTCMF<U 
44 WRITE (61*950) ILOOP* RANGE 
950 F0RMAT(7H IL00P*I4» IX6HRAN6E«E15,6) 

46 CONTINUE 
C 

HTa *H(J(2.l)*n*HQ(3»l)*T2*HQ(4*l)*T3 

* *HQ (5 * I ) *T4*hQ (6 * I ) *T5*HQ (7*1 ) *T6 

D T ■ * OQ ( 2 , 1 ) # T I ♦ DQ ( 3 , 1 ) *T2 * DO ( 4 . 1 ) *T3 

* *0Q(5#l)*T4*0Q(6»l)*T5*0Q(7*l)*T6 
CALL ASERVO(lNOEX«HT*DT*KHfKO.MHE*MOE»MhfMD#MDH,MOD) 

vhsT«whe 

VDST*wDE 

DHIN*WDH 

OOIN«WOD 

I H*«H«RADSEC/ ( XKHST*XKHCI ) 
lnsMO«RAOSEC/(XKOST*XKDCl) 
c#*** temporarily skip the CORRECTIONS.,. 

IF(1.EQ,1)30,31 
31 CONTINUE 

c perform the encoder corrections 

I«lH/360 

F*IH/360..i 

lHaIH*IHE2*IHE(I*l)#(IHE(I*2)-lHE(I«l) )«F/360, 

I.ID/360 

F* ID/360,. T 

in*ID+lOE2*IDE(I*l)*(IDE(I*2)-IDE(l*l) )«F/360, 

c perform the mount corrections 

IalH/10800 

FsIH/lOgOO,.! 

lHaIM*lHM(ia)^(lHM(I*2)-lHM(I*l)r*F/10800, 

IrlD/lOaOO 

F»ID/I0a00,-I 

lDsID*IDM(I*l) *- (lDM(l42)-IDM (l4l) )*F/I0a00, 

c***» Temporarily SKIP the corrections... 

30 CONTINUE 

hour-ih/raosec+hafac 

DECL IN ■I0/RA0SEC*DEFaC 

HOURnhOUR*RAOSEC 

DECLiN«OECLIN*RADSEC 

Range«Range*]'*^°* 

INOEST«INOEX*l 

CALL L703(IN0EST*RANQE*H0UR»0ECLIN*M0DE) 

DO 778 lMs*l,3 
ha<ims)*ha(ims)/radsec 

DEC ( IMS) *DEC (IMS) VRADSEC 
R(IMS)*R(1MS)/100, 

778 CONTINUE 

SH*StN( HA(D) 

SD*«SIN(DEC(1)) 

CHaCOS( HA(l)) 

CO*COS(DEC(l) ) 

HA(1)»HA(1)-HAFAC 


A-3 



DEC(l)*DECa)-DEFAC 


IF(MOCE.LE.O)3,4 

3 VwST*XKHST*XKHCl«HAa) -MH 
VnST«XKDST*XKDCl*OEC < 1 )-M0 

4 CONTINUE 

CALL LT04(|NDEX»VHSTfVOSTiOHINfOOIN>KH»KOtCOI 

Hft(l):«HA(l)*HAFAC 

0EC<1)*DEC(1)*0EFAC 

CALL LTOl (INDEX»SHiCH»SO»CD) 

IF (ILOOP.EOtO) 100*33 
33 IF (IL00P.E0.32) 100,1000 
100 CONTINUE 
WPITE(61,X0) 

in FORMAT( /7X5HINDEX7X7HSIN{HA)8XTHC0S(HA)8X8HSIN(DEC)TX8HC0S{DEC) ) 
WR I TE ( 6 1 , I 1 ) INDEX , SH , CH , SO , CO 
11 format (XI ll, 4X4 (1H/E14. 6)) 

WR I TE ( 6 1 , 9 0 1 ) RANGE , HOUR , OECL IN 
901 FORMAT (X3(1H/El4. 6) ) 

WRITE(61,12) ( (MX»NX*1,3) (NKBita) 

1 ? F0RMAT(/3(3X6HRAN6E(Il,lH)9X3HHA{n,lH)9X4H0EC<Il,lH)5X)) 

WRITE (61, 900) (R(MX) ,HA(MX) ,DEC (MX) ,MX«1 ,3) 

900 FORMAT (X9(1H/E14,6)) 

WRITE(61,907) ( (MX,NX,MXsl,3) ,NX«l,3) 

907 FORMAT (/9(3x3HXT(,M,IH,M,1h|5X) ) 

WRITE (61,900) {(XT (MX,NX) ,MX«1,3) ,NX*1,3) 

WRITE (61,908) ( {MX.NX,MX*1,3) ,NX*1,3) 

908 FORMAT (/9(3X3HXL(,n»lH, II, 1H)5X> ) 

WfllTE(61,900) ( (XL(MX,NX) ,MX«l t3) ,NX»l ,3) 
write (61,909) ( (MX,NX»1,3) *MX«1,3) 

9 0 9 format (/3(5X4hRM0(I1,IH)9X3HEL(II,IH)9X3haZ(11,IH)6X) ) 
WRITE(61,900) ( (POLAR(MX,NX) ,MX«1,3) ,NXiil,3) 

1000 continue 

STOP 

end 
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SUBROUTINE SETOATA(OUM) 

This routine prepares the constant array -con- 

common CONZ»CON(lOO) ♦TZ»T(50) »XT(3t3) »XL(3i3) iPOLAR(3»3) 
COMMON IXT,IXTA,IxTB 

COMMON R(3) »HA(3) *0EC(3» ♦S{3»3ia) lOSOfS) ♦RRf3»2) 

COMMON NSAM»KSA»KS8 

COMMON/OATA/RK(4»50) fHK{4#50) ,OK(4»50) 

C 

DATA( { (RKU fJ) »I*1»4) iJ*lf5) « 



1.1239E 00* 

-8,0322e 00, 

1.5762E 

01, 

-9.1377E 

00* 


1.9054E-Y*!# 

-2.8l30E.Ol, 

-S,3558E-01, 

6.9002E- 

>01* 


7.7n22E-02, 

7.5374E-01, 

-2.721SE 

00, 

1.9947E 

00, 


-7,2046E-03, 

1,4559E 00, 

-4.1077E 

00, 

2.7796E 

00, 

4* 

-6»5360E-02» 

1.R6S4E 00, 

-4,7963E 

00, 

3.U42E 

00) 

OaTA(( (RK(I«J) * lslt4) tJ«6f 

10) « 





-1.0067E-nly 

2.0223E 00, 

-4,8899E 

00, 

3.0676E 

00 , 


« 1 • I635E— ^ ^ t 

1.9667E 00, 

-4,4905E 

00, 

2.7093E 

00, 

ik 

— 1 • l'S63E-01 ( 

1.7387E 00, 

-3.7006E 

00, 

2.1084E 

00, 


-l.t>l73E-0lt 

1,3785E 00, 

-2, 6226 E 

00, 

1 . 3344E 

00, 

# 

-7.7882E-02, 

9.2620E.01, 

•1.3S85E 

00, 

4.S633E-01) 

DaTA( ( (RKCItJ) »I«1|4) ♦J*ll 

»l5)* 





-4, 730 If. 02, 

4,2l88E-0l, 

-1,0930E.02, 

•4.5633E.01, 


.1,32 12E-02, 

-9.4320E-02, 

1.3180E 

00, 

•1,3344E 

00, 


2,U61E«02, 

.5, 8227E— 0 1 , 

2.52S8E 

00, 

-2,1084E 

00, 

« 

5,2S94E-02, 

•1,0019E 00, 

3,5i04E 

00, 

-2.7093E 

00, 


7*7fi65E-02, 

-1.3130E 00, 

4,1693£ 

00, 

-3.0676E 

00) 

DATA ( { «RKU* J) 

«I»lt4) ,Jal6 

, 2 0 ) a 





9,3T40E*O2, 

-1.4755E 00, 

4,4002E 

00, 

-3,U42E 

00, 


9.7022CiR02, 

-1,4493e 00, 

4.1009E 

00, 

-2.7796E 

00, 


8,4462E^02, 

-1.1942E 00, 

3,l690E 

00, 

-l,9947E 

00, 


5*2g46E»02f 

-6.7016E-OI# 

1.5021E 

00, 

«6.9002e< 

•01, 


••2t2584E»0l» 

3,5646E 00, 

-1.1223E 

01, 

9,1377E 

00) 


DATA((<RK(IfJ) 

t 1*1 »4) , J*2i 

,25) a 





« 1 • ^^995E*0 1 9 

1,7565E 00, 

-5.6398E 

00, 

4,7230e 

00, 


^ 1 » Gg28Ew 0 1 9 

1,6648E 00, 

•S.lOSlE 

00, 

4.0555E 

00, 


** ^ 7 G \ f 

1,7480e 00, 

.5,IS69E 

00, 

3,9l76E 

00, 


^ 1 • 2g 34E*» G 1 9 

1.8693E 00, 

-5,3554E 

00, 

3,9277E 

00, 


- 1 *3823EwOI 9 

1.9761E 00, 

-5,5327E 

00, 

3.9482E 

00) 


DATA ( ( (RK <1, J) tl*l f 4) , Ja26 

, 3 0 ) a 





mX 0 ^5 A 0 0 1 9 

2,0458E 00, 

.5.6208E 

00, 

3,9240E 

00, 


•> 1 • A»93se*G 1 9 

2,0722E 00, 

-5,6007E 

00, 

3.8389E 

00, 


• U5n38E»Gl9 

2.0584E 00, 

.5.4el2E 

00, 

3,6972E 

00, 


*»l9^897E-Gl 9 

2,Q127e 00, 

•5,284lE 

00, 

3.5130E 

00, 


••1 *4S83e»G1 f 

1.9446E 00, 

-5,0354£ 

00, 

3.3025E 

00) 

DaTA< ( (RK(IiJ) il*l»4) »J*3l 

*35) a 





-1.4161E-01» 

1.8632E 00* 

-4.7586E 

00* 

3,0806E 

00* 


*1 ,3£86E— Ol t 

1.7758E 00* 

-4.4724E 

00* 

2.8587E 

00* 


-1.3196E-01, 

1.8879E 00, 

•4.i903E 

00, 

2*6446E 

00* 


-1.2721E-01, 

1.6030C 00, 

•3.9210E 

00, 

2*4434E 

00* 

« 

-1.2278E-0I, 

1*5233E 00, 

-3.6695E 

00, 

2.2575E 

00) 

DATA( ( (RKd tJ) 

»I*1*4) *Ja36 

♦ 40) a 





-1.1877E-01* 

1.4500E 00* 

•3r4381C 

00* 

2.0877E 

00* 


.1,1s23e«01, 

1,3835e 00, 

-3.2276E 

00, 

1.9340E 

00, 


.l,i2l5E-0l, 

1,3240£ 00, 

•3.0374E 

00, 

1.7956E 

00, 


.1,09S4e»01, 

1.2709E OO, 

-2.8663E 

00, 

1.6712E 

00, 

« 

.I ,0734e-01 , 

1.2240E 00, 

-2,7l27E 

00, 

1.5597E 

00) 


OaTA( < (RKdf J) 

,45) a 




# 

-1.0S54E-01# 

1.1825E 00, 

•2.5750E 

00* 

1.4597E 

00* 


-l,U409E-0I, 

1.1460E 00, 

•2.4S15E 

00* 

1.3700E 

00* 
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« -1,0295E-01, iai39C 00i •2.3406C OOt 1.2893E OOi 

« .UOSOTE-Ol, 1,08SSC 00, •2,2409E 00, 1.2166E 00, 

« •1,0142£«01, 1,060Se 00, •2,l509g 00, 1.1S09E 00) 

OATA( { (RK Uf J) »I«lf 4) f J**6f 50)« 

* -l,O096E*Olf 1,0383E OOi •2,0694E 00« 1.0913E 00» 

« •.l,0065E-0lt 1.018SE 00* •1.99S2E OOt 1.0372E00* 

« -lt0048C*0l« 1*0007E 00* •1*927SE OOt 9*B774E-01t 

* -1.0040E-0lt 9*8464E-0lt •1.8653E 00* 9t424lE-0it 

* -l.0039E»0lt 9.6993E-01* -ItSOTee 00* 9*0070E-0U 

OATA( { {HKd* j) »i«l*4) *J«l*S)» 

« 1.1239E 00* -8.0322E 00* 1.5762E 01* -9.1377E 00* 

« 1,9054E»01* -2,8l3QE*01* <i5,3SS8E-01 * 6.9002E-01* 

* 7.7022E-02, 7,5374E-0l* «2,72i5E 00, 1.9947E 00, 

« •7,2046E-03, 1,4S59E 00* •4.1077E 00, 2.7796E 00* 

« •6.5360E.02* 1,86S4E 00* •4^7963E 00, 3,1U2E 00) 

DATA( ( (HKdt J) *I«1*4) *Ja6*10)ii 


« .1,0067E«01, 2.0223E 00* 

« .1.1635E-01, 1,9667E 00* 

» .1,1563E-01, 1,7387E 00, 

* -l,0173E-0l* 1.3785E 00, 

« •7,7882E«02* 9,2620E.01* 

OATA ( { (HKdt J) *I«1»4) tJsll 
« .4.7301E-O2, 4,2168E«Ol( 

« .1.3212E.02, -9,4320E»02* 
2,1161E*02* wS , R227E*0 1 * 
« 5,2594E*>02, .1,0019E 00, 

* 7,7865E-02, -I.SUOe 00, 

DATA( ( (HKdtJ) *1*1*4) *J«l6 

* 9,3748E-02. -1,47S5E 00* 

* 9,?022E-02* -1.4493E 00* 

* 0.4462E-O2* -1.1942E 00* 

* 5#2e46E-02* -6*7016E-0l* 

* -2.2584E-01* 3*5646E 00* 

OATA( ( (HKdtJ) *1*1 *4) *J«21 

« -1.0995E*0l* 1.7565E 00* 

« .l.OeaSE-Ol, 1,4648E 00* 
« •1,I723 e* 01, l,74eOE 00, 
« .1,2834 e-01, 1,8693E 00* 

« •l,3a23E«01, 1,9761E 00, 

DATA( ( (HK d .J) *I*1»4) , Ja26 
« •1,4540E.01, 2,04S8E 00, 

« *1,4938E*01, 2,0722E 00* 

« «1,So38E«01* 2,0S84E 00, 

« «1.4897E«01, 2,0127E 00* 

« -1,4583E-01* 1.9446E 00, 

DATA (( (HKdtJ) *1*1,4) ,J*3l 
« «I.4161 e«01, 1,8632e 00, 

« .1,3686E*01* l,77S8E 00, 

« .1.3196E.01, 1.6879E 00, 

« ,1,2721e»01, 1,6030£ 00, 

« .1,2278E»01, 1.5233E 00, 

Data ( ( (HKd*j) *i*l*4) ,j*36 

* -1,1877E«01* 1.4500E 00* 

* -1,1523E»01* l,3835E 00* 

* -1.1215E-01. 1.3240E 00* 

* -lt09S4E-0l» l»2709E 00* 

* -1.0734E-01* 1*2240E 00* 

OATA ( ( (HKd*J) *1*1*4) *J*41 

* •.1.0554E-01* 1.182SE 00* 

« .1,0409e.OI, 1,1460e 00, 

* ,l,0295g,0l, l,1139g 00, 

* .l,0207g.0l, l.OsSSE 00, 


A 


•4,8899E 

00, 

3.0676E 

00, 

•4,490SE 

00, 

2.7093E 

00, 

•3.7006E 

00, 

2.1084E 

00, 

•2.6226E 

00, 

1.3344E 

00, 

•1,3S8SE 
IS) a 

00, 

4.S633E" 

•01) 

•1,0930E« 

• 02, 

•4 .5633E" 

• 01* 

1.3180E 

00, 

•1 .3344E 

00* 

2.5258E 

00, 

•2.1084E 

00, 

3.51 04E 

00, 

-2.7093E 

00, 

4,1693E 

20). 

00, 

-3.0676E 

00) 

4.4002E 

00* 

•3.1142E 

OOt 

4.1009E 

00* 

-2.7796E 

OOt 

3*1690E 

00* 

*l*9947E 

00* 

1*S021E 

00* 

*6t9002E-01. 

•1«1223E 
25) a 

01* 

9* 1377E 

00) 

•5.6398E 

00* 

4.7230E 

00* 

•5.1031E 

00* 

4.0S55E 

00* 

•5.1S69E 

00, 

3,9l76E 

00* 

•5.3554E 

00* 

3.9277E 

00* 

•5.5327E 

30). 

00, 

3.9482E 

00) 

•S.6208E 

00, 

3,9240E 

00* 

.5.6007E 

00, 

3,8389E 

00* 

•S.4812E 

00, 

3,6972E 

00* 

’••5.2841E 

00, 

3,5130E 

00* 

•5,03S4E 
35) a 

00, 

3.302SE 

00) 

•4.7586E 

00, 

3.0806E 

00, 

•4,4724E 

00, 

2.8587E 

00* 

•4,1903E 

00, 

2.6446E 

00* 

•3.9210E 

00, 

2.4434E 

00* 

-3,6695E 
40) a 

00, 

2.2575E 

00) 

•3 .438 IE 

00* 

2.0877E 

00* 

-3.2276E 

00* 

1.9340E 

00* 

•3.0374E 

00* 

1*7956E 

OOt 

•2t8663E 

00* 

1»6712E 

OOt 

•2*7127E 

45)* 

00* 

1 *55972 

00) 

-2.5750E 

00* 

1.4597E 

OOt 

•2.4515E 

00, 

1.3700E 

00, 

-2,3406e 

00, 

1,2893e 

00, 

■2,2409E 

00, 

V.2166E 

00, 


6 
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« -1»014?E*01* 1,0605E 00* •2.1509C 00* 1.1509E 00) 

DATA( ( (HK(I*J) *I*1*4) *J«46*50)* 


* -1.0096E-01* 1.0383E 00* •2.0694E 00* 1.0913E 00* 

* -l,0065E-f>l* 1.018SE 00, •l,9952£ 00, 1,0372 e 00* 

* -1,00418E-01* l,n007E 00, «1,9275E 00, 9*8774E-01* 

« «1.0040e» 01, 9,g464E«01, 01,86532 00, 9,424lE*01, 

* -1,0039E-01, 9.A993E-01, -1,8078E 00, 9,0070Eo01) 

DATA( ( <0K(I*J) *1*1*4) *J*1*5)* 

« 1,1239E 00, -8*0322e 00, l,5762c 01, o9,1377E 00* 

* 1,9054eo01* o2,8130e*01* o5,3558E**01, 6,9002E*01* 

* 7,7n22Eo02* 7.5374E-01* o2,7215E 00, 1.9947E 00* 

« -7,2046E«03* 1.45S9E 00* o4,1077E 00, 2,7796E 00, 

« -6,S360e*02* 1.A654E 00, »4,7963e 00, 3.1142E 00) 

DATA ( < (DK ( I ♦ J) *1*1 *4) *Ja6*10)« 

« -l,0067E-01* 2.0223E 00* -4,8899E OO* 3.0676E 00* 

* -1.1635E-01* 1.9667E 00* -4.49052 00* 2.7093E 00* 

* -1.1563E-01* 1.7387E 00* -3.7006E 00, 2.1084E 00* 

« -1.U173E-01, 1.3785E 00, -2*62262 00, 1*33442 00* 

* -7.7882E-02* 9*26202-01* -1*35852 00* 4*56332-01) 

DATA( ( (OK < I, J) *1*1*4) * Jail *15)* 

* -4.7301E-02* 4*21882-01* -1*09302-02* -4.5633E-01* 

« .1.3212E-02, .9.4320E.02, 1,3180E 00, «1,3344E 00* 

« 2.1161E-02* -5*a 227E-01* 2.52582 00, -2,10e4E 00. 

« 5,H594E-02, -1.0019E 00, 3.5104E 00, •2,7093E 00, 

* 7.7865E-02, -1.3130E 00, 4,16932 00, -3,06762 00) 

DATA ( ( (OK < I, J) *1*1*4) *j8l6*20)* 

* 9.3748E-02* -1,47552 00* 4.4002E 00* -3,1142E 00* 

* 9.7022E-02* -1.4493E 00* 4,10092 00* -2,77962 00* 

» 8.4462E-02, -1*1942E 00* 3*16902 00* -1*99472 00, 

* 5.2P46E-02, -6.7016E-01* 1*50212 00* -6*90022-01, 

* -2.2584E-01, 3*56462 00* -1*12232 01* 9*1377E 00) 

DATA( ( (OKd.J) ,1*1*4) *J*21*25)* 

* -1*09952-01* 1*75652 00* -5*63982 00* 4*72302 00, 

« -1.08282-01* 1.6648E 00, .5,103lE 00, 4,05552 00, 

* -1,17232-01, 1,74802 00, -5,15692 00, 3.91762 00, 

« .1,2834e-01, 1,86932 00, .5,35542 00, 3,92772 00, 

« -1,38232-01, 1,97612 00, -5.53272 00, 3,96822 00) 

DATA( ( (OK(I,J) *I«l*4) ,Ja26,30)s 
« -1.45402-01, 2,04562 00, -5.62082 00, 3.92402 00* 

« -1.49382-01* 2*0722E 00* -5.60072 00* 3*83892 00* 

* -1*50382-01* 2*05842 00* -5*48122 00* 3*69722 00* 

« -1*48972-01* 2*01272 00* -5*28412 00* 3.51302 00* 

« -1*45832-01* 1*94462 00* -5*03542 00* 3*30252 00) 

DATA( ( (DK(I*J) *lal,4) *Ja3l*35)a 

* -1*41612-01* 1.8632E 00, -4,75862 00* 3*08062 00* 

« -1,36862-01* 1,7758E 00, -4,4724E 00, 2.85872 00, 

« -1.31962-01, 1,68792 00, -4,19032 00, 2,64462 00, 

« -1.27212-01, 1,60302 00, -3,92102 00, 2.4434E 00, 

* -1,22782-01, 1,52332 00, -3,66952 00, 2*25752 00) 

DATA( ( (DK(!,J) *lal,4 ) ,Jb36,40)s 

« -1*18772-01, 1.4500E 00, -3, 43812 00, 2,08772 00, 

* -1,15232-01, 1,38352 00, -3,22762 00, 1.93402 00, 

* -1.12152-01, 1,32402 00, -3,03742 00, 1,79562 00, 

« -1,09542-01, 1,27092 00, -2,86632 00, 1,67122 00, 

* -I, 07342-01, 1, 22402 00, -2,71272 00, 1,55972 00) 

DATA( ( (OK(I,J) »Ial,4) ,J*41,45)s 

* -1,0554e-o1, 1, 10252 00, -2,57502 00, 1,45972 00* 

« -1,04092-01, 1*14602 00, -2.45152 00, 1,37002 00* 

* -1*02952-01* IM139E 00* •2O406E 00* 1*28932 00* 

* •1*0207E«01, 1*08552 00* •2*2409E 00, 1*21662 00* 

* -1*0142E*01, 1*06052 00. -2*15092 00* 1*15092 00) 

DATA ( ((DK(I*J) *1*1*4) *J«46*50)* 
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-1.0096E-01. 

1.0383E 00 t 

•2.0694E 

OOf 


-1.0065E-01f 

1.018SE 00. 

•lr99S2E 

OOf 


-1»0048E»01, 

1*0007E OOf 

-1.9275E 

OOf 

# 

-1.0040E-01. 

9»8464E-0l. 

•lt86S3E 

OOf 


•1.0039E-01, 

9.6993E-0lf 

•1.8078E 

OOf 


C0NZ»2t 

C0Na)a. 0001314 

CON( 2)«,000292 

CON< 3) ■2.7182818 

CON( 4) ■•000001 

C0N( ,5) *1,0 

CON( 6) ■*015625 

C0N( 7) ■•000244141 

CON( 8) ■.0078125 

CON (9) ■-22.32509986 

CONUO) ■2283.338565 

CON Ul)«2056. 1279001 

CON (12) ■.99999986 

CON ( 13 ) ■3.59575851 lE-04 

CON ( 14 ) ■a .2212200927E-04 

CON ( IS) ■-3.596324088E-04 

CON (16) *,99999992 

CnN(17)^1.7553777069E-04 

CON ( 18 ) ■-S . 2205886434E-04 

CON ( 19 ) •-1 . 7565359575E-04 

CON (20) ■.99999993 

CON (2U ■.8780847482 

CON (22) *.4785051463 

CON(23)s100, 

CON (24) *20, 

CON (25) ■600. 

C0N(26)^120. 

CON (27) ■640000. 

CON (28) ■3200000. 

CON (29) ■.170 
Con (30) *.034 
CON ( 31) *55600 • 

CON (32 ) *5880. 

CON (33) ■.3 
CON (34) ■a 19.1 
CON (35) ■4700, 

C0N(36)«6.8 
Con (37) ■782, 

CON (38) ■20. 

CON (39) ■•015625 
CON(40)«300. 

CON (41) *18. 

CON (42) *6, 

CON(43)«18. 

CON(44)*10, 

CON(45)*.087 

CON ( 46 ) ■S . 1415926536/ (180. *360 0 . ) 

NSAM 9 O 

KSA*3 

KSBal 

RETURK 

END 

3200 FORTRAN OIASNOSTIC RESULTS • FOR 
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31/32/3300 FORTRAN (2,2)/MSOS 04/25/68 

SUBROUTINE ASERVO ( INDEX .HT,DT»lHt ID»MHE#MDEtMH»MDfMDH, MOD) 

LASER tracker angle SERVO SIMULATOR 

This program calculates the ANGLE»ANGLE ERROR AND ANGLE RATE 
quantities generated from the hardware EQUIPMENT 
iNPUrs ARE- 

1. IH- HOUR ANGLE SERVO DRIVE OUTPUT 

2. ID- DECLINATION ANGLE SERVO DRIVE OUTPUT 

3 ht- target hour angle 
4.- target declination angle 
outputs ape- 

1. MHE- STAP TRACKER HOUR ANGLE SIGNAL 

2. MpE- star tracker DECLINATION ANGLE SIGNAL 

3. MOH- tachometer HOUR ANGLE SIGNAL 

4. MOD- tachometer DECLINATION ANGLE SIGNAL 

5. MH- HOUR ANGLE 

6. DECLINATION ANGLE 

IFdNnEX.EQ.O) 1»2 

1 PTa3,l415926S36 
DEGRACsPI/180. 

TP2al,/1024, 

WW£s400, 

XKHA»XK0A=,3 
BHslOO. 

0Da2O. 

XJH«600, 

XjDal20, 

XKHSTaXK0STa782, 

XKHTaxKDTa20, 

XKHCiaXKDCla8l9,l 
XKHC2aXKDC2a4700, 

XKHC3a,170 
XKDC3a,034 

iHalQsO 
HTMC=CTMCaO, 

HYTMCaDYTMCaO, 

DHlaDDlaHlaDIsO, 

C 

2 CONTINUE 

DO 7 KaO,15 

HXTMCeWME*{XKHA*IH*XKHC3-HTMC) 

HTMCaHTMC*TP2* (3,»HXTMC-HYTMC) ».5 
HYTMCsHXTMC 

0HlaQHl*TP2«(HTMC-DHI*BH)/XJH 

HlaHI^DHI^TPE 

C 

DXTMCaWME* (XKDA»ID*XKDC3«DTMC) 

0TMCsPTMC*TP2* (3,*DXTMC-DYTMC) *.5 
DYTMCaDXTMC 

DDlaOCI*TP2« (OTMC-DDI*BO) /XJD 
0IaDl*0DI«TP2 
7 CONTINUE 
C 

HINaXKHST*XKHCl*HI 

MHaHiN 

MHEaXKHST*XKHCl*DEGRAD*HT-HIN 

MDHaXKHT*XKHC«*DHI 

OINaXKOST«XKOCl*OT 

MpaDIiN 
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MnE»iXKOST*XKDCl*DEQRAD«OT«OXN 

MDD*XkOT*XK0C2*DDI 

RETURk 

END 
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31/32/3300 FORTRAN (2,2)/MS0S OA/25/68 


SUBROtTlNE RSERVO(lNDEX#IfRT.JltJ2»J3,Il*I2f 13) 

Range servo simulator 

THIS PROGRAM CALCULATES THE RANGE CHANNEL PJASE SIGNALS GENERATED 
FRCM THE RANGE SERVO EQUIPMENT 

inputs ARE- 

1. J1 

2. J2 

3. J3 

these are the range channel ERROR SIGNALS 

4. RT- RANGE OF TARGET IN METERS 
OUTPUTS ARE- 

1. II 

2. 12 

3. 13 

THESE ARE THE RANGE CHANNEL PHASE SIGNALS 

IFUNCEX.EQ.O) lt3 

1 IF(I.EQ.0)2,3 

2 TPlsl/20000, 

XKAalO. 

WFalbO. 

ElsE2sE3aO. 

FPl=FB2»FB3a0.0 

Jlsj2sJ3aO 

3 CAaXK««WF*TPl 
CRal,«TPl«WF 
ClaRT«,00l2 
C2aCl*32, 

C3aC2«32, 

DO 7 K*0»79 
FBI=FeUTPl» Jl»0,384 
F02*F82*TPH J2*0,384 
F83sFe3*TPI* J3«0,384 
ElaCA*(Cl-FBl) ♦Cn* El 
E2 sCa*(C2-F82)*CR* E2 
E3sCa*(C 3*F83) ^CB* E3 
7 CONTINUE 
IlsX02.*El 
I2*102,«E2 
I3 b102,*E3 
RETURN 
END 

3200 fortran diagnostic results • FOR RSERVO 


NO errors 
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31/32/3300 FORTRAN (2.2)/MS0S 04/25/68 

SU8R0UTINE LTOl (INOCXtSH* CHiSOf CO) 

VERSION -0- 

PROGRAM LTOlt CALCULATE LASER TRACKER COORDINATE CONVERSIONS 

inputs are- 

1. T-TEMPORARY STORAGE 

2. SO«SINE OF DECLINATION ANGLE 

3. COsCOSINE OF DECLINATION ANGLE 

4. SH.SINE OF HOUR ANGLE 

5. CHbCOSINE OF HOUR ANGLE 

6. C0n*LIST of CONSTANTS 

7. INCEX-COUNT OF BASIC TIME INTERRUPTS- 
(InCRAMENTS every 64TH OF A SECOND t STARTS AT ZERO) 

6. RwRANGE vector (FROM ESTIMATOR- IN METERS) 

9. HA-HOUR angle VECTOR (FROM ESTIMATOR- IN RADIANS) 

10. DEC-OECLINATION ANGLE VECTOR (FROM ESTIMATOR- IN RADIANS) 

11. Place to store xl»xt»polar coordinates. 

OUTPUTS ARE- 

1. XL-LAUNCH CENTERED COORDINATES 

2. xT-TRACKER CENTERED COORDINATES 

3. polar* vehicle polar coordinates 

COMMON CONZfCONaoO) fTZ»T(50) fXT(3»3) .XL (3»3) fPOLAR (3, 3) 

COMMON IXT.IXTA.IXTB 

COMMON R(3).HA(3).OEC(3)#S(3»3,2),DS(3,.5),RR(3,2) 

COMMON NSAM.KSA.KSB 
lE(INCEX.EO.O) 1.2 

1 I*Jal 

GO TO 10 

2 IF(ANc(lNDEX.l).EO.O)3.4 

3 I»1 
J*2 

60 TO 10 

4 1*2 
J*1 

10 00 5 k»1.3 

5 RR(K.I)*R(K) 

TZ*SD«SH 

T(l)«SO*CH 

T(2)*C0*SH 

T(3)*C0*CH 

T(4)*cON(22)«SD 

T(5)*c0N(22)#c0 

T(6)*C0N(21)*S0 

T(7)»C0N(21)»C0 

T(8)«CEC(2)«DEC(2) 

T(9)*HA(2)*HA(2) 

T{10)*T(S)+T(9) 

T(11)«C0NZ*HA(2)*DEC(2) 

T(12)*ha(2)*T(2)*0EC(2)*T(1) 

T(13)*HA(3)*T(2)*nEC(3)«T(l)*T(10)#T(3)*T(ll)#T2 

S(1.1.I)«C0N(21)*T(3)*T(4) 

S(2.l.I)«CON(22)*T(3)«T(6) 

S(3.1.I)*-T(2) 

S(lt2.I)« 0EC(2)*T(5).C0N(21)#T(12) 
S(2.2.I)*-oEC(2)*T(7)-CON(22)#T(12) 

S(3.2,i).dec(2)*TZ-HA(2)*T(3) 

S(1*3.I)*deC(3)*T(5)-T(8)*T(4).T(13)*C0N(21) 

S(2.3.I)*T(8)*T(6)-0EC(3)4T(7)-T(13)#cON(22) 


A- 12 



F-7180-1 


S(3»3.I)*OEC(3)*TZ-HA(3)*TC3)*TUO)*T(2)*Tai»«Ta> 
TZ*SQBTF<S<2iif I)«S<2»l»lMS<3»lil)*S(3*l»in 
ir{TZ,LE.C 0 N{ 4 ) >6»7 

6 T(h*T(2)«T(3)«0, 

60 TO 8 

7 T(4)»RR(1»1)*C0N(1) 

T(5)acON(5)*E)(PF{T(A)*S{l»l»I) ) 

T(5)aT(4)#cON(2)/T(5) 

T(2)«T<5)*S(l*l»I)*S(2,lfl) 

T<1)»-T(5)*TZ*TZ 
T ( 3) aT (S) *S ( 1 * 1 » I ) *S (3, 1 f I ) 

8 DS(l»5)aOS(U4) 

DS(1»4) *05(1*1) 

OR ( 1 *1 ) aT < I ) 

0S(2».5)*0S(2»4) 

DS (2f4) *0S(2f 1) 

0St2*l)*T(2) 

0S(3t5)*0S<3»4) 

OS(3»4) *05(3*1) 

DS(3»1)*T{3) 

0S(1*2)«(DS(1*1)-05(1*5) )*C0N{8) 

OS ( 2 . 2 ) a ( DS ( 2 , 1 ) .05 { 2 , 5 ) ) *CON (8 ) 

OS ( 3 ♦ 2 ) * ( OS < 3 ♦ 1 ) .05 ( 3 , 5 ) ) *CON ( 8 ) 

OS < 1 * 3 ) a ( OS ( 1 » 1 ) *ns ( 1 » 5 ) .CONZ*OS (1*4)) *C0 n (7 ) 

DS (2*3) a (08(2* 1) *08(2*5) .CONZ*OS (2*4) ) «C0N(7) 

DS ( 3 » 3 ) a ( DS ( 3 » 1 ) *05 ( 3 * 5 ) -CONZ*OS (3 * 4 ) ) #C0N ( 7 ) 
XT(l»l)aRR(l*g)«S(l*l*J) 

XT(2*1 )*RR(l*g)*S(2,liJ) 

XT(3*l)aRR{l»g)»S(3*l*J) 

X T { 1 * 2 ) *RR ( 2 ♦ j ) *S ( 1 • 1 * J ) *88 ( 1* J ) *S (1*2 »g ) 
XT(2*2)aRR(2,g)*S(2*l»J)*flR(lij)*S(2*2tg) 

X T ( 3 ♦ 2 ) *RR ( 2 ,g ) «S ( 3 , 1 , J ) *RR ( 1 , J ) *S ( 3 » 2 ,g ) 

XT(l*3)i«RR(3ig)*S(l*l»J)*RR(2»J)4S(l*2»g)#C0NZ*RR(l»J)«S(l»3,j) 
XT(2»3)sRR(3*g)*S(2*l»J)*Rfl(2*J)4S(2*2,g)*eONZ*RR(l*g)*S(2,3,j) 
X T ( 3 » 3 ) aRR ( 3 , J ) ♦$ ( 3 » 1 1 J ) *88 ( 2 ♦ J ) •$ ( 3*2 ,g ) *CONZ *RR ( 1 » J ) *S ( 3 . 3 ♦ j ) 

XU I » I ) *CON (12 ) *XT (1*1) *CON (15) *XT (2*1) *CON { 18 ) *XT (3* I ) *CON (9) 
X(, (2*1) aCON ( 13 ) *XT ( 1*1) *CON ( 16) #XT (2* I ) *CON ( l9) *XT (3*1 ) *CON ( ® > 
XU3f 1)*CON(14)»XT(1»1)*CON{17)*XT(2*1)*CON(20)*XT(3*1)*CON(11) 
XL ( 1 * 2 ) *C0N ( 12 ) a XT (1*2) *CON ( 15 ) *XT ( 2 * 2 ) *CON ( 1 8 ) *XT (3 * 2 ) 
XL{2»2)aC0N(l3)*XT(l*2)*C0N(16)aXT(2*2) *CON (l9) *XT (3* 2) 
XL(3*2)*C0N(14)*XT(1»2)*C0N(17)*XT(2*2) ♦CON(20) #XT (3*2) 

XL ( 1 *3) aCON ( 12) aXT (1 *3) *CON ( 15) #XT (2*3) *CON ( 18) *XT (3*3) 

XL (2*3) aCON ( 13) «XT ( 1 *3) *CON { 16) *XT (2*3) *CON ( l9) #XT (3*3) 
XL(3»3)aCON(14)axT(l*3)*C0N(17)*XT(2*3)*CON{20)#XT(3»3) 

TZ«XL(l,l)aXL(l»l) 

T(l)*XL(2»l)*XL(2*l) 

T{2)aXL(3»l)«XL(3,l) 

T(3)*T(1)*T(2) 

T(4)*SQRTf (T(3) ) 

T{5)aTZ*T(l)*7(2> 

T(6)*XL(1»1)*XL{1.2) 

T(7)aXL(2*l)*XL(2,2) 

T(8)*XL(3*1)*XL(3»2) 

T(9)«T(7)*T(8) 

P0LAH(1,1)*SQRTF(T(5) ) 

POL AR ( 2 , 1 ) *A TANQD ( XL ( 1 * 1 ) * T ( 4 ) ) 

P0LAfi(3,l)aATANQ0(XL(3»l) ♦XL(2,1) ) 

PQLAR(l*2)a(T(6)*T(9) )/POLAR(l»l) 

POLAR (2,2) « (XL ( 1 ,2) 4PGLAR ( ^ » ^ ) -XL ( * » ^ ) *P0LAR (2,1)) 
*/(P0LAR(l*l)aT(4) ) 
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POLAR (3.2) ■ (XL (2# I ) *%L (3i2) -XL (3.1 ) *XL (2.2) ) /T (3) 

POL AR (i 1 3 ) ■ ( XL ( 1 f 2 ) *XL (1 # 2 ) *XL ( 2 • 2 ) ♦XL (2 » 2 ) *XL ( 3 • 2 ) ^XL O » 2 ) 
•♦XL(l.l)*XL(l*3)*XL(2.1)*XL(2.3)*XL(3tl)*XL(3t3) 

•-POLAR (2.1 ) •POLAR (2.1)) /POLAR (1.1) 

XX*T (♦ ) * (POLAR ( 1 . 1 ) «XL ( 1 #3) -POLAR (3. 1 ) •XL ( 1 . 1 ) ) 

YY-POLAR (2,2 ) • (POLAR ( 1 . 1 ) •? (9) •POLAR (2. I ) •! (3) ) 

POLAR (2,3) . (XX-YY) / (POLAR ( 1 * 1 ) *T (3) ) 

POLAR ( 3 , 3 ) ■ ( XL ( 2 . 1 ) *XL ( 3 . 3 ) -XL ( 3 . 1 ) •XL ( 2 . 3 ) -CONZ*POLAR ( 2 , 3 ) •! ( 9 ) ) 

•/T(3) 

RETURN 

END 

3200 Fortran diagnostic results - Eor ltoi 


NO ERRORS 
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31/32/3300 FOflTRAN (2.2)/MS0S 04/25/68 

SUBROUTINE LT03aNDEX»RANQE#HOUHfOECLIN#IFUQ) 

program lto3— target dynamics estimator/extrapolator 

COMMON CONZ,CCN (lOO) iTZ.T(50) iXT(3f3) fXU3t3) tPOLAR (3*3) #IXTf 
X IXTAfIXTB 

COMMON R(3) «HA (3| ,0EC(3) »S (3,3,2) ,03(3,5) ,RR(3,2) ,NSAM,KSA,KS8 
COMMON/DATa/ RK(4,50) ,HK(4,50) ,DK(4,50) 
dimension aLPH(5,3,3) ,AK(A,50,3) ,CAI-PH(S,3) ,M0DE{3) ,Z(3) 
equivalence (RK(l) ,AK<1) > 

IK*INCEX-NSAM*(IN0EX/NSAM) 

IF(IK,EQ, 0)1,2 

1 iKsNSAM 

2 JlKaIN-(NSAM/3)*(lK/(NSAM/3) ) 

IF <JIk.EQ.0>3»4 

3 J!KaNSAM/3 

4 IF (INDEX.LE.NSAM)5,6 

5 lF(IK,LE,NSAM/3)8,7 

6 IF(JIK.EQ.I) 11,12 

T IF(IK.LE.2*NSAM/3)9,10 
H MODEd)*! 

MODE(2)aMODE(3)*0 
GO TO 12 
9 MODE (I) *2 
MOOE(2)*l 
MODE (3) sO 
60 TO 12 

10 M0DE(l)a3 
MODE (2) *2 
MOOE(3)=1 
GO TO 12 

11 MODETaMOOE (3) 

M00E(3)bMCDE(2) 

MOOE(2)aMOOE (1) 

MOOElUaMODET 

12 DO 13 Kal,3 
iFdFLAG.GT.O) 15,14 

15 IF(K.EQ,1) 17*16 

16 IF(K,E(3.2) 18,19 

17 ZNaRANQE 
60 TO 20 

18 ZNaHOUR 
60 TO 20 

19 ZNaDECLIN 

20 Do 21 IEalt3 
MOOOPaMODE(IE) 

IF(MOnOP,EQ,0)2l,23 

23 NPTaJlK*(MOOOP«l)*NSAM/3 
TIEaFL0AT(NPT)/64, 

IF(NPT. £0,1)24,25 

24 ALPH(5,IE*K)aZN 
DO 22 Jal,4 

22 ALPH(iu,IE,K)bO. 

25 IF (NPT,LE,2«NSAM/3) 26,27 

26 NBRaNPT/KSA 
NN8R«NPT-KSA«N8R 
IF (NNBR,E(3,0)28,21 

28 DO 34 Jal,4 

34 ALPH <ijtIE«Kl ■aLPH <J»1E»k1 ♦AK(J,N8R*k1*(ZN-ALPH<5»IE»K11 
GO TO 21 
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27 NBRR»<NPT-2«NSAM/3)/KSB 

NNSRR« (NPT«2«NSAM/3) •KS8«NBRR 
IF(NNBRR*EQ.0)29«30 

29 XPI*ALPH(A,IE*KJ 
DO 31 I*2t4 
iNOaS.I 

31 XPI*ALPH(lND*IEfK)*XPJ*Tie 
NGNaNBR«NBRR 

00 32 Ial*4 

32 ALPH(I*IE*K)aALPH(I»IE*K)*AK{I»N0N»K>*(ZN*XPI-ALPH<5*If iK) ) 

30 00 33 laltS 

33 CALPH(IiK)»ALPH(I.IE«K> 

21 CONTINUE 

36 IF(IFLAQ,6T,0)14»13 
14 TKaFtCAT(2*NSAM/3*JIK)/64, 

Z (1 ) aCALPH (S»K) ♦CALPH( I f K) *TK* (CALPH (2,X) ♦TK* (CALPH (3,K) ♦TK*CAUPH 
X (4,K))) 

Z<2)acALPH<2*K)*TK*(2,*CALPH(3tKUTK*3,*CALPH(4»K)> 

Z(3)ag,*CAUPH(3fK)*TK*6**CALPH(4fK) 

IF(K#EQ.1)40,41 

40 DO 42 l*1|3 
A2 R(L)«Z(L) 

GO TO 47 

41 IF(K.EO. 2)43*45 

43 DO 44 L.al*3 

44 H4(L)aZ(L) 

60 TO 47 

45 00 46 Lal»3 

46 OEC(L)*Z(U 

47 IF(IFLAG.GT.O) 13*48 

48 ZNaZ(I) 

60 TO 20 

13 CONTINUE 
return 

END 

3200 FORTRAN DIAGNOSTIC RESULTS • FOR LT03 


NULL statement numbers 
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31/32/3300 FORTRAN (2,2)/MS0S 04/25/68 

subroutine UT04( Index. VHST fVDSTfOHINf ODIN* IHOUTilOOUTtCO) 

C PROGHV LT04, CALCULATE LASER TRACKER ANGLE SERVO CONTROL* 

I^PUTs ARE- 

1. CONaLIST OF CONSTANTS 

2. INCEXsCOUNT OF BASIC TIME INTERRUPTS- 
<InCRAMENTS every 64TH OF A SECOND. STARTS AT ZERO) 

3. HA»HOUR angle VECTOR (FROM ESTIMATOR- ARCSECSI 

4. dec* OECLINATION angle vector (FROM ESTIMATOR-ARCSECS) 

6. VHST*RAW hour angle star tracker INPUT (GAIN 8ITS) 

7. VOSTaRAW DECLINATION ANGLE STAR TRACKER INPUT (GAIN BITS) 

(VHST OR VOST DIVIDED BY XKHST*XKHCI GIVES RADIANS) 

8. DHINaRAW HOUR ANGLE TACHOMETER INPUT (GAIN BITS/SECONO) 

9. OOINxRaw DECLINATION ANGLE TACHOMETER INPUT (GAIN BiTS/SECOND) 
<UH1N OR ODIN DIVIDED BY XKHT*XKHC2 GIVES RADIANS/SECGNO) 

OUTPUTS ARE- 

1, IHa HOUR ANGLE SERVO DRIVE OUTPUT (BITS) 

2. ID* OECLINATTON ANGLE SERVO DRIVE OUTPUT (BITS) 

COMMON CONZ*C0N(I(»0 ) .TZ.T (50) .XT {3.3) .XL (3.3) .POLAR (3.3) 

COMMON IXT» IXTA.IXTB 

COMMON R(3) .HA(3) .DEC(3) »S(3.3,2) |DS(3.5) ,RR(3.2) 
common NSAM.KSA.KS8 
equivalence (CON ( 27) .XKHl ) 

EOUIVALENCE(CON(2R) .XKDl ) 

E ou I V alence ( con ( 3 1 ) , ehrl ) 

£OUIVaLENCE(C0N(31 ) .EDRL) 

E ou I V A LE N CE ( CO N { 3? ) . EH TL ) 

EQU I valence ( con ( 32 ) .EOTL) 
equivalence (CON ( 34) .XKHCl) 
equivalence (CON ( 34) .XKDCl) 
equivalence {C0N(36) .XKHR) 

EQUIVALENCE {C0N(3A) .XKOR) 

EQU I V ALENCE ( con ( 37 ) . XKHST ) 

EQUIVALENCE ( CON (37 ) * XKOST ) 
equivalence (CON (38) .xkhT) 

E Q U I V A L E N C E ( C 0 N ( 3 A ) . XK 0 T ) 

EQUIVALENCE(C0N(39) .TP) 

EQU I V ALENCE ( CON (41), WHCP ) 

EQU I valence (CON (41 ) .WOCP) 

EQU I VALENCE (CON (42) .WHIP) 
equivalence (C0N(42) ,wOIP) 

EQUIVALENCE(C0N(43) .WHIR) 
equivalence (C0N(43) ,WDIR) 
equivalence (CON (44) .ALPHA) 

EQUIVALENCE (CON (4A) .ASTOR) 


IF(INnEX.EO.O)21.25 
C initialization IF INDEX IS ZERO 
21 HITaHIVaO 

XVHSTaEHlaEHR3aXEHR2aEHR2aO, 

XVOSTaE01aF.OR3aXEDR2aEOR2aO. 

process the hour angle 

25 IF (HIT, NE.”.) 1.3 
1 XaVHST-XVHST 

WHCPLaXKHloTP/ (ABSF(X) *1.) 

IF (WMCPL.6E.WHCP)3.4 
3 WHCPlaWHCP 
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60 TO 5 

4 WNCP1"WHCPL 

5 EHRC* (EHl*WHCPl*VHST)/CD*XKHST*XKHCl*HA(2) 

CALL LIMIT <EHPCfEHRL*HIV) 

C 

IF(HIV«EQ.0.)2«9 
9 Eh2»0, 

2 ehri«ehrc»xkhr*ohin 
mout*ehrUeh2 
C 

CALL LIMIT (HOUTfEHTLiHIT) 

IP(HlT.EQ.0.)6fl0 

6 IF(HIV.EQ.0,)7*8 

7 EH l*EH 1 ♦ , 5#WHCP*WH I P*TP* < 3 • *VHST-XVHST ) 
XVHST-VHST 

8 EHR2*EHRI-EHR3 

EHR3*EHR3t.S*WHIR*TP*(3.*EHR2-XEHR2)/ALPHA 

XEHR3sEHR3 

EH2»( ALPHA-1, )*EHR3 

10 CONTINUE 

PROCESS THE DECLINATION ANGLE 
IF(D1T.NE.0,)U|13 

11 X.VOST-XVDST 
W0CPL*XK01*TP/ ( ABSF (X) *1 . ) 
IF(W0CPL.GE.WDCP)13»U 

13 WDCPIbWDCP 
GO TO 15 
U W0CP1«WDCPL 
C 

15 EnRC«EDl+VlDCPl*VDST ♦XKDST*XK0C1*0EC(2) 

CALL LIMIT (EORCfEORLiDIV) 

IF (OIV.EQ.O.) I2fl9 

19 Ed2*0, 

12 EnRl*EORC«XKDR*ODlN 
D0UT*EDR1*E02 

C 

C A LL L I M I T ( DOU T , EOTL f D I T ) 

IF (OIT.EQ.O.) 16,20 

16 IF(OIV,EQ,0)l7il8 

1 7 EO 1»ED 1* , 5*i«DCP*W0IP*TP* ( 3 , *VDST-XVOST ) 
XVDST«VDST 

10 EDR2*EDR1-E0R3 

EDR3«EDR3*,5*W0IR*TP«(3.*EDR2-XE0R2>/ALPHA 

XE0R2«E0R2 

E02sULPHA-1,)«E0R3 

C 

20 CONTINUE 
IHOUTbHOUT 
lOOUT-OOUT 

return 

END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR LT04 


NO ERRORS 
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31/32/3300 FORTRAN (2,2)/MS0S 

SUBROUTINE LIMIT (VALUE»B0UN0»FLA6) 

LIMIT-SUBROUTINE OF LTOA 
FOR LaSER TBACKER PROGRAM 

IF ( VALUE. 6E. BOUND) 1*2 

1 VALUE«B0UN0 
FLAO-l, 

SO TO 10 

2 IF(VALUE*B0UND.6E,0.)3»4 

3 FLA6*0, 

GO TO 10 

4 VALUE«-BOUNO 

flag*i. 

10 CONTIKUE 
RFTUMN 
END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR LIMIT 


NO ERRORS 


04/25/68 
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31/32/3300 fortran (2,2)/MS0S 04/25/68 

SUBROUTINE LT05 ( INOEXi lMAX3t Uiia* 13* Jl t J2t J3*RAN6E*RGNRATE) 

LT05-RANQE SERVO ROUTINE 

This program computes the range and range rate and the range servo 

ERROR SIGNALS 
inputs ARE- 

1. INDEX«C0UNT of BASIC TIME INTERRUPTS* 

(INCREMENTS EVERY 6ATH OF A SEC0N0#STARTS AT ZERO) 

H, IM4X3«C0UNT OF RANGE INTERRUPTS SINCE LAST INOEX- 

(STARTS AT ZEROtMAX SIZE IS 3tRESETS WHEN INDEX CHANGES) 

3. H» I2» I3«THE THREE INPUTS FROM THE RANGE CHANNELS- 

OUTPUTS ARE- 

I. Jlf J2» J3»THE three OUTPUTS TO THE RANGE CHANNELS. 

2. range* the calculated range in meters. 

3. RGnRATE* The calculated range rate in meters / SEC. 


1F(1NDEX.EQ.0)10»V2 
10 IFaMAXS.EQ.O) lltl2 

C initialization if index and IMAX3 ARE ZERO 
n XEE1*XEE2*XRROOT*XROD*EE2*0. 

XRRa RAN6E«eOO. 

c 

12 IF(IABS(I3) .LE.SIO) lt2 

1 EE la 13 
60 TO 5 

2 IF mars (12) .LE.510)3i4 

3 EElaI2«32 
Go TO 5 

4 EElaM*1024 

5 RR 00 Ta ( EE 1 * 3 0 . ♦ EE 2 ) * 1 2 . 

EE2aXEE2*EEl*3.-XEEl 

xeei»eei 

XEE2aEE2 

RGNRATEaRRDOT*O.000l04l66 
RDDa (PROOT-XRRDOT) *2.6 
X*(ROO*XRDD)*.5 

RRsXRR* (RRD 0T*3.-XRRD0T) *0,000 166666 
RAN6EaRR*,00l25 *x*0, 00002547 

e 

XRRDOTaRROOT 

XROOaROO 

ICHNG*(RR*XRR) 

XRRaXRR«ICHNG 

JlaICHNQ*0. 009765625 

J2aICHNG*0. 031250 

J3aICHNG 

RETURN 

END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR LT05 


NO ERRORS 
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POLFITW 


Program to compute the weighted least-squares fit of a polynomial to data 
in arrays AX and AY. The fit starts at IMIN, includes every ( IS) th point 
until IMAX. 

INPUTS 

1. AX- -Array of values of independent variable. 

2. AY--Array of values of dependent variable. 

3. IMIN--Location in arrays AX and AY at which fit is to begin. 

4. IMAX- -Location in arrays AX and AY at which fit is to end. 

5. IS--Interval between those points in arrays AX and AY which are 

to be included in fitting. 

6. NDEGl, NDEG2, NS — Every (NS) th polynomial from degree NDEGl to 

degree NDEG2 is fitted to the data. 

7. NPRINT = -l--No printout 

= 0 — Print coefficients and bounds on error at fitting points 

= +1^ — Print coefficients, error bounds, and table with AX and 
AY arrays, values for dependent variable computed from 
polynomial, and the error at each fitting point. 

Uses subroutines CHOLSQR, MATMP, and GINV, and function R. 
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31/32/3300 FORTRAN (2.2)/MS0S 04/25/68 

SUaROUTINE POLFlTW(AXfAY«IMINtIHAXf IStN0E6l«N0E62fNSfNPRlNTtC) 
dimension AX(1) »AY(1) tA(200) tU ( 100) fAFLAQ ( 10) tATEMP(lO) «C(10) tX(2 
10) «Y (20) tYE(8*20) «YA(20) fSINV(400) 

COMKOK STEMP(200) 

EXTERNAL R 
IF(IMIN) 2l*20t21 

20 IMIN*IMIN*1 

21 NR*(IPAX-IMIN)/IS ♦! 

00 1 I*1»NR 

IXY ■ IMIN*(I-1)*IS 
X(I)«AX(IXY) 

I Y(I)*4Y<IXY) 

DELX-X {2)-X(D 

NCl*NCE6l*l 

NC2*NCEQ2*1 

DO 30 NC«NCltNC2»NS 
YMIN«n, 

YMAX«Or 
DO 2 I«l,NR 
8X»1. 

DO 2 io«ltNC 
INDaI« ( J.1)«NR 
A(INO)«BX 

? R y M n V f T ) 

CALL CHOLSOR(R.SINV«NR.DELX) 

CALL l»^ATMP{SlNV,NP»NRtA»NC»STEMP) 
call filNV(STEPP»U.AFLA6tATEMPtNR»NC) 

CALL VATMP(SINV»NP.NR»Y»lfYA) 
call MATMP(YAf IfNPfSTEMPtNCfC) 

IF(nPRINT) S0,5l,5l 
51 NCM1*NC-1 
PRINT 3.NCM1 

3 FORMAT (/IX,16HC0EFFICIENTS F0Rtl3.63HTM 0E(3REE POLYNOMIAL FIT IN F 
lORM Y«C(0) ♦ C(1)X ♦ C(2)X**2*ETC.) 

PRINT 4, (C(I) *I*1.NC) 

4 format (6(E20,9) ) 

50 CONTINUE 

DO 5 I*1»NR 
YA(I)«C(NC) 

Do 5 ij«2*NC 
IND«NC*1-J 

5 YA(I)bC(IND)*YA(I)*X(I) 

IP»NC*NCU1 

Do 6 IsltNR 

6 YE(IP.I)»Y(I)-YA(I) 

DO 7 I*l|NR 

IF(YMIN-YE(IP»I) ) 9,9.8 
0 YMIN*YE(IP,I) 

60 TO 7 

9 IF(Y£(IP»I)-YNAX) 7,7,11 

II YMAX*YE(IP,I) 

7 CONTINUE 
IF(NPRINT) 30,45,40 

40 print 10 

10 FORMAT(10X,5H X(I),15X,5H Y(I),15X»12H APPROX Y(I),8X,30H ERROR IN 

1 FIT*Y (D-APPROX Y(D) 

DO 13 I«1,NR 

PRI^T 12,X(I) ,Y(I) ,YA(I) ,YE(IP,I) 

12 FORMAT(4<E20,9»3X) > 

13 CONTINUE 
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45 PRINT 46»YMIN»YMAX 

46 FORMAT (/1^»39HERR0R IN FIT CONFINED TO THE INTERVAL lX.eiT ,9 , IM, 

l»lX»El7,9f IX, IH)/) 

30 CONTINUE 
RFTURn 
END 

3200 fortran diaqnostic resulTS • for polfitw 


NO ERRORS 
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GAIN 


This program computes the gain arrays to be used in SUBROUTINE LT03. 
INPUTS 

1. NSAM--For the current version of LT03, NSAM =50. It is the total 

nvimber of columns of the gain array, or the total number of 
measurements actually processed. 

2. NBA TCH- -Presently set at 20. It is the number of measurements 

processed in the batch processing phase. 

3. SIGMA, TAU--The measurement noise is assvimed modeled as (SIQ1A**2) 

*EXP( -TAU*ABS( T) ) . For batch processing, set T = 3/64. 
For sequential processing, set T = 1/64. 

4. NPRINT- -Equals 0 implies no printout. Equals +1 results in printout 

of the gain array. 


OUTPUT 


1. AK--Duinmy variable for external gain array. 
This program calls CHOLSQR, R, MATMP, and GINV. 
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31/32/3300 FOSTRAN (2.2>/MS0S 04/25/68 

S U a R OU T 1 NE 6 A I N ( AK f N S A M , NB A TCH I QM A » T AU * NPR I N T ) 

DIMENSION PUe) »U<16) tAFLAQiA) »ATEMP U> »H (4) »AK(5»150) *AKW(4) f8 (4) ** 
1»F(4) iO(16) *ALPH(4) »HH(240) tS|NV(400) i STCMP (80) f TEMP (80) 

external R 

T*3./64. 

DO 1 I*1»NBATCH 
BX»l. 

XTaFLCAT(I)*T 
DO 1 ij»l»4 
iNQsI* (J-l)*NBATCH 
HH(INC)sBX 
1 8X»BX«XT 

CALL CHOLSQR(RiSINV»NBATCHiT»TAU) 

CALL MATMP(SlNV,NBATCH.N8ATCH»HMt4,STEMP) 

DO 20 laltNBATCH 
DO 20 Jsl,4 
lN'DaI*(J«l)#N8ATCH 
INOTMJ^ n-l )*4 

20 TEMP(INOT)aSTEMP(TNO) 

call MATMP(TEMPf4,N8ATCHtSTEMP*4fP) 

CALL GINV(P»U»AFLA6f ATEMP,4,4) 

CALL 6INV (STEPP, U,AFLAGiATEMP,NBATCH»4) 

DO 21 lalfNBATCH 
DO 21 Jal,4 
lNDaI*(J-l)*N8ATCH 
iMOTaj* (l»l)*4 

21 TFMP(lNDT)aSTEMP(TNO) 

CALL VATMP(TEPP,4,N8ATCHfSINV,NBATCH»STEMP) 

DO 3 1*1,4 
DO 3 ioal,NBATCH 
iN'DaU (J«l)«4 
3 AK(I,ij)bSTEMP(INO) 

DO 6 lal,16 
6 P(I)a(SIQPA#*2)«P(l) 

Tal,/64, 

NRATCRlaN8ATCH>l 

DO 210 KaN8ATCHl,NSAM 

KTaK*40 


Do 211 1*2,4 *» 

211 H(I)a(FL04T(KT)*T)**(I*l) 

CALL PATMP(P»4,4,H,l,8) •# 

CALL PATMP(H,l,4,B,l,F) #* 

DO 212 1*1,4 •« 

212 Akw{I)*8(I)/(F(1)*SIGMA**2) 

DO 213 J*l,4 ** 

213 AK (J,K)*AKW(J) 

CALL PATMP(AKW,4,1,H,4,U) ** 

00 216 1 * 1,16 
216 U(I)*-U(I) 

DO 214 1*1,4 ** 

IND*I+(I«1)«4 ** 

214 U(IND)*1,*U(IND) 

CALL PATMP(U,4,4,P,4,0) *« 

00 215 1*1,16 *# 

215 Pm *0(11 ** 

210 CONTINUE 

IF(NPfiINT,EQ*0) 8,9 
9 PRINT 4 , ( (aK(I,J) ,1*1,4) ,J*1,NSAM) 

4 FORMAT (lx, 2 hAK, 4 £ 19 . 8 ) 
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8 ReTUftN 
END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR 


NO ERRORS 


GAIN 
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CHOLSQR 

This program computes SINV, the inverse of the lower triangular square 
root of the positive-definite matrix R. NR is the size of R. It is assimied 
that R is generated by the function subroutine R (with DELX s D) . 


31/32/3300 fortran (2,2)/MS0S 04/25/68 

SUBROUTINE CHOUSQR (RiSlNVtNR,DELX,TAU) 

DIMENSION SINV(1> f 0 <20»20> »H <20»20> 

00 1 IJ*1»NR 

JmIsJmI 

SUM»0, 

D02 

^ SUM * SUM+G(JfK)**2 

G(J»J)aSQRT(R(J»J»DELX*TAU)«SUM) 
jpi »ij*i 
DO 1 1*JP1»KR 
SUMaO, 

00 3 Kal..JMl 

3 SUMsSUM*Q(I.K)*G( J»K) 

6 ( I f J ) SR ( I f J t DELX . T AU) -SUM 
1 6(ItJ)«G(I*J)/G(J« J) 

004 lal»NR 

4 Hn.I)sl,/6(I»I) 

DO 5 I*2»NR 
iMlal-l 

DO 5 g*ltlMl 

SUMaO, 

DO 6 KaJtIMl 

6 SUM a SUM *Q ( I »K) *H (K» J) 

5 H(IiJ)a«SUM/G(I»I) 

00 7 Ial,NR 

00 7 Ijal (NR 
IND*I* (J-1)*NR 

7 SINV<IND)aH(I»J) 

return 

END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR CHOLSQR 


NO ERRORS 
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R 


This function subroutine generates the correlation matrix R for samples 
from a noise process with autocorrelation function exp(-'IAulT|) . The element 
R(I,J) of R is returned for a sample interval of D seconds. 


31/32/3300 fortran (2.2)/MS0S 04/25/68 

function R (If J*D|TAU) 

A»EXP(-TAU*D) 

NEXP«IABS(I-J) 

RaA**KEXP 

return 

END 

3200 FORTRAN 01 AGNOSTIC ReSUlTS • F0« R 


NO ERROHS 
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GINV 


This program forms the generalized inverse of a matrix. 

INPUTS 

1. A, NR, NC — A is the matrix to be inverted, with NR rows and NG 

columns. The external array must be single subscripted 
and stored by columns. 

2. U- -Bookkeeping array with NC**2 locations. 

3. AFLAG, ATEMP — Temporary storage arrays, each with NC locations 

The transpose of the generalized inverse of A is stored back in A upon 
return. This program calls DOT. 
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31/32/3300 fortran (2.2)/MS0S 

SUBROUTINE 6INV ( A*U» AFLAGt ATEMPvNRf NC) 

dimension A(1) tU(l) tAFLAS(l) tATEMPdl 
DO 10 I*lfNC 
DO 5 g«l*NC 
IND1*I*(J«1)*NC 
IND2«I*(I-1)*NC 
5 UdNOU-O, 

10 U(IN02)al« 

FAC*D0T(NR»A»1»1) 

FAC»1.0/S0RTF(FAC) 

DO 15 1*1 rNR 
A(n*A(l)«FAC 
DO 20 1*1 »NC 
20 U(I)*ua)«FAC 
AFLAG(1)»1,0 
T0Ll*l455l9l5228.F*20 
T0L«TCLl*T0Ll 
DO 100 J*2,NC 
OOTl=COT{NR»A*J*J) 

JM 1 3 Jw 1 

DO 50 L*1»2 
DO 30 K*1«JM1 
30 ATEMP ( K ) *DOT (NRt A * J f K ) 

DO 45 K*l*JMl 
DO 35 I*l*Nfl 
IndI*1* ( J-1)*NR 
IN02»I* (K-1) *NR 

35 A ( INDl ) aA ( INOl ) -ATEMP (K) *A ( IN02) *AFLAQ (K ) 

DO 40 Isl.NC 
IND1«I+ < J-1)*NC 
IN02*I4 (K-U*KC 

40 U ( INDl )aU( INOl) -ATEMP (K)*U<IND2) 

45 CONTINUE 
50 CONTINUE 

D0T2*C0T(N'R«A»J»J) 

IF(OOT2/dOT 1-TOL) 55,55,70 
55 DO 60 Ial,JMl 
ATEMPn)*0, 

DO 60 Kel,I 

INOlaK4(I-l)*NC 

IN02*K*(J-1)*NC 

60 ATEMP(I)aATEMP(I)*U(IN0l)*U(IND2) 

DO 65 1*1, nr 
IN 0l«l4 (J-1)*NR 
A ( INOl ) *0, 

DO 65 K*l,JMl 
IN02*I* (K-l)*NR 

65 A ( INDl ) *A ( INDl ) -A ( IN02) *ATEMP (K) *AFLAG (K) 

AFLAG(J)*0, 

FAC*DCT(NC,U,g»J) 

FAC«1./SQRTF(FAC) 

GO TO 75 
70 AFLA6(J)*1, 

FAC*1./SQPTF(D0T2) 

75 DO 80 1*1, NR 
INO*I* (J-1)*NR 
80 A(lNO)aA(IND)*FAC 
DO 85 1*1, NC 
IN0*I+(J-1)*NC 
85 U(IN0)*U(IND)*FAC 
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100 continue 

DO 130 J«l»NC 
DO 130 I alt NR 
lNDlal*(sJ»l)*NR 
FAC-0. 

DO 120 K-JtNC 
IN02-I* (K»1>*NR 
IN03-J4(K-1)*NC 
120 FAC«FAC*A<IKD2)*UaND3) 

130 A(IN01>*FAC 
RETURN 
END 

3200 Fortran diagnostic results - ^or 


NO ERRORS 


6Inv 
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DOT 

Called by GINV, this function subprogram forms the dot product with 
double precision acctimulation. 


31/32/3300 fortran {2,2)/MS0S 04/25/68 


FUNCTION D0T<NR»A*JC»KC) 

TYPE CFLOATU) DBLfDBLA»DBl.B 
DIMENSION A(l) 

OBL»0, 

DO 5 I*l»NR 
IND1«I*(JC-1)*NR 
IND2*I* (KC»1)*NR 
D8LA*A(1NDI) 

DBLB-/MIN02) 

5 D81.*D8L*D8LA*DBLB 
00T*DBL 
return 

END 


3200 FORTRAN DIAQNOSTIC RESULTS - FOR DOT 


NO ERRORS 
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MATMP 

This program forms the matrix product G = AB where A is L x M and B is 
M X N. In the calling program the actual arrays corresponding to the dummy 
arrays A, B, and C must be single subscripted and stored by coltimns. Double 
precision accumulation of all inner products is used. 


31/32/3300 fortran (2,2)/MS0S 04/25/68 

SUBROUTINE WATMP (A»Lf MfB*N»C) 

DIMENSION All) *B(l) #C<1) 

type CFL0AT(4) OBL*OBLlfOBL2 

00 10 I*1»L 

DO 20 K«1»N 

IN01»I*(K-1)*L 

DRL*0. 

DO 30 J«ltM 

INd2»I*(J-1)«L 

IM03aij* 

0BL1*A(IND2) 

0BL2*R(IND3) 

30 D8L»0PL*DBLl*0aL2 
20 C(INDl)«D8L 
10 CONTINUE 

return 

END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR MATMP 


NO ERRORS 
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POLORT 


This subroutine orthonormalizes the polynomials 1, X, X , and X on the 
interval [0,T], The orthonormal polynomials Pq, P^^, Pg, and Pg are printed 
out . 


31/32/3300 fortran (2,2)/MS0S 04/25/68 

PROGRflM POLORT 
Ts90t/64« 

ALPHI.T/2, 

T0P2«T*(T*( (AUPHl*42)/2.*T*(-2,*ALPHl/3.*T/4.) ) ) 
BnT2«T*(ALPHl**2*T* (-ALPH1*T/3,) ) 

ALPH2bTOP2/BOT2 

BETA2*T*<-ALPH1/2,*T/3.) 

A*»ALPH1-ALPH2 

8BALPhl*ALPH2*BETA2 

TnP3«T*(T*(B*B/2.*T»{2f*A*B/3.*T*( (A*A42,*B)/4.*T*(2 ,*a/5.*T/6.) ) ) 
1 ) ) 

B0T3*T* (B*B*T* U*fl*T* ( (2#*B*A*A) /3.*T* { A/2.*T/5i ) ) ) ) 
ALPH3*TOP3/bOT3 

BETA3b(T*(T«(-ALPh 1*B/2.*T*( <8«ALPHI*A) /3,*T* ( (A-ALPHl) /4,*T/5. > ) ) 
1))/B0 t 2 
AUA*ALPH3 
BlaB-ALPH3<»A-BETA3 
CaALPH3*B*BETA3*ALPHl 
PRINT ltT«ALPHl*A«8iAltBlfC 

I FORNAT {1x»74H0RTHONORMALIZATION of the polynomials (l,X»X«*2fX**3) 
I ON The interval (O.f fFl0,6»lH)//lX»9HP0(X) • 1 /Ix»13hp 1 (X) *X- 
2{fEl8,8»lH)/lXf 16HP2(X) • X4*2 * (iEl8.8»8H)* X ♦ ( tEl8#8»lM) /lX» 1 
36 hP3(X) m X*«3 * (iEl8,e«nH)« X*42 ♦ (•ElB,8t8H)« X ♦ («El8.8tlH) 
4) 

END 

3200 FORTRAN DIAGNOSTIC RESULTS - FOR POLORT 

NO ERRORS 
LOAD»56 

CONTROL TR 03 
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APPENDIX B 

COMPUTER CONTROL OF ANGLE SERVOS 
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SECTION 1 


INTRODUCTION 

Conventional closed loop control of angular pointing servos 
for telescopes and steerable antennas is accomplished by special 
electronic controllers of analog or mixed analog/digital nature. 

A computer program may replace the electronic controller entirely 
and control the servos in a sampled data process. This may re- 
sult in a significant cost saving when a digital computer is used 
as part of a telescope or antenna system, and computer time is 
available. It may also be possible to improve the control char- 
acteristics by computer programming. 

A particular application is reported here of computer con- 
trol to the pointing of a laser tracking telescope. This work 

has been carried out in part under NASA contract No. NAS8-21089 

2 

and applied to a cassegrain telescope and mount to be delive»^.-d 
to NASA by Goerz Optical Co. Inc. Besides making use of available 
computer hardware to eliminate the need for an electronic con- 
troller the primary objective for this application was to: 

a) improve tracking performance by the use of rate feed- 
forward signals available as a result of other data 
processing 

b) to provide additional control for optimum performance 
of the servos under conditions of servo saturation. 

Control equations are developed to achieve the above objec- 
tives and a final program is written. The results of an extensive 
digital simulation study are also reported to demonstrate the con- 
trol stability under all conditions and to support the theoretical 
estimates of the tracking performance. Simulation results of 
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linear and saturating step responses are given and accurate track- 
ing performance is illustrated in tracking simulations of slowly 
moving targets and the tracking of a Saturn V launch. 


1-2 
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SECTION 2 

CONTROL LOOP DESCRIPTION 

A simplified block diagram of the angle servo control loop 
is shown in Figure 1 for one of the two axis. Angular control of 
the telescope is accomplished by use of electric torque motors 
driving the telescope directly (no gear reducers) around the two 
axis of the equatorial mount. A solid state, high output imped- 
ance amplifier which is controlled by the output of a digital-to- 
analog converter provides a drive with an accurate translation of 
digital input to axis torque. 

The motion of the telescope is measured by three sensors for 
each axis. The fir^t is an angle encoder of high accuracy (22 
bits) mounted on each axis. The second sensor is a tachometer 
also mounted directly at each axis provides an analog (voltage) 
signal proportional to axis rates. The third sensor is an electro- 
optical star tracker system providing an analog signal proportional 
to the difference between a telescope axis reference (or boresight) 
angle and the target angle during target tracking. The output of 
the latter two sensors are converted to provide all-digital signals 
for the telescope motion. 

Accurate positioning of the telescope is aecomplished by sam- 
ple data computer control. The main routines of the computer pro- 
gram pertinent to a single axis angular control are indicated in 
the simplified flow diagram to the left of the dotted line in 
Fig. 1. The processes carried out in this program can be summa- 
rized in the following 5 classes: 

1. Signal summation and compensation which are carried out 
for stable and accurate control in the linear dynamic 
range of the angle servos. 
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Figure 1 . Block Diagram of the Heir Angle Servo and Computer Control Program 
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2. Compensation for well-defined transducer sensitivity 
changes. 

3. Compensation for servo saturation to improve control 
characteristics during large signal transients. 

4. Correction of angular (encoder) measurements by means 
of calibration tables and interpolation techniques. 

5. Angular position and rate estimation and smoothing for 
improved readings in the presence of noise. 

The processes in the first three classes above are normally 
carried out by conventional analog servo controllers for antennas 
and telescopes. In 2) and 3), however, we are able to take signi- 
ficant advantage of the flexibility in computer programming due 
to the nonlinear characteristics of the processes involved. All 
of these processes will be discussed in detail in Sections 4 and 

6. The process of 4) and 5) above may be considered separate and 
independant of the servo control loops and are considered in 
Reference 1. 
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SECTION 3 

LINEAR CONTROL SYNTHESIS 

A simplified block diagram of the linear part of the angle 
servo control loops is shown in Fig. 2. The axis drive consists 
of a power amplifier and a linear, direct drive torque motor. The 
amplifier has a high gain current loop resulting in a flat fre- 
quency response characteristic between input voltage and drive 
torque within the amplifier -mo tore bandwidth This charac- 

teristic of the axis drive is described by the transfer function 


K 


Gd(s) = 


;a 

j 


w, 


ME 


s + W, 


MEi 


( 1 ) 


The sampled data control process to be carried out by the 
computer introduces a phaseshift in the control loop of 


0 


s 


-Tpco 


( 2 ) 


where Tp is the sampling interval. The computational burden on 
the computer limits the value of Tp to no less than (1/64) second. 
With the axis drive bandwidth set conservatively at W^ = 400 1/sec, 
a closed rate loop cross-over frequency of = 54 can be obtained 
with a phase margin of 35°. Thus 


3|C 

This is probably the minimum phase margin for acceptable rate 
loop stability. 
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Pointing errors that develop due to torque disturbances such 
as friction and wind loads are reduced substantially when integral 
compensation is used in the rate loop. When these disturbances 
are expected to be of relatively minor importance, such as in this 
case, the compensation may be of the form 


Gg(s) 


= + ”lR/a 


(4) 


which differ from the compensation with of = oo by a finite low fre 
quency gain and reduced additional phase lag for a given value of 
The latter ratio is usually picked at about 3 for near 
optimum compromise between rate loop stability and compensation 
bandwidth. 

Compensation for the position loop is introduced to obtain a 
finite positional error during target acceleration. It can be 
z'eadily shown that this error is 


■ W, W 


IP CP 


(5) 


where is a constant target acceleration. For near optimum con- 
trol the following choice is usually made: 





( 6 ) 


Substituting (6) into (5) we may express the required minimum 

value for W__ in terms of maximum target acceleration and maximum 
OK 

allowed error as follows 
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^CR ^ 1^27. 

When rate feedforward is used in the absence of position loop com- 
pensation it can be shown in a similar way that the requirement 
on W„„ is now 




where it is assumed that the rate command is an accurate deriva- 
tive of 0^. 

The accuracy of the estimated target rate to be used for rate 
feedforward is not known. Also, slowly varying load torques and 
offset disturbances are present in the loops. For those rea- 
sons it is advisable to include position loop integral conpensa- 
tlon in addition to rate feedforward. Using the expected maxima 
of = •13°/sec2 and s 1/3600 degrees In (7) and 

(8) we get the following range for Wp_ : 


37 < < 112 


(9) 


where the lower value is to be applied when a perfect rate command 
is available. It is clear that the value of W^p = 54 rad/sec se- 
lected for stability reasons puts a high requirement on the rate 
feedforward signal if the 1 sec tracking accuracy is to be met. 

It should be noted, however, that the above calculations are very 
conservative since 
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a) both rate feedfoirward and integral compensation can be 
expected to combine to improve the tracking accuracy 

b) the .13°/sec^ acceleration is the peak expected value 
and does not last long enough for the servo error to 
reach its corresponding final value. 

These conclusions are bom out by the simulations reported in 
Section 7. 

Using integral -plus -proportional compensation for the posi- 
tion loop and the relation of (6), we get with = 54 the follow- 
ing : 


K G (s) = 18 
P P 




s+6 


( 10 ) 


and with n = 10 for the rate loop; 




54 


_J_ s-i-18 

Ta 


( 11 ) 
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SECTION 4 

LINEAR SERVO MODEL 

A detailed model of the hour angle servo is shown in Fig. 3. 

An H-subscript is used to denote hour angle constants and variables. 
With these changed to D-subscripts, Fig. 3 is also the model of the 
declination servo. The only significant difference between the 
two servos, other than differences in the constant values, is the 
"constant" which is unity for the declination angle. For the 
hour angle servo this is a gain which corrects for the variation 
of the star tracker sensitivity with the declination angle and 
takes the form 

which is exact for small | zracking errors. All other constants 
and variables for the two angle servos are listed in Tables 1 and 
2 respectively. 

The servos may be conveniently separated in 3 parts: 

a) the axis power amplifier and motor (or "torquer"), 

b) the angle sensors and A-to-D converting equipment, 

c) the computer control program. 

The control program shown in Fig. 3 is an exact model of the 
linear part of the desired control equations and has been used to 
facilitate the design of the computer program as outlined in the 
Appendix. Some approximations have been made. First, the power 
amplifier and motor combination is given a single order electrical 
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lag only. This is a simplification which is valid and use- 

ful when 



(13) 


where is the lower frequency (pole) in the exact amplifier- 

motor transfer function and is the cross-over frequency (or 

approximate bandwidth) of the servo rate loop. Second mechanical 

resonances of the telescope mount are also neglected in the model. 

This second approximation is considered valid because the lower 

2 

resonance for the type of mount to be used is generally at fre- 
quencies of 500 rad/sec or higher. This value is sufficiently 
high so that the resonance has no significant influence on the 
closed rate loop stability or other performance characteristics 
when the rate loop bandwidth is about 50 rad/ sec. 
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TABLE 1 



ANGLE SERVO CONSTANTS 


Hour Angle 




Axis viscous damping 

100 ft-lb-sec 

^HRL 

Limit on Rate Command 

5.56 * 10^ 

®HTL 

Limit on Torque Command 

5880 


Axis Inertia 

600 ft-lb-sec^ 

■Su 

Motor torque per D/A bit (10 bits) 

0.3 ft-lb/blt 


Constant for computation of W„„„_ 

HCPL 

6.4 * 10^ 

^2 

Star tracker sensitivity correction 

1/cOS 

’%C1 

Star tracker conversion gain 

819.1 bits/volt 

\C2 

Tachometer output conversion gain 

4700 bits/volt 

*SlC3 

Rate loop gair| constant 

0.17 

*=HR 

Rate feedback gain 

6.8 

^SlST 

Star tracker sensitivity 

782 volts /rad 

*Sn 

Tachometer sensitivity 

20 volts /rad/ sec 

T 

P . 

Computer sampling period 

1/64 sec 

T 

^MCL 

Motor torque limit 

300 ft -lb 

W 

HCP 

Position Loop cross-over freq. 

18 rad/ sec 

W 

HIP 

Position Loop compensation freq. 

6 rad/sec 

W 

HIR 

Rate loop compensation freq. 

18 rad/sec 

a 

Rate loop integral constant 

10 

^HL 

Rate limit 

. 087 rad/sec 
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Declination Angle 

Same as for hour angle except 


Bjj = 20 Ib-ft-sec 
J = 120 Ib-ft-sec^ 
Kjj = 3.2 « 10® 

■ 4,2 = 1 - 

Kpc3 =0.034 
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TABLE 2 

ANGLE SERVO VARIABLES; BOTH AXES 



^DTOA 


^R2) 

^R3) 

^RC 


E 

ST 

TC 




MC 


W, 


GPL 


W 


CPI 


A 

-• 

e 


6 

o 


6 

o 


e 


T 


Position loop compensation value 
Rate loop compensation value 
Final control program output (digital) 
Rate loop error value 

Rate loop compensation variables 
Rate command 

Star tracker output in volts 
Torque command 

Telescope rate converted to bits/rad/sec 
Telescope position converted to bits/sec 
Motor torque 
Computed limit on 
Modified W„p 

Estimated target rate in rad/sec 
Telescope angle in rad 
Telescope rate in rad/sec 
Known target angle in rad 
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SECTION 5 

CONVERSION PARAMETERS AND GAIN CONSTANTS 

The conversion of the star tracker and tachometer outputs 
introduces random errors with an average magnitude of 1/2 least 
significant bit. In the case of the star tracker, this noise is 
introduced directly at the point of angular measurement and should, 
therefore, be considerably smaller than the desired accuracy. If 
the linear range of the star tracker to be covered by the converter 
is a-i^c-sec and the conversion accuracy desired is 40^ arc/sec, 

we have 



(14) 


where n is the number of|>its of conversion. Substituting 3600 

and 1/4 arc/sec for Ad in (14) and solving for 

2 we get 


arc/sec for 


2 " ^ 1440 


(15) 


Thus, for a 1/4 arc/ sec accuracy the converter must be at least 
14 bits. With a full range of the star tracker output of + 10 V 
the conversion gain with 14 bits is 

^C1 “ = 819.1 bits /volt. (16) 

The total gain between the tracking error and the computer input 
is therefore 
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Kjjst ^ci ~ ^ 819.1 = 6.4 X 10^ bits/radians (17) 

The effect of the noise introduced by the conversion of the 
tachometer output is difficult to estimate accurately, but since 
this noise is being inserted after the position loop gain the 
error resulting should be 




(18) 


where A0 is the conversion noise and where it is assumed that both 
servo loops are well damped. It has been found^ that the maximum 
angular rate required for tracking purposes is well below 5 /sec. 
Using + 5°/sec as the desired range for the rate measurement and 
a conversion accuracy of 14 bits in (18) gives 

< .12 arc/sec 

■ I 

which is below the desired accuracy of 1/4 axe/ sec.. 

With a tachometer sensitivity given at = 20 volts/rad/sec 
the 14 bits conversion yield 


^C2 "■ 


16383 

10°/sec X ^ X 20 volts/rad/sec 


4700 


The rate feedback signal must be inserted into the control pro- 
gram equations with a gain which is exactly equal to the gain ap- 
plied to the angular error signal, Eq. (17). The rate feedback 
"equalizing" gain is therefore 
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j, ^ ^ST ^C1 

*SlC2 


(19) 


Substituting values in (19) we get 


The least significant bit in the digital-to-analog conversion 
determines the smallest discrete value of torque being applied to 
the axis inertia load. In the sample time Tp this value of torque 
will result in an angular motion of 


A0 = A 



"^MCL 

t 2 

1 


2 

2'’-l 

J 


( 20 ) 


where N is the bits of D/A conversion and torque 

limit. The value of A0 determines the accuracy with which the 
axis inertia can be controlled and should again be less than the 
desired 1/4 are/sec, Thus, if (20) is solved for 2^ we get 


2N+1 


T 

. MCL 
^ JA0 



( 21 ) 


Substituting values of J for the declination angle servo we get 
the upper requirement 

N 

2'-^ > 250 
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which is satisfied if N = 8. A 10 bit conversion can be had with 
little additional cost and is used here. Thus 

N = 10 for D/A conversion 

giving 

- 0.3 ft-lb/bit 


The closed loop gain of the rate loop is given by the pro- 
duct of all gains in the loop and must equal thus 


^C3 




*SlA ^T %C2 


or with values 


( 22 ) 


■^C3 = 
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SECTION 6 

COMPENSATION FOR SERVO SATURATIONS 

An axis drive of this type is generally limited by two types 
of saturations. Of these the more important is the torque satura- 
tion characteristics of the electromagnetic torque motors. The 
inherent saturation of the torque motors is very gradual, but a 
very precise limit must be built into the motor drive circuits to 
prevent motor demagnetization. This gives an abrupt torque (or 
acceleration) limit which is very detrimental to servo stability 
and requires some form of compensation. 

The rate of the servo drive is also limited either by the 
motor and the power amplifier or by the linear range of the angular 
velocity sensor. A set of limits, therefore must be designed into 
the servo control circuits (or conttol program) which prevent the 

servos from exceeding t^ese limits without affecting loop stabil- 

i 

ity. The methods designed to implement compensations for both 
types of saturations in a computei^ control program is discussed 
in separate paragraphs below. 


6.1 TORQUE SATURATION COMPENSATION 

The design goals for the saturation compensation are not only 
to prevent a deterioration of servo stability when an abrupt 
torque limit is introduced, but also to minimize the time required 
to bring the servos out of the saturated condition. These goals are 
accomplished by modifying the control program such that the servos 
operate with full torque applied, first of one polarity, then 
switched to the opposite polarity until the servo error is brought 
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to near zero. If the switching time is chosen properly, the servo 
will then enter the linear operating region at approximately the 
correct velocity. 

The pointing servos have integral compensation in both rate 
and position loops. These compensations are effective only during 
linear servo operation and must be removed from the control equa- 
tions during the saturated mode of operation. The method and con- 
ditions under which this removal is done is discussed at the end 
of this paragraph. The integral compensation equations are re- 
moved at the instant of time, say T = 0, when the torque command 
has reached the + or - limit level. A simplified block diagram of 
the servo in this condition is shown in Figure 4. 

Let us assume that over the time interval of interest the 

• m 

target velocity, 0 ^, is approximately constant (note that 0 ^ can 
be either the feedforward signal or, if this signal is not used, 
the output of the integral compensation) . To obtain a stable 
operation during the torque limited period of a transient, we may 
seek to adjust such that the velocity command to the rate loop 
does not exceed the velocity that can be obtained with the maximum 
torque command. Let the corresponding maximum acceleration be A 
and denote the adjusted value of by W^pj^ then 

9g(t)Wcpj_(t) < A t (23) 

where it is asstimed that 

9j(0) = 9 q(0) (24) 
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Taking the derivative of (23) we get 


For the transient condition it is highly probably that 


0 \j » 0 w 
e CPL e GPL 


(25) 


(26) 


We shall make this assumption and attempt to show its validity in 
simulations (see Section 7). We then have that 




A 


(27) 


Equation (27) is valid for both + and - since A changes sign as 
• 

0^(t) changes sign. For palues of the drive torque below the 
limits the system is linear and W^„ is determined from the linear 

Lit 

analysis (see Section 3). During torque saturation the values of 
W_„ must be modified to W _ such that 

uir Crl 


_ (”cp ”cP-”cPL 

' rCPL "cP > «CPL 


(28) 


For the absence of rate feedforward the velocity command to 
the rate loop is obtained from the position loop integral compen- 
sation. When rate feedforward is used the compensation equations 
correct for the inaccuracies in the rate feedforward signal only. 
In either case, however, when the drive torque reaches the limit, 
this compensation must be removed from the control program to 
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prevent a build-up of rate command during the saturated mode of 
operation. At the end of this period the compensation should be 
reinstated with an initial condition equal to the condition at 
the time it was removed. The integral compensation in the rate 
loop provides a torque command during acceleration. This compen- 
sation must also be removed when the torque limit is reached and 
reinstated with the last computed integral value as the initial 
condition. During the saturation period this value is used to- 
gether with the proportional component to determine the torque 
requirement . 


6.2 RATE LIMIT IMPLEMENTATION 

The servo rate may be limited simply by limiting the rate 
command to the servo rate loop. It is necessary under this con- 
dition also to remove the position loop integral compensation to 
prevent a build-up of th| rate command. The compensation in the 
rate loop could remain in tact during this mode of operation. 
However, if the servo is in the saturated mode of operation when 
the rate limit is reached, the response of the rate loop is rela- 
tively slow and the veloeity would overshoot significantly. This 
overshoot can be avoided by resetting (to a low level) the inte- 
gral part of the torque command momentarily upon reaching the rate 
limit. 
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SECTION 7 

SIMULATION STUDIES 

The angle servo control performance was studied in a digital 
simulation which included all linear characteristics of the model 
of Fig. 3 and the nonlinear parts of the control program (see the 
Appendix) . In addition to the nonlinearities built into the con- 
trol program the simulation also included the sampled data process, 
the conversion processes and the axis friction- stiction charac- 
teristics. The purpose of the simulation study was to show that 
satisfactory performance can be achieved under the assumptions of 
low-noise star tracker output and under all possible transient 
conditions. Specifically, the following performance characteris- 
tics where studied: 

1. Small and large signal step responses. 

2. Tracking of slo^ .y moving targets. 

3. Tracking a Saturn V Vehicle during launch (worst case). 

4. The significance of rate feedforward. 

The response to a small step change in the angular command 
is shown in the plots of Fig. 5. The 0.057 degrees (1 mr) step 
is too small to cause torque saturation and the control program 
operates entirely in the linear range. The linear response is 
distorted, however, by the axis friction- stiction which was set 
at 5 and 8 ft-lb respectively. Adequate stability and satisfac- 
tory speed of response is indicated. The response to a step an- 
gular command of 0.57 degrees (10 mr) is shown in Fig. 6. This 
command is sufficient to cause torque command saturation which 
activates the part of the control program that compensates for 
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Figure 5. Small Signal Step Response 
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this saturation. It is evident that the compensation yields a 
near optimum response for this input since the new commanded posi- 
tion is reached with practically no overshoot under the condition 
of maximum acceleration followed by near maximum deceleration. 

When the step command is increased the angular rate of change 
reaches the rate limit. Figure 7 shows a set of plots for a step 
input command of 1.43° (25 mr) . In this case, maximum accelera- 
tion takes place until the rate limit of 5°/sec is reached. The 
rate is limited to this level without overshoot and remains there 
until maximum deceleration is switched on. In this case there is 
a relatively small overshoot at the end of the transient, but no 
instability is indicated. Both types of saturation also occur 
when a large step change of velocity command is applied. This 
type of response is shown in Fig. 8 which again indicates fast 
recovery and good stability. 

The performance of tpe servo during tracking of a slowly mov- 
ing target is illustrated by the plots of Fig. 9. In this run the 
target velocity is varied according to the equation given below 
the plots. The dynamic errors in this example are negligible so 
that the tracking errors developed results from the combined effect 
of the discrete errors in the conversion processes and the friction 
stiction load. The effect of the latter is quite evident in the 
stop-and-go characteristics at very low velocities in Fig. 9. Peak 
errors are found at the point where the velocity changes sign. 

The value of the errors are here approaching 1 arc/sec. 

The tracking of a Saturn V launch has been simulated by ap- 
proximating the expected launch trajectory with a 6 -order poly- 
nomian^ of the form 

0 = Cjj + C^t + Cgt^ + ... Cgt® (29) 
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For the hour angle the coefficients used are: 


= -8.1293125 * 10 “ 
Cg = -2.0993992 ♦ lO"^ 
Cg = 8.1421283 * 10’^ 

= -8.8581258 * lO"® 

G. = 1.9413534 ♦ 10‘® 

5 

Cg = 1.2218085 * 10"® 


which give the following significant peak values: 

== -3.258 deg/sec 0^^ = -0.1269 deg/sec^ 


The above values of rate and acceleration, which are considered 
typical worst case, would result in a peak dynamic error of ap- 
proximately 0.0013 degrees (4.17 arc/sec) without rate feedforward 
and less than 1/2 arc/se| when perfect rate feedforward is used. 
The plots of Fig. 10 show the target angular trajectory, angular 
rate and tracking errors when no rate feedforward is used. The 
peak error developed is close to the 4 arc/sec expected. The 
simulation is repeated in Fig. 11 but with 85 percent of the exact 
rate feedforward signal applied. In Fig. 12 the exact rate feed- 
forward is applied 100 percent. The improvement of the tracking 
error is better than what would be expected. This unexpected im- 
provement is due in part to the fact that the estimate of the ac- 
celeration error is based upon a servo using rate- feedforward 
without the simultaneous action of integral networks. A perfect 
rate signal is not expected, but the plots indicate that slowly 
varying deviations of up to 15 percent of the true target angular 
rate would still yield a performance improvement of 5 to 1. 
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SECTION 8 

SUMMARY AND CONCLUSION 

The linear control equations for the telescope angle servos 
were derived by neans of conventional control theory. These equa- 
tions involve linear integration which could make the equations 
for numerical solutions lengthy and complex if many digits of ac- 
curacy is required. However, the control equations are part of 
a closed loop and inaccuracies are automatically compensated for. 

It has been found sufficient, therefore, to use numerical inte- 
grating techniques which involve the previous and present sample 
(or calculated) values only. 

The flexibility of programming has been used with advantage 
to implement torque and rate limits and to solve the control equa- 
tions for optimum performance during operations when these limits 
are applied. Simulatyms of the complete angle servo indicate 
unusually quick and stable responses in both linear and saturated 
conditions. 

The accuracy of A/D and D/A conversion required to hold the 
discretizing noise down to a level which does not deteriorate the 
tracking accuracy was found to be well within the state-of-the-art. 
An A/D conversion of 13 bits total and a D/A of 9 bits was found 
adequate for an accuracy of 1/4 arc/sec in this application. The 
overall noise introduced by the conversion and computation pro- 
cesses was found by simulation to be approximately 0.5 arc/sec 
peak-to-peak. 

When a good target angular rate signal is available a very 
significant improvement may be realized by so called rate feed- 
forward. Simulations indicate that dynamic errors may be reduced 
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by better than 10:1 if a true rate signal is available. Some 
indication of how much the rate signal may deviate from the value 
was indicated by introducing a fixed 15 percent error in the sig- 
nal used for feedforward. Under this condition a 5:1 improvement 
of the dynamic errors was realized. 
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APPENDIX 
CONTROL PROGRAM 


A flow diagram of the hour angle servo control program is 
given in Fig. 13. The blocks at Level (1) simulates the target 
and the hardware that provides the signals to the computer. Sim- 
ilarily, the block at Level (7) is a set of equations simulating 
the hour angle drive mechanism. (The latter is discussed in Sec- 
tion 4) . The branching at Level (2) is determined by the mode of 
operation which have been picked as follows: 

mode = 1; Automatic Tracking 

mode = 0; Programmed Control 

mode = -1; Reacquisition 


In the case of automatic tracking the star tracker output signal 
is used directly in the control program. In the other two 

Ho i 

modes is calculated or estimated in programs which are de- 

Ho i 

scribed in Ref. 1. 


The test at Level (3) will branch the program into a set of 
calculations which are designed to compensate for the limit im- 
posed on the torque output level of the drive motors. This com- 
pensation, which will optimize the servo response under the satu- 
rated condition (see Section 6), determines the value of the posi- 


tion loop gain 


The command signal to the servo rate loop 


is computed at Level (4) and tested to determine if this command 


exceeds the maximum rate for which the servo is designed. The 


gain correction at the beginning of Level (4) is included to 
compensate for the change of star tracker sensitivity to hour 
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angle errors as the declination angle changes. This cross coupling 
does not exist for the star tracker sensitivity to errors in the 
declination angle making unity for all hour angles. 

The computation of the torque command E„rpr> which follows at 

rliO 

Level (5) depends to some degree on the result of the tests at 
Level (4). Thus, if a rate limit has been reached the part of the 
command, caused by the rate loop integral compensation, is 

reset to zero. This eliminates most of the torque command in 
transient conditions (see Section 6.2) and stops the drive accel- 
eration quickly. The torque command is tested next at Level (5) 
and both rate and position loop compensation equations are by- 
passed if the torque exceeds the limits. If the rate limit is 
reached the torque command will drop abruptly as explained above 
and the test at Level (5) will exclude the position loop compensa- 
tion equation to prevent build-up of the rate command. The torque 
command is finally multiplied by a gain constant and then converted 
to an analog (voltage) signal for control of the power amplifier. 
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APPENDIX C 
THE RANGE SERVO 

C.l SHORT DESCRIPTION 

A simplified block diagram at the range servo loop is shown in Figure C-1. 
The input signal to the loop, is an IF signal, generated in the hetero- 

dyne receiver, with a phase shift proportional to the target range, R^. A 
locally generated IF signal f with a phase shift proportional to the 

computed range, Rq, is multiplied with ii^ ^ phase detector circuit to 

provide dc signal, with an amplitude proportional to the difference in 

the phase 0^, between snd The detector output is filtered and 

converted to a binary signal for "sampled data" processing by the computer. 

The computer process pertaining to the range loop can be separated into 
three parts for a simplified description. The first part is a computation 
which provides proportional plus integral compensation for the servo loop to 
effectively eliminate velocity errors, minimize acceleration errors, and as- 
sure good stability. The second part calculates the range rate, Rq, and 
range, Rq, output words. The third process provides the feedback signal which 
enables a sufficiently acdurate control of the phase shift of the locally 
generated IF signal. This is accomplished by controlling, with a digital word 
from the computer, the number of pulses which are to be added to- a train of 
clock pulses. The resulting pulse train is then divided down to yield the 
local IF signal. The later is a 3 to 4 bit binary signal which is used in a 
technique of gain switching to provide the desired multiplication for phase 
detection. 

Two other channels of the range servo are indicated in Figure C-1. Each 
channel has a separate phase control, IF generator, phase detector, filter, 
and sample-hold circuits. Switching between channels is accomplished by soft- 
ware programming only, as elaborated on later in this appendix. The follow- 
ing discussions pertain to the fine channel only unless otherwise stated. 
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Figure C-1 . Range Servo Block Diagram 
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G.2 PHASE CONTROL OF THE LOCAL IF SIGNAL 

The local IF signal is generated by counting in a 4 bit register a pulse 
train of frequency Ngf^^,^^ where Ng is the maximum number of counts (up to 16) . 
If the pulse train is generated in synchronism with the heterodyne generated 
IF signal, the two signals will have a constant relative phase shift. Control 
of the phase shift can be accomplished by addition of a controlled number of 
pulses per second to the pulse train. One method of implementing this is shown 

in Figure C-2 where the control pulse train, f„, is generated by counting a 

«■ 

preset number of pulses of the pulse train, f^^, in response to the computer 
feedback word. In the following is a development of the parameters of the 
circuit of Figure C-2. The frequency of the local IF generated by the circuit 
of Figure C-2 is 


"IFfb 


V2 


1 + 


^0 2 + 


«2 


n 


•1 2U 


(C-1) 


where aQ is the sign bit, a^^ the most, and the least significant bits 

of the computer feedback wo’"d. To give an approximately equal range 


around f_„ we set 
IF 


R 


V2 


1 i 
^ 2 


= f 


IF 


(C-2) 


where a„ = 1 and a. to a - = 0. 
0 1 n-1 

is therefore 

IFfb 


The maximum difference between f_„ and 

IF 


2N.N„ 
max 1 2 


(C-3) 


C-3 



F-7180-1 



Figure C-2. Block Diagram of Phase Control and Local IF Generator Circuits 

(One Channel Only) 
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which corresponds to a maxiimim phase rate-of-change of 




2ir 


max 


(C-4) 


max 


Each pulse of fr,/2 corresponds to a phase increment where 


= 


2ti 

V2 


(C-5) 


But since 


2irf, 


^ = 2R 


M 


(C-6) 


where f^ is the laser beam modulation frequency and c the speed of light, we 
also have that 


A^ = 2AR 


2irf 


M 


(C-7) 


where AR is the corresponding range increment. Equating Eqs. (C-5) and (C-7) 
we get 


^1 ^ ^2 “ 2ARf 


(C-8) 


M 


Equations (C-3) , (C-4), and (C-6) yield 


M 


'1 ^^2 


= 4(R) -t 


(C-9) 


Substituting Eq. (C-9) into Eq. (C-2) and solving for (R)„^^ we get 


[R] 


^1F'= 


max 6 f , 


(C-10) 


M 
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The computer range feedback word is provided every t„ = • 5 — seconds dur- 

S ^IF 

ing which time interval the most significant bit of the word should provide 


the required • Since 


max 


[«',J 


oT1"1aA 

2 A9 


fb-'max 


(C-11) 


we find, by use of Eqs. (C-7) and (C-IO), that 

2"=IV2 


(C-12) 


To hold errors in the range servo down to 0.5 cm or less it is expected 
that AR must be about 0. 5 cm also. Substituting this value for AR and 
fj^ = 29.3 X 10 ” Hz into Eq. (C- 8 ) we get 


= 1024 or 2 


10 


Worst case values of (R) for a Saturn V launch have been found to be 

'max 

about 360 m/sec. Using a value of f_„ =i 256 in Eq. (G-10) yields a value for 

(R) of 
max 

(R) = 437 m/sec 

' ' max 


which is sufficient to cover the peak expected velocity. The value of f. 
from Eq. (C-2) is 


R 


f = 2 /3 per see 


Equation (C-12) then finally yields the following closest (integer) value 
for the number, n, of the computer feedback word: 

n = 10 


SinRilations of the range servo have indicated that AR may be significantly 
larger. 
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C.3 GENERATION OF THE COMPUTER FEEDBACK WORD 

The computer feedback word must be proportional to the rate of change of 
the range word. The integration which generates the phase of the local IF is 
taking place in the 4 bit counter of the circuit of Figure G-2. The feedback 
word cannot be taken directly from the range rate word Rq, however, since this 
would require an integration for the calculation of the range word Rq which 
would be outside the control loop. Instead, the feedback word is generated by 
comparing a stored value of the range Rqj^ with each computed value of R^. A 
feedback word is then generated which is proportional to the difference be- 
tween the two range samples when the difference exceeds AR. R^^^ is then up- 
dated and the process repeated. In brief the following steps are followed: 

1. Rq = Rq + f(Rq) tg I 

2. IDR = (Rq - Rqj^)/^ } 

3. + AR*IDR I (C-13) 

In step number 1 the range word is computed from samples of the range 
rate. In step number 2 the fixed point computer feedback number is computed 
which is equal to the number of increments AR. In step ntimber 3 the value of 
Rqj^ is updated if IDR # 0 giving a new value of Rq^ which has been incremented 
exactly the same amount as the phase circuit. Any drift in the calculation of 
Rq^ due to truncation errors in IDR will be corrected in subsequent calcula- 
tions of IDR. The accuracy of the feedback loop is now assured if the number 
of increments A<j> in one sample period is equal to IDR. This was accomplished 
by proper choice of f^^, N^Ng, and n[Eqs. (C-2), (C-8), and (C-12)]. 

C.4 LOOP COMPENSATION AND DYNAMICS 

Proportional plus integral compensation has been chosen for the range 
loop to eliminate velocity errors. The errors due to range acceleration only 
for the servo can be estimated with good accuracy from 
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R 


( 0 . 0) 

xp cp 


(G-14) 


where is the lead corner frequency in rad/sec of the proportional-plus 

integral compensation, and w the cross-over frequency for the servo loop. 

cp 

Good loop stability is obtained when 


(t). < o' w 

xp — 3 cp 


(C-15) 


Using this upper value of and solving Eq. (C-14) for we get 



(C-16) 


The peak value of R 
acceleration of R^ = 1 cm 


2 

expected is about 10 m/sec . 
we should therefore make 


For an error due to 


0 } > 55 rad/ sec 

cp - 

Phase lag introduced by the sampling time interval will cause significant 

deterioration of stability if co is made larger than 1/3 of the sampling fre- 

cp 

quency. This limit puts an upper bound on of about 85 rad/sec. The low 

pass filter after the phase detection introduces further lag and should have 

a corner at least three times higher than (o . Assuming an upper value of 

0 ) =70 rad/sec will result in a filter corner frequency of 210 rad/sec which 

cp 

would allow little filtering action in the 4 msec sampling period. This dif- 
ficulty can be avoided, however, by the use of sample-hold circuits which are 
synchronized to the heterodyne IF output. 
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G.5 MULTICHANNEL SWITCHING 

Switching between three phase detectors operating with sensitivities 

5 

separated by 2 is accomplished by multiplexing, at each sampling time, all 

detector outputs and converting. The coarse channel output is tested first 

-5 

for the presence of a minimum sample value of 2 times the peak linear de- 
tector output. If this or a greater value is detected the coarse channel is 
selected with appropriate gains for generation of range rate and range. If 
not, the medium channel is tested in a similar manner, bringing the selection 
to the fine channel if the fine detector output is within its linear range. 

(It is estimated that each detector will have a useful linear range of about 
H- 1/2 rad of phase difference which, for the fine detector, corresponds to 
about 40 cm range errors . ) 

C.6 SYNCHRONIZATION AND TIMING 

The computer controlled range servo has been designed with particular 
emphasis on asynchronized operation of the servo hardware and the computer 
I/O and peripheral equipment. This is accomplished by using the phase counter 
of Figure C-2 (with anti-coincidence circuits wherever required to avoid 
pulse splitting) , the multiplexer, and the sample and hold computer input 
equipment. The use of the phase counter assures that the required phase 
shift is provided for with sufficient smoothness within each sample period. 

The computer multiplexer provides the channel switch capability and the sam- 
ple and hold provides an input independent of the sample and hold circuit in 
the range servo hardware. 

The timing of the computer input and subsequent output is not critical 
so long as the interval between input samples and output feedback words re- 
main fixed. If the delay between computer input and output approaches one 
full sample time (~4 msec) however, the stability of the loop will deteriorate. 
This delay should therefore be held to a minimum. 
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APPENDIX D 

COORDINATE CONVERSIONS 

In giving a unified treatment of all coordinate conversions used, both 
in the real-time program and in the program to obtain tracking data, all 
coordinate systems will be referred to an equatorial inertial coordinate sys- 
tem. In particular this coordinate system will have its 1-axis directed to 
the North Pole, its 2- and 3-axes lying in the equatorial plane, and its 2- 
axis passing through the point of intersection between the launch site meridian 
and the equator. The system is depicted in Figure D-1 where 

geocentric declination of launch site, 

geocentric declination of tracker site, 

difference in geodetic latitude between launch and 
tracking sites, 


difference in longitude between launch and tracker sites. 


= 


6 i 

t 




and 


AX i 


Note that AX is measured positive in an easterly direction. 

A topocentric coordinate system is one centered at a point on the geoid 
with its 1-axis in the direction of the local vertical, its 2-axis directed 
south along the local meridian, and its 3-axis completing a right-handed 
orthogonal system. The basis vectors for any topocentric system may be ex- 
pressed in the reference inertial system. To do this let T^ and T^ denote the 
matrices of basis vectors for topocentric coordinate systems at the launch 
and tracker sites, respectively. The results are 


If = 


sin 
cos <l> , 


-cos <j)j^ 0 

sin 0 

0 1 _ 


(D-1) 
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Figure D-1. Reference Geocentric, Equatorial, Inertial Coordinate System 
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and 


= 


s in ^ . 


cos cos AX 

cos sin AX 


-cos <l>. 


sin cos AX -sin AX 

sin sin AX cos AX 


(D-2) 


where 

A 

— geodetic latitude of launch site, 

= geodetic latitude of tracker site. 

The geocentric radius vectors of the two topocenters are 


and 




sin 




-f 

■ ^tf 

cos 






_ 0 

- 





sin 

6 

t 


-1 

^t 

^tt 

cos 

6 

t 

cos 

AX 



cos 

6 

t 

sin 

AX 


(D-3) 


(D-4) 


Let X be the coordinates of any point (actually, the target position) as mea- 
sured in the reference geocentric coordinate system, the coordinates of 

the same point as measured in the launch-site topocentric system, and X^.^. the 
coordinates of the same point as measured in the tracker-site topocentric 
system. Then there obtains both 
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Equations (D"5) may be solved to yield either 


(£t - ^ hht) 


(D-6) 


which takes tracker topocentric to launch topocentric, or 




(D-7) 


which does the reverse. By introducing 


A T 

Ar= l[(r^ 


^a) = 


sin (p^ sin + cos cos <p^ cos AX - 
-sin (p cos (p„ + cos sin (p „ cos AX 

t X t X 


cos sin AX 


and 


(D-8) 


rp ^ rpT rp 

-It - -£ -t 


(D-9) 


Eqs. (D-6) and (D-7) are simplified to 


If At 


(D-10) 


and 



(D-11) 


In order to complete the coordinate conversion required of the real-time 
system program, it is necessary to detail two more coordinate conversions: 
from equatorial to topocentric at the tracker, and from rectangular to spherical 
at the launch site. The later conversion follows readily from consideration 
of Figure D-2; in detail it is 
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(D-12) 


Conversion from equatorial to topocentric occurs in two steps. The first, 
from Figure D-3(a), (b), is 



sin 


DEC 

X = r 

cos cos 0,,. 

— e 

DEC HA 


%EC V 


(D-13) 


which gives rectangular equatorial coordinates. Then, from Figure D-3(c), 


^tt 


sin cos ^ 

t t 


-cos </». 


sin <f) 




(D-14) 


The sequence of operations to convert from equatorial tracking data to 
launch-centered spherical coordinates is given, in order of performance, by 
the following equations: (D-13), (D-14), (D-10), and (D-12). 

To complete the conversion required to render MSFC-supplied trajectory 
data into tracking data, it is necessary to consider one more coordinate con- 
version. In Figure D-4 a launch-azimuth aligned coordinate system is shown, 
relative to the launch- site topocentric system, where 



launch azimuth measured from north. 
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The basis vectors, f^, can be expressed in the launch topocentric system, re- 
sulting in the matrix 


10 d 


1! 

0 

s in 

Krr 

AZ 

“COS 

^2 


0 

cos 

AZ 

sin 

M 

1 


(D-15) 


Now any point in a launch-azimuth-aligned system with coordinates is ex- 
pressible in the launch topocentric system as X^^ = — f" trajectory 

data this launch- azimuth- aligned coordinate system is translated to the geocen- 
ter by using the geocentric radius of the launch topocenter. Hence, the 

basis vectors of the geocentric, launch-derived coordinate system in the 
reference inertial system are given by 

G = (D-16) 


since the axes of the geocentric, launch-derived system are parallel to those 
of the launch-azimuth aligned system. Let X denote the vehicle position as 
given in the MSFC- supplied trajectory data . Then the vehicle position in 
the reference inertial system is 

X = G X = T,F« X (D-17) 

- - -g -g 

Use of Eq. (D-5) gives 


r, = X = Xg 


Solving for X^^ yields 


= tJ (liZi X^ - 


(D-18) 


(D-19) 
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Conversion of MSFC-supplled trajectory data to tracking data is accomplished 
by use of Eq. (D-19) followed by the inverse of Eq. (D-14), 


X = 
— e 


sin 
cos 0, 


■cos 0^ 
sin 0, 


O' 

0 

1 


2tt 


(D-20) 


and the conversions: 


= 


+ ^l2 * 


(D-21) 


"dec = 


■1 Pel 


(D-22) 


and 


. -1 

•s in 


X 


e3 


il- 


2 2 

^e2 ^3 


’ ="e2 > “ 


■HA - ^ (X^3) • 90 , = 0 

-1 PesI 


-sgn (X^ 3 ) “ I 90 + Gos 


V 


x^ + x^ 

^e2 ^e3 


, X^2 < 0 


(D-23) 


Equations (D-21), (D-22), and (D-23) follow from Figure D-3( a) and (b) . 

Expressions for the conversion of the time derivatives of all quantities 
may be easily obtained, since all of the transformation matrices are constant. 

In the coordinate conversion expressions there appear certain constants 
associated with launch site and tracking site locations. The launch site 
constants were obtained, or computed, from MSFC-supplied data for the S-504 
trajectory, and are 
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= 6,373,343.9 meters 
= 28.608420 degrees 

and 


6^ = 28.446962 degrees 

The tracker site used to develop the nominal tracking data was chosen to nearly 

coincide with that possible site which was closest to the launch pad. The 

constants for this specific case are in the heading of Figure 3-2. When an 

actual tracking site location is chosen and surveyed, the constants 

6^, and ^ will have to be entered into the coordinate conversion subroutine. 

The calculation of r^^ and 6 should be made via the standard conversion from 

tt t 

geodetic parameters to geocentric, using as a basis the Apollo reference 
geoid. 
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APPENDIX E 

TARGET DYNAMICS ESTIMATOR 

As discussed in paragraphs 2.5 and 3.6 the estimation procedure consists 
of the least-squares fitting of three overlapping cubic polynomials to 
successive data spans of 90 measurements (at a 64/second rate). In this ap- 
pendix only the principals of operation in fitting one cubic are discussed, 
since the implementation of the overlapping feature is a programming problem. 

Let the measurements be Z^, i = 1,..., 90. Denote by the vector of 
measurements 




(E-1) 


2 3 

which are to be batch processed. In fitting the cubic cKq + a^t + o^t + ofgt 
to the measurements Z^, the mathematical problem is that of finding a param- 
eter vector o' , where 



% “l "2 “3 


so that the set of equations 



or 


(E-2) 


(E-3) 


-B ^ -B 


(E-4) 
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has a solution, in some sense. In particular it is desired to find that 
solution of Eq. (E-4) , which minimizes the weighted error function 


'B 



(E-5) 


where R is the correlation matrix for the measurement noise. It has been as- 
sumed that the measurements are corrupted by additive noise which is stationary 
and has correlation function 

p(t) = <P‘ e“^l*^l (E-6) 


Consequently, the elements of R are 


r. , 
11 


1 


r. . = r, , 

ij 


exp 



(E-7) 


and are computed by the function subroutine R (see Appendix A) . The choice 
for ^ is given by the standard Gauss-Markov formula 



T -1 



(E-8) 


Basically the batch-processing procedure is to compute the gain matrix 



T -1 



T -1 
«B ^ 


(E-9) 


before the measurements are taken, and, while the measurements are being 
taken, to compute the matrix -vector product in Eq. (E-8) one column at a time. 



vector ^ has no importance or interpretation. 
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In the subroutine GAIN, ^ is not computed by means of the direct formula, 
Eq. (E-9), because of notorious numerical difficulties with taking the indi- 
cated inverse. Rather, the generalized inverse is used, which yields better 
numerical results. The generalized inverse can be directly applied to the 
solution of Eq. (E-4) only when the weighting matrix in Eq. (E-5) , is the 

identity matrix. This can be accomplished by factoring the positive-definite 

T 

matrix R as R = S where S is lower triangular; the factorization is by 
means of Cholesky decomposition. Then R~^ = (^'^) where is calculated 

directly from by a back substitution procedure. Cholesky decomposition and 
back substitution are performed in the subroutine CHOLSQR (see Appendix A) 
which yields S given R. Now Eq. (E-5) can be written as 



and Eq. (E-8) as 


(E-10) 



giving, as an equivalent expression for Eq. (E-9), 


(E-11) 



(E-12) 


The operation (O^j i.e., forming the generalized inverse, is performed by 
the subroutine GINV. 


As mentioned in paragraph 4.2 the subroutines GAIN, CHOLSQR, GINV, and R 
may be used for an off-line batch processing of data. To facilitate this off- 
line processing, in the case of performing polynomial fits, the subroutine 
POLFITW is made available. This subroutine performs all of the bookkeeping 
associated with a weighted-least-squares fit and has selectable options for 
printout of results and errors. However, POLFITW fits only combinations of 
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2 3 9 

1, t, t , t , etc., up to t . For higher degree polynomials or for longer 

time intervals it is sometimes advisable to work with a set of polynomials 

which are orthonormalized on the fitting interval. The subroutine POLORT is 

2 3 

supplied to orthonormalize 1, t, t , and t on any interval beginning at 
zero. Hence, in the search for better data smoothing results in off-line 
processing, orthonormalized polynomials may be used to minimize error con- 
tributions from purely numerical sources. 

In the sequential estimation phase of each estimator subunit, the param- 
eter vector estimate ^ is modified as each new measurement is made available. 
In order to sequentialize the Gauss-Markov formula, Eq. (E-8) , it is neces- 
sary to assume the noise on the measurements to be white. In effect, R is 
assumed to be the identity matrix. The sequential Gauss-Markov formula is 
derived in many places in the literature; only the results are given here. 
Using the subscript k to index the measurements and parameter estimates within 
the sequential estimation phase, the updating formula is 

where 


Kfe = Sk.l ^ (Sk Sk.l ^ 

(E-14) 


(i - Sk 

(E-15) 


1 '^k '^l] 



and 


r = mean- square value of measurement noise 
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The times t^ start at t^ = 61/64 and go to t^^ = 90/64. To start requires 
only 


A 

% 


A 

% 


and 

Eo = «b) ^ 

from the batch-processing phase. The gain vectors, are also computed 

within the subroutine GAIN via Eqs. (E-14) and (E-15). The gain values are 
outputted by GAIN in a single array (having the dummy variable designation 
AK(I,J) within the subroutine) with the first 20 columns giving and the 
remaining 30 columns giving ^ for k = 1,..., 30. Note that in the batch 
processing phase not every measurement made available at the 64 sec"^ rate 
is actually used. The vector is constructed by using every third measure- 
ment beginning at 3/64 second and ending at 60/64 second. This is not done 
to reduce computation (since the time must be available an 3 rway every third 
measurement) but to reduce the estimation error due to correlated noise by 
spacing samples further apart. In the sequential estimation phase the re- 
quirement for output every 1/64 second does not allow a similar practice. 
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APPENDIX F 

THE EFFECTS OF RANDOM ATMOSPHERIC TURBULENCE 
ON LASER BEAM PROPAGATION 
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1. INTRODUCTION 

The work reported here was undertaken in support of the de- 
sign of a digital computer program for the real-time control of a 
laser tracker. The tracker automatically follows a moving target 
and provides measurements of target position. In addition to con- 
trolling the tracker the computer is to be programmed to smooth 
the raw position data, to estimate target velocity and accelera- 
tion, and to perform certain coordinate transformations. Random 
turbulence of the atmosphere affects laser beam propagation in 
various ways. In designing the control, smoothing and estimation 
programs, a quantitative statistical description of the atmospheric 
turbulence noise in the tracker target measurements would be most 
valuable, The study reported here had as its purpose ascertain- 
ing whether such a description was available or could be derived 
from available data. 

Theoretical predictions of turbulence effects are discussed 
in Section 2 of this report. The theory treats a number of effects 
as separate and distinct. In practice, however, it is not always 
possible to measure one such effect in isolation. Some of the 
published experimental measurements of laser beam disturbances 
caused by the atmosphere are summarized in Section 3 and compared 
both with the theory and with one another. The conclusions drawn 
from this work are presented in Section 4. 

2. THEORETICAL DISCUSSION 
2 . 1 BACKGROUND 

The effects of random atmospheric turbulence on laser beam 
propagation have been treated from a theoretical viewpoint by 


1 
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12 3 

Beckmann and Hodara . Both assume the validity of Tatarski's 

blob theory, and they use methods of analysis based on geometric 

optics. 

The atmosphere is characterized as having a refractive index, 
n, given by 

n = 1 + Nq X 10"® + 4n 

where 

Nq = the average refract ivity 

and 

M = the random variation in n. 

The average refractivity, Nq, is a function of atmospheric pres- 
sure, temperature and composition. Its consequent spatial varia- 
tion gives rise to a steady (dc) refraction over a given path. 

The random variable, ^(r,t) is a function of both space and time. 
Its effects must be evaluated statistically and it is these ef- 
fects which are of particular interest here. The spatial and tem- 
poral correlation functions of are defined as 

'‘l2('^l''^2) ^ ^(ti) ^(t2) 
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where 


<•> denotes ensemble average 


and 

• denotes time average. 

It is usual to assume that the statistical distribution of M is 
the same at all points in the medium (homogeneity) and that C is 
therefore a function only of p = r„ - r, • I.e. 


Ci2(p) = Mr+p)> 


This correlation function may be characterized by two parameters, 

C, i(0) and L . The first is seen to be the variance of ^ 

11' c 

The second called the correlation distance is the value of p at 
which 

h2(l-c) = 


where a is a predetermined constant. If the meditim is isotropic, 
i.e., if is the same for all directions in space, then it may 
be assumed that C^g(p) functional form 


Ci2(p) 




2 
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Assuming that An( t) is a stationary process, it follows that 
^ 12 (^ 1^*^2 ) " ^ 12 ^'^^ 


= An( t) An( t+r) 

i.e. the temporal correlation function of An depends only on t. 
Like the spatial correlation function R may be characterized by 
R^^(O) and a correlation time, T . If An is ergodie then 

= G^^(O) = <An^> = An^ 


There is reason to believe that the conception of the atmosphere 
as consisting of frozen inhomogeneities (blobs) , whose linear di- 
mensions are proportional to L^, is a meaningful one. This con- 
cept permits establishment of a simple relation between L and T . 

c c 

If the component of relative velocity between the atmosphere and 
the laser beam transverse to the beam is v.^, then 



This transverse velocity component may be due to wind, beam motion 
or both. 

Consider now the effects on the beam, as observed at the re- 
ceiver of a laser communications or tracking system, induced by a 
random turbulent atmosphere as described above. These effects 
will include the following: 
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a) The angle of arrival of the beam at the receiver 
will fluctuate (quivering) . 

b) The location of the beam center will vary (beam . 
wander, spot dancing) . 

c) The cross-section of the beam will vary in size and 
shape (breathing) . 

d) The transit time will vary causing phase fluctuations. 

All these effects result from refraction of the beam as a whole. 

The remaining effects are due to differential refraction of the 
rays making up the beam which causes wavefront distortion. 

e) The image formed by the receiving optics will be 
blurred. 

f) The intensity distribution across the beam will 
fluctuate due to reinforcement and cancellation 
( boiling) . 

g) The phase distribution across the beam will fluctuate 
due to transit time variations. 

From the point of view of a particular laser system several 
of these effects may act in combination and some may be of little 
consequence so far as system performance is concerned. The par- 
ticular system of interest here is a laser tracking and ranging 
system. Its function is to measure the angular position and range 
of a target. The system consists of a cw-laser transmitter and a 
receiver telescope mounted coaxially in servo controlled gimbals. 

A retroreflector on the target returns the transmitted beam to 
the receiver aperture. The receiver optics focus the reflected 
light on a star tracker tube whose two output signals are pro- 
portional to the target image displacement from the receiver bore- 
sight. These signals are used by the gimbal servos to track the 
apparent direction of arrival of the received light. Low frequency 
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intensity fluctuations (fading) induced by the atmosphere are 
minimized by means of automatic gain control in the receiver. The 
transmitted beam is modulated at one or more frequencies and this 
modulation is detected in the receiver. The phase shift of the 
received modulation relative to the transmitted modulation is mea- 
sured in the receiver and provides a measure of the target range. 

For this particular system the atmospheric turbulence effects 
of primary interest are fluctuations in angle of arrival and re- 
ceived phase, since these affect directly the system measurement 
data. Of secondary interest in designing the AGO system are ran- 
dom fluctuations in received intensity. The ability to predict 
these effects would materially aid the design of i) the tracking 
system servo loops, ii) the AGG system, iii) the ranging system 
phase comparator and iv) the algorithms for smoothing the mea- 
surement data. In the following section the theoretical analyses 

1 2 

of these particular effects as derived by Beckmann and Hodara 
are presented and compared. 


2.2 THEORETICAL RESULTS 

Both Beckmann and Hodara use the methods of geometric optics 
in their analyses. The condition under which these methods are 
applicable is 

L » VlX 
c 

where L^ is the correlation distance of atmospheric inhomogene- 
ities 
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L is the pathlength 
X is the wavelength. 

This condition implies that diffraction effects of turbulence blobs 
of width are negligible, i.e., that all points on the optical 

path are in the near field of the farthest turbulence blob. For 
optical wavelengths and practical atmospheric pathlengths this 
condition is usually met. 

2.2.1 Quivering 

The angle of arrival at a receiver may be resolved into two 

orthogonal components A0 and A0 transverse to the direction of 

X y 

beam propagation (z-direction) . The mean square values of these 
components are according to Beckmann 

<a4> = IJ. <A.2> l 

<A8p> = (i L 

For an isotropic medium = L = L = L and 
^ cx cy cz c 

<A0^> = <A0^> = 11 <An^> L/L 

X y c 


where 

jU = 2\Itt 
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The analogous results derived by Hodara are identical to the above 
except that the constant, fi, is 1 instead of 2^. The variance of 
the total angular fluctuation may be taken to be the sum of vari- 
ances of the components. 

<A0^> = <Ael> + <^l> 

X y 

These relations give the variance of the angle of arrival 
of a laser beam at the receiving end of the path. The apparent 
direction to the source will fluctuate with the angle of arrival. 
The effect of this angular variation on the target direction de- 
tector and the pointing servos must be considered in the system 
design. 

2.2.2 Beam Wander 

Just as random variation in the index of refraction causes 

quivering of the apparent angle of arrival of the received beam, 

it also causes the entire beam to wander about an average position 

in the receiver plane. The theoretical variance of this transla- 

2 

tional wander as derived by Beckmann is L /2 times the variance of 
the angular fluctuation, e.g. 

<Ax^> = <A$l> L^/2 
= }i <Ac?>'L^ 

cz cx 

Hodara assumes without proof time the rms wander is L times the 
rms quivering. I.e., 
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Ax = L A0 
rms xrms 

= L 

2 

This implies that the variance of Ax will be L times the variance 
of A0 , which agrees with Bechmann's result to within a factor of 

V5F. 

The effect of beam wander on a tracking system will depend 
on the relative magnitudes of the beam diameter at the receiver, 

r 

the receiver aperture diameter and the rms beam wander. If the 
beam diameter plus the beam wander is small compared to the aper- 
ture, then all of the beam will be collected all of the time. In 
this case the usual type of angle of arrival sensor (star tracker 
tube or image dissector) will sense the beam wander and its effect 
will be indistinguishable from fluctuations in angle of arrival. 

If the receiver aperture is small enough compared with the beam 
diameter that it is always in the beam, then the beam wander will 
cause apparent intensity modulation. However, assuming the angle 
sensor to be insensitive to this modulation, the beam wander will 
not affect angle of arrival measurements. If the beam diameter, 
the aperture diameter and the rms beam wander are comparable in 
magnitude, then the beam wander may cause the beam to entirely 
miss the receiver part of the time. It must be supposed that this 
situation would degrade system performance to an intolerable ex- 
tent and it should, therefore, be avoided by suitable system de- 
sign. 


9 



RR-542 


2.2.3 Breathing 

Random fluctuations in beam cross-sectional area, called 
breathing, cause intensity modulation effects to be observed if 
only part of the beam is collected by the receiving aperture. 

Both Beckmann and Hodara give conditions on the desired signal 
modulation index, m, such that the signal modulation power exceeds 
the breathing modulation power. The derivations of these condi- 
tions involve many assumptions and approximations. Furthermore, 
they are not readily verifiable via experiment because those modu- 
lation effects caused solely by breathing cannot be separated from 
those due to other causes such as beam wander and boiling. In 
any case the results are as follows. 

According to Beckmann the condition which must be satisfied 
is (for the isotropic case) 



(B) 


Hodara gives two relations for L/L in different ranges. They 
are 


2 iii 


vl/4 


< m < 1 for LA- « 10 
— — c 


(HI) 


and 
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V2/< 


L‘ 


- 1/2 


< m < 1 for L/L » 10 
— — c 


(H2) 


Numerical evaluation of these three expressions for 1 < L/L < 10 

2-12 ~ c — 

and > = 10 , a value which agrees well with experimental 

measurements over paths near the earth's surface, shows that ex- 

3 

press ions B and HI give values of m < 1 for L/L < 2 x 10 and 

c “ 

that these values agree within a factor of ten for 
2 X 10^ < L/L < 2 X 10^. Expression H2 gives values of m < 1 

^ A 

for 1/L^ > 2 X 10’ . This particular expression is difficult to 
reconcile with physical intuition, since it implies that modula- 
tion effects due to breathing decrease with increasing pathlength 
when the path is very long. Moreover, for such large values of 
L/L^ the near field condition, L^ »-^iX is easily violated. 

If nothing else, these theoretical results indicate that 
breathing effects alone can produce deep fading in the received 
signal. Therefore, the laser receiver system should include an 
AGC which is wideband enough to minimize the relatively low fre- 
quency AM induced by the atmosphere. 


2.2.4 Phase Fluctuations 

Variations in the index of refraction along the transmission 
path produce random fluctuations in the velocity of propagation of 
the laser wavefront. If a coherent plane wavefront is assumed at 
the transmitter aperture, then the phase at any point in the re- 
ceiver aperture will vary randomly about a value which increases 
linearly with time. This may be termed temporal incoherence. In 
addition, the phase across the receiver aperture will be non-uni- 
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form at any instant of time. This will be referred to as spatial 
incoherence. In what follows only the isotropic case is considered. 

The first of these effects is found by Beckmann to produce 
a phase variance given by 


where 


k = ci)/c = 2tt /X 

CO = laser frequency 

X = laser wavelength 

This expression also yields the phase variance of a modulation 

signal at co if k is replaced by k 
m m 


where 



CO /c = 2ir/X 
xti xn 


Hodara's result is identical to Beckmann's except it does not con- 
tain the factor, The second effect may be characterized in 

terms of the cross -correlation function of the instantaneous phase 
at two points in the aperture plane separated by a distance, X. 
Bechmann's result in this case is 





2 
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and Hodara's result differs from this only in the omission of the 
factor, V^. 

Random phase variations in the received signal cause spuri 
ous frequency modulation. If we define 


Aco = 


d M 
dt 


then the variance of this frequency modulation may be found. The 
derivation involves the introduction of the inhomogenity correla- 
tion time, T^, and the results are given in terms of and v^ 

where v„ = L /T . Beckmann finds 
i c c 


<Aa)^> = 2\lTr <An^> Vrp L/L 

i c 

Hodara's result is the same except the factor 2'^ is omitted. 


3. EXPERIMENTAL RESULTS 

The theoretical results presented in the previous section 

express the mean square values of the various effects of interest 

in terms of pathlength, L, wavelength X, and the parameters of 

the correlation functions of the random component of the index of 

2 

refraction, <^ >, L and T . A number of investigators have per- 

c c 

formed experiments designed to measure atmospheric effects on a 
laser beam. In this section those results which are most relevant 
to laser tracking and ranging systems will be presented and dis- 
cussed. 
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3.1 FUNDAMENTAL ATMOSPHERIC PARAMETERS 

As pointed out in Section 2.1 atmospheric turbulence effects 

on laser beam propagation can be characterized in terms of the 

2 

variance of the index of refraction, <An >, and a correlation dis- 

2 

tance L^. Hodara gives empirical formulas for these parameters 
as functions of altitude, h in meters. These are presented below. 

<An^> = 10 exp( -h/1600) 

= 0.4h/(l + 10 ‘2h) 

2 

The first of these indicates that > is maximum at the surface 

(h=0) and decreases rapidly with altitude. L on the other hand 

c 

is minimum at the surface and increases toward an assymptotic value 
of 40 m with increasing altitude. Hodara does not state how these 
relations were obtained nor the conditions under which they are 
valid. However, they will be assumed for purposes of the present 
discussion to give results which are at least approximately cor- 
rect. 


3.2 CORRELATION FUNCTIONS AND SPECTRAL DENSITY FUNCTIONS 

The design of computer programs for smoothing and estimation 
of laser tracker data would be materially benefitted if we could 
develop a statistical model of the measurement noise. The theo- 
retical results of Beckmann and Hodara can be used to predict the 
noise variance, but the variance alone is not a complete descrip- 
tion of the noise process. What is required is a correlation 
function or a spectral density function. There is some theoreti- 
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4 

cal and much empirical evidence indicating that an exponential 
cosine correlation function is a suitable choice for the modelling 
of atmospheric turbulence effects. Such a correlation function 
has the form 


R( r) = A e ^ cos c r 
The corresponding spectral density function is 




2 1,2 2\ 

2 A K 


0 ) + \k +c / 

7T 

4 

Cl) + 

„/, 2 2] 2 L2 2] 2 

2 \ k -c / Cl) + \k +c ) 


2 2 

If 3c > k , W( co) has a maximum at 


/ 2 L2 2)^^^ 

(jt) = \k + c / [2c - (k + c / 


1/2 


In cases where the spectral density decreases monotonically with 
frequency from a maximum value at o) = 0, these models may be sim 
plified by taking c = 0. In this case 

R(t) = A 

and 


W(o)) 


2 A k 1 

2 , 2 
IT Cl) + k 
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The parameters A and k of these functions have the following in- 
terpretation. Let the random process they describe be |x(t)} . 
Then by definition 

~ x( t) x( t+r) 

from which it follows that 

R(0) = — — - 5 - = var x = A 

x(t)2 

Thus, the parameter A is the mean square value of the random pro- 
cess. The parameter k may be viewed as the inverse of the corre- 
lation time constant. The parameter c is the angular frequency 
of the periodic component, if any, of R( r) . 


3.3 QUIVERING 

5 

Kurtz and Hayes have performed a series of experiments in- 
tended to measure random fluctuations in the angle of arrival of 
a laser beam. Two paths were used, one 3200 m long; the other 
165 m long. The receiver employed a Questar telescope with 9 cm 

aperture and a star tracker tube. The transmitter was a He - Ne 

o 

laser operating at 6328A. The receiver was mounted 2 m above 
ground. The elevation angle of the receiver was about 4° for both 
paths. The star tracker tube produced two outputs proportional 
to the focused spot horizontal and vertical displacement from a 
central reference position. These outputs were recorded and then 
analyzed to obtain amplitude spectra in the range 2 - 150 Hz with 
a 1 Hz filter bandwidth. After filtering to a 15 Hz bandwidth 
these outputs were also processed by a probability density analyzer. 
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The results of 55 measureinents over the 3200 m path and 52 
measurements over the 165 m path were processed. No significant 
difference was found between measurements taken in the vertical 
and horizontal directions. The computed characteristics of the 
two sets of data are presented in Table 3-1. 

TABLE 3-1 


ANGLE OF ARRIVAL MEASUREMENTS 

Path Length 

3200 m 165 m 


Mean rms 0-15 Hz (jLtrad) 

15 . 7 

9.2 

Limiting rms 0-15 Hz ( jurad) 

26.0 

17.0 

Limiting peak 0-15 Hz ( /irad) 

48.3 

28.6 

Mean rms 0-15 Hz (jLtrad) 

31.4 

18.4 


Typical frequency spectra derived from the recorded measure- 
ments are shown graphically by Kurtz and Hayes. While these vary 
somewhat in shape, they are all characterized by a maximum spectral 
density at low frequency (below 20 Hz) and decreasing density with 
frequency above 20 Hz. At 150 Hz (the upper limit of the data pre- 
sented) all of the spectral densities are at least 10^ of their 
maximxim values, which indicates a significant amount of noise power 
above this frequency. No attempt is made by the authors to compare 
their experimental spectral densities with theoretical results nor 
to fit them with analytic expressions. It has been found, however, 
that many of their measured densities can be approximated quite 
closely by a spectral density function of the form (see Section 3.2) 
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W(w) 


2 A k[ 1 

2 ,2 

7T + k 


Typical parameter values for the 3200 m case are 

A s 10^ ( firad) ^ 
k = 200 rad /sec 

Values which give a good fit to data for the 165 m case are 

A = 5 X 10^ ( jurad) ^ 
k = 24 rad/sec 

In Section 2.2.1 the theoretical variance for angular quiver- 
ing in the x- and y- directions was given as 

<A0^> = <A0^> = II <An^> L/L 

where 

' 2\jTr ( Beckmann) 

, M = ) 

1 ( Hodara) 

The variance of the beam wander is (see Section 2.2.2) 
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<Ax^> = a. L^/L 

c 

where 

■\/jr ( Beckmann) 

a - 

1 ( Hodara) 

2 

The theoretical formula for <A0 > = A predicts a linear dependence 

on L. This is borne out by the experimental results, since 

o _■] o 

A( 3200) s 20 A(165) . If we assume a value for <Ai > of 10 , we 

may compute using the experimental values of A. Using Beck- 
mann's forimila this yields 

L s 1 m 
c 

for both experimental cases. It should be noted here that this 
empirical value of does not agree well with Hodara 's predic- 
tions (see Section 3.1 and Table 3-2) for the two cases which are 
considerably larger. Of course, this discrepancy may be apparent 

rather than real, since there is no reason to believe the assumed 

2 

value of <An > is either accurate or constant independent of the 
particular path. 

It is clear that the experimental data of Kurtz and Hayes 
can be fitted very closely with spectral density functions which 
correspond to exponential or exponential cosine correlation func- 
tions and the predicted linear dependence of A on pathlength, L, 

is corroborated by the experimental evidence. However, too little 

2 

is known about the dependence of <Aa >, L and k on atmospheric 
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and equipment parameters to model the quivering correlation or 
spectral density functions quantitatively under general conditions. 

The best fit values of A and k were used to compute rms quiv- 
ering in the frequency range 0-15 Hz by ntraierical integration of 
W{ o)) . The result for the 3200 m case was approximately 20 /jtrad 
and for the 165 m case 8 /xrad. These values agree well with the 
experimental results which were 15.7 and 9.2 fivad respectively 
(see Table 3-1) . The theoretical rms quivering and beam wander 
for the two experimental pathlengths have been computed using the 
formulas derived by Beckmann and Hodara. The results are given 
in Table 3-2. 


TABLE 3-2 

THEORETICAL QUIVERING AND WANDER 

Average path height, 
h = 2 + L sin 4°/2 (m) 


L 3200 m 
115. 


165 m 
8.75 


Correlation distance 

L (h) = 0.4h/(l+0.01h)(m) 20.4 3.2 

c 


Index variance 

<M^> = 10 exp( -h/1600) 

Rms quivering, AQ (jurad) 

( Beckmann) 
(Hodara) 

Rms beam wander, Ax(cm) 

( Beckmann) 

( Hodara) 


0.93 X 10"^^ 


1.00 X 10“^^ 


22.8 

12.1 


5.16 

3.87 


13.5 

7.16 


0.16 

0.12 
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Comparing Tables 3-1 and 3-2 we see that for each pathlength 
the two theoretical values for rms quivering bracket the measured 
value in the band 0-15 Hz. Thus, predictions based on theory are 
of the same order of magnitude as the measurements. However, the 
theoretical values include contributions at all frequencies and 
the measured values even in a band limited to 150 Hz are greater 
than predicted. This discrepancy is not surprising, since the 
theory predicts quivering due only to random refraction of the en- 
tire beam, whereas the experimental receiver is sensitive to other 
effects and it introduces some noise itself. The star tracker 
tube measures the position of the beam image in the focal plane of 
the receiver optics. This position will change not only as a re- 
sult of deviations in the angle of arrival of the entire beam but 
also as a result of phase front disturbances within the beam 
(boiling). Furthermore, if the entire beam is collected by the 
receiver aperture, then the star tracker tube will be sensitive 
to beam wander too. It is not possible to say whether this last 
effect is significant in the experiments described here, because 
the characteristics of the laser transmitter are not given by 
Kurtz and Hayes in sufficient detail. What we need to know is 
the transmitter aperture or beam divergence. Assuming the laser 
to have a beam divergence of approximately 0,7 mrad (typical of 
many commerical units) then the half-power beamwidth at the re- 
ceiver aperture is given by 

D = 0.7 X 10"^ L 
r 

= 11.5 cm for L = 165 m 
= 224 cm for L = 3200 m . 
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It may be seen that under this assumption the received beam over 
the 165 m path is about equal in size to the receiver aperture 
(9 cm) . It is possible, therefore, that for this pathlength beam 
wander may contribute to the measured angle of arrival. Over the 
long path, however, the postulated received beamwidth is so large 
compared to the predicted beam wander and the receiver aperture 
that sensitivity to beam wander should be very small. 

To summarize this discussion it may be said that the experi- 
mental results obtained by Kurtz and Hayes show a measured rms 
quivering of the same order of magnitude as predicted by theory. 
However, in practice other random effects will contribute to the 
noise in angle of arrival measurements so that theoretical values 
based only on quivering effects may be expected to be somewhat 
low. No attempt was made in these angle of arrival experiments 
to correlate the measured spectral densities with atmospheric con- 
ditions such as pressure, temperature and crosswind velocity al- 
though some information on weather conditions during each experi- 
mental run is reported. More experimental work is needed to de- 
termine the dependence of the noise model parameters, A and k, on 
atmospheric conditions, pathlength, path elevation, motion of the 
path relative to the atmosphere and equipment aperture sizes. 


3.4 PHASE FLUCTUATIONS 

The results of a phase fluctuation measurement experiment 

6 

have been reported by Buck”, The apparatus used was an equal-arm 
Michelson interferometer with round-trip arm lengths of 48.8 m. 

The reference arm was enclosed in a tube. The mixed measurement 
and reference beams were detected at two points separated by a 
distance X/4 to obtain quadrature measurements of the random phase 
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fluctuation at the laser source frequency (X = 0.6328 ji) . Data 
were taken at night under calm conditions over a horizontal path 

1.4 m above ground. The rms phase fluctuation was found for about 
16 samples each 1.7 sec long. The average measured fluctuation was 

2.5 rad. 

The phase fluctuation predicted by theory has been calculated 
using Beckmann's formula 

with 


k = 27 t/X = 0.99 X 10*^ rad/m 
<6n > = 10 

L (h) = L (1.4) = 0.55 m 
L = 48.8 m 


The result is 

AA = 68.8 rad 
^rms 

The theoretical and experimental results do not agree well. The 

reason for this large discrepancy is not know. Buck points out 

that his experimental data are quite limited in number. Also, 

2 

the value of <h > used to predict A0 in the above calculation 

^rxas 
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was based on Hodara's empirical formula. If the experimental con- 

2 ”15 

ditions were so calm that > was as small as 10 , then the 

theoretical prediction of M would be 2.17 rad, which agrees 

"^rms 

well with the measured value. 

The only conclusion that may be drawn from this discussion 
is that insufficient experimental evidence is available to either 
confirm the theory or discredit it. In any case, since in the ap- 
plication of interest the target range is measured in terms of 
phase shift in several modulating frequencies, it is the atmo- 
spheric effects on the modulation phase rather than the laser car- 
rier phase that are of interest. If one considers the following 
worst case conditions 


m ma x 
X 

m min 
k 


m max 


= 30 mHz 
= 10 m 


= 27T/X 


m m in 


= 0.628 rad/m 




max 


max 


L max 



= 20 km (round trip) 
= 40 m 


then the theoretical phase fluctuation is 


Ad =0.75 mrad . 
^rms 
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In this example a change in phase of 2 jt rad indicates a change in 
target range of 5 m. The phase-to-range scale factor is thus 
5/27T m/rad. The theoretical worst case phase fluctuation corre- 
sponds, therefore, to an rms range error of only 0.06 cm. If the 
phase noise due to atmospheric turbulence actually encountered in 
practice is no worse than this prediction, it may be concluded 
that atmospheric effects on the laser range measuring system will 
he small compared to equipment induced noise. 

3.5 INTENSITY FLUCTUATIONS 

In a laser ranging and tracking system the fundami^ntal vari- 
ables which are measured are angle of arrival and modulation phase 
shift. To a large extent the measurement instruments are inher- 
ently or by design (through automatic gain control) insensitive 
to intensity fluctuations. However, such fluctuations are of in- 
terest in the design of the AGC system and in evaluating the sec- 
ond order effects which they will have on the primary measurements 
Some of the experimental observations of intensity variations will 
therefore, be discussed in this section. 

3*5.1 Atmospheric Modulation 

The extent to which the atmosphere causes modulation of the 

6 

intensity of a laser beam has been investigated by Buck” and by 

7 

Subraraanian and Coll ins on . Buck defines a fluctuation index, a 

n 

as the standard deviation of the received intensity signal nor- 
malized to its mean. 
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The fluctuation index was measured as a function of receiver aper- 
ture for pathlengths in the range 0.55 - 145 km. The index was 
found to decrease from a value of about 1 at an aperture of 1 mm 
to values less than 0.02 for apertures which were large compared 
to received beamwidth. No systematic dependence on pathlength 
was observed in these experiments. Buck attributes this modula- 
tion to boiling within the beam. However, even if the intensity 
distribution across the beam were smooth, some fluctuation would 
be expected due to beam wander. 

The results reported by Subramanian and Collinson are in 
terms of percent modulation as a function of receiver aperture. 
They alsa found a decrease in modulation with increasing aperture 
to the point where essentially all of the beam was collected. 
However, they observed a relatively large (around V^) and constant 
residual modulation for apertures larger than the received beam. 

These experiments indicate that beam boiling and wander 
cause large spatial and temporal fluctuations in the intensity of 
the received beam. If the receiving aperture is small compared 
to the extent of the beam, severe intensity modulation of the re- 
ceived signal results. Enlarging the receiving aperture reduces 
this effect, but only to the point where the aperture collects 
substantially all of the beam power. Decreasing the size of the 
transmitted beam by increasing the colliiBator aperture can reduce 
the atmospheric modulation of the received beam just as increasing 
the receiver aperture does. However, there is a limitation to 
this approach. Buck has found experimentally that a transmitting 
aperture of 11 cm yields a minimum beamwidth over paths a few 
kilometers in length. The minimum diameter of the beam is about 
four times the theoretical Airy disc size for this transmitting 
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aperture. This effect is attributed by Buck to spherical aberra- 
tion in the transmitter optics. Since the measurement technique 
involved the direct exposure of photographic plates for 30 sec- 
onds, the measured bearawidths Include the integrated beam wander. 
This probably accounts for most of the discrepancy between experi- 
mental and theoretical beamwidth. 

3.5.2 Spectral Density of Modulation 

Reports of experimental measurements of the spectral den- 
sity functions of atmospheric laser intensity modulation have been 

6 8 
given both by Buck and by Gilmartin and Horning . However, 

neither of these experiments appears to have been well controlled, 

and the results are presented in sufficiently different forms as 

to make comparison difficult. Generally, both investigations 

yielded power spectral densities which decreased with increasing 

frequency. The frequency range from 0.2 to 500 Hz was covered. 

Buck' s measurements were made in Colorado at an average 
elevation of 2 km above sea level. The transmitter was on a mesa 
60 m above the general terrain. Measurements were made over dis- 
tances of from 0.55 to 145 km. The average spectral slopes in 
the range 1 - 10 Hz and 10 - 20 Hz were compared as functions of 
path length and receiver aperture with a fixed transmitting aper- 
ture of 15 cm. No obvious correlation was found between these 
slopes and pathlength. Both slopes did, however, become more 
steeply negative with increasing receiver aperture. Apparently, 
no crosswind velocity data were taken. 

The experiments performed by Gilmartin and Horning used a 
folded optical path at a desert location. The transmitter aper- 
ture was 5 cm and the receiver aperture was 32 cm. In one setup 


27 



RR-542 


both receiver and transmitter were on a ridge 247 m above the 
desert floor and an array of from 12 to 25 corner cube reflectors 
each 5 cm in diameter on the desert was used to return the beam. 
One way path lengths of from 3.2 to 25.6 km were used. In the 
other setup both ends of the range were on the desert floor sep- 
arated by distances of from 0.4 to 4.15 km. Measured power spec- 
tral densities W(f) were plotted in the form fW(f) vs In f. 1110 

function, fW(f), has a maximvim value at some frequency, f . The 

3 ^ 

formula for this frequency as derived by Tatarski is 

f = K Wrjr/yJiA 

Jr 

where 

v^ = crosswind component of velocity 

L = path length 
X = wave length 

Gilmartin and Horning found good correlation between their mea- 
surements and this predicted value of f^. Their experimental 
value for K was 2.5 + 32^ for the ridge-floor paths and 2.8 + 34^ 
for the floor-floor paths. They found further that all of their 
measured spectral densities could be closely fitted with a log- 
Maxwellian curve, viz 


fW(f) 


^J^f (In / 3 )^ 


exp 


(In 


The parameter j3 was found experimentally to have a value of about 

12 . 
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Buck has not taken crosswind velocity into account and he 
reports no significant variation with path length in measured 
spectral densities. This finding can be reconciled with Gilmartin 
and Horning's only if one supposes that in Buck's experiments the 
crosswind velocity increased fortuitously as the square-root of 
path length. Buck did find, however, that for a given transmitter 
aperture and path length the receiver aperture size affected the 
spectral density shape. Gilmartin and Horning on the other hand 
used fixed transmitting and receiving apertures and a variable 
(but not precisely specified) retroref lector aperture. Thus, it 
would appear that their results may depend on the aperture effect 
found by Buck as well as the crosswind effect to which they at- 
tribute them. It must be concluded that neither group of investi- 
gators has controlled or reported all of the parameters which a 
comparison of their experiments shows may be of significance. 
Clearly more carefully designed experiments need to be performed, 
if these ambiguities are to be resolved. 


4. CONCLUSIONS 

The purpose of this section will be to draw from the preceed- 
ing discussion of theoretical and experimental results those con- 
clusions which are of relevance to the problem of designing com- 
puter programs for use in conjunction with the planned laser 
tracking and ranging system. The raw data input to the computer 
from the tracker consists of two orthogonal angular pointing er- 
ror measurements and a range measurement. What is desired for 
computer program design purposes is a model of the statistical 
characteristics of the noise in these measurements caused by ran- 
dom turbulence in the atmosphere. 
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The primary causes of noise In the angle of arrival of the 
laser beam at the receiver are quivering and boiling. As was 
pointed out earlier, beam wander and Intensity modulation can also 
contribute to this noise depending on the size of the receiver 
aperture In relation to received beamwldth and on the characteris- 
tics of the error sensing device and associated AGC. 

Theoretical predictions and experimental measurements of angu- 
lar quivering have been compared In Section 3.3, It was found 
that the experimentally measured spectral density functions for 
3200 and 165 m fixed pathlengths could be fitted In many cases 
with a model having an exponential correlation function of the 
form 


R(r) = A 

where 

A “ 3L ( jurad) ^ 

L = pathlength (m) 
k “ 200 (rad/sec) for L - 3200 
s 24 (rad/sec) for L - 165 


Attempts to correlate this empirical result with theoretical pre- 
dictions were not very successful. The parameter A Is the vari- 
ance of the random quivering which Is given theoretically by 
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A = <A0^> = <A0^> 

X y 

= juL L A-^ 

where 

fi = 2^ ( Beckmann) 

= 1 (Hodara) 

2 

> = variance of index of refraction 

L = correlation distance of inhomogeneities 
c 

2 

According to Hodara both <<2n > and are functions of altitude, 
h. Thus, over a constant altitude path the thoery predicts a 
linear dependence of A on L. However, for the two paths in ques- 
tion h was not constant, since they were both elevated about 4° 

2 

above horizontal. If <An > and L are evaluated using Hodara' s 

c . 

formulas and an average value of h over the path, values of A are 
obtained which are much smaller than the best fit experimental 
values. This difference is due in part to the fact that the ex- 
perimental measurements include noise from other sources in addi- 

2 

tion to quivering. It seems quite probable, too, that both <Ani > 
and L are strongly dependent on atmospheric conditions which vary 
with both time and space. Representing them as functions of alti- 
tude only must be a simplification based on some average condition 
of the atmosphere. 
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If one considers the spectral density function 


2 A k 

1 1 

7T 

2 

Ct) 



the parameter A is the total area under the W(a>) - curve i. e., 

00 

/ W( co) dw = A 
o 

The parameter, k, determines the shape of the W( co) - curve. In 
particular* 

W(0) = 2A/7Tk 

and 

W( k) = A/i7k = W(0) /2 

Thus, at CO = k the spectral density function is one half of its 
zero “frequency value. In view of the theoretical model of atmo- 
spheric turbulence discussed in Section 2, the parameter k must 
depend on correlation distance, L^, wind velocity, v,j,, and tempera- 
ture and pressure gradients along the path. However, the nature 
of this dependence is unknown. More experimental work must be 
done, before it will be possible to predict k quantitatively. 
Qualitatively Kurtz and Hayes’ experiments show that k increases 
with pathlength and is widely variable with time for a fixed path. 
(Values in the range 125 - 375 rad/sec were obtained in fitting 
W( ct)) to different experimental spectra for the 3200 m path.) 


32 



RR-542 


The theoretical formulas for angular quivering may well be 
correct, but because of the large discrepancies between what they 
predict and the experimental results available at present, they 
must be considered unproven. The difficulty may simply be due to 
lack of precise knowledge of the variance and correlation distance 
of the random fluctuations in the atmospheric index of refraction. 

It is clear, however, that angular quivering effects are of suffi- 
cient magnitude to significantly affect the operation of the tracker 
pointing system and data processing algorithms. For purposes of 
designing computer programs for use with the tracker all that can 
be done pending more thorough experimentation is to assume the 
angle noise may be represented by an exponential correlation func- 
tion with fixed and hopefully conservative parameters. 

As pointed out in Section 3.4 the theoretical phase angle 
noise at the modulation frequencies used in the laser tracker 
ranging system is very small. Because the phase comparison cir- 
cuits used to measure phaseshift over the laser path are sensitive 
to amplitude noise, atmospheric intensity modulation effects can 
produce noise in the range measurements. Such intensity modula- 
tion noise may be minimized by careful design of the optical and 
electronic components of the tracker. The choice of transmitter, 
re tr ore flee tor and receiver aperture sizes is important (see 
Section 3.5.1). Probably even more important is the design of the 
AGC circuitry. This should be capable of minimizing the effect 
of fading which has been found to have a frequency spectrum from 
dc to several hundred Hertz. The capabilities of the AGC system 
also have bearing on the problem of distinguishing deep fades from 
loss-of-track. The automatic reacquisition system must be able 
to discriminate between these two conditions. 
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To summarize this discussion the following observations are 

made : 

1) With respect to angular position measurements, atmo- 
spheric quivering will introduce significant noise. 

While this effect cannot as yet be predicted quanti- 
tatively with precision, a method for approximating 
its correlation function has been presented based on 
the available theoretical and experimental results. 

2) So far as can be determined from the literature on 
the subject, range noise due to phase shift in the 
atmosphere will be negligibly small at the planned 
modulation frequencies. 

3) Range noise is more likely to result from sensitiv- 
ity of the phase measuring circuitry to intensity 
modulation of the laser beam by the atmosphere. The 
extent of this effect depends on the parameters of 
many systems such as the AGG circuit and the various 
pre- and post-detection filters used in generating 
phase error signals. The fact that the range mea- 
surement is derived from three different phase mea- 
surements each of which is subject to noise further 
complicates any noise analysis. To develop a quan- 
titative model of this noise would require detailed 
information on the characteristics of the range mea- 
suring circuitry. This information is not available 
for use in connection with this program. The best 
that can be done is to design the range data pro- 
cessing algorithms on the basis of reasonable as- 
sumptions concerning and approximations to the range 
noise model. 

4) As a final general comment, it is evident that a 
great deal more experimental work of a carefully 
controlled nature needs to be done, before the de- 
pendence of turbulence noise on the measurable state 
of the atmosphere and on laser system parameters can 
be established in an accurate, quantitative fashion. 

For purposes of the present project the information contained 
in this report will be utilized as fully as possible in designing 
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the computer programs for control, smoothing and estimation. 
Actual field testing of the resulting laser tracking system may 
well reveal a need for modification or refinement of certain sys- 
tem parameters to optimize noise rejection. Since the data pro- 
cessing is being done digitally such changes in parameters can be 
made very easily at the software level. 
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